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PHYSICS 



CONNECTING EXPERIMENT WITH THEORY: A MODEL-INDEPENDENT 
PARAMETERIZATION OF NEUTRINO OSCILLATIONS 

DORIS JEANNE WAGNER 

Dissertation under the direction of Professor Thomas J. Wcilcr 
Many experiments are currently looking for evidence of neutrino mass in the form of neutrino os- 
cillations. Oscillation probabilities are non-linear functions of the neutrino mixing matrix elements, 
so most comparisons of data to theory are based on simplifying models of the mixing matrix. We 
begin this dissertation with a review of neutrino interactions and a few of the popular models describ- 
ing neutrino masses and mixing. Next we present our model-independent description of neutrino 
oscillations and derive the predictions of various models in terms of our new "box" parameterization. 
Finally we use our boxes to find mixing matrices consistent with existing neutrino data. As more 
definitive data becomes available, these solutions will probably need to be adjusted; when such a 
need arises, our box notation will provide a convenient method for finding new solutions. 
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CHAPTER I 



INTRODUCTION 



Neutrinos, they are very small. 

They have no charge and have no mass 
And do not interact at all. 
The earth is just a silly ball 

To them, through which they simply pass, 
Like dustmaids down a drafty hall 

Or photons through a sheet of glass. 

They snub the most exquisite gas, 
Ignore the most substantial wall, 

Cold-shoulder steel and sounding brass, 
Insult the stallion in his stall, 

And, scorning barriers of class, 
Infiltrate you and me! Like tall 
And painless guillotines, they fall 

Down through our heads into the grass. 
At night, they enter at Nepal 

And pierce the lover and his lass 
From underneath the bed - you call 

It wonderful; I call it crass. 

- COSMIC GALL by John Updike Q 

Updike's description of the neutrino, while quite poetic, could use a bit of updating to be sci- 
entifically accurate. Neutrinos do travel great distances through matter before interacting, but to 
say they "do not interact at all" is incorrect. Bernstein, the author of reference [Q, also considers 



iFrom COLLECTED POEMS (Alfred A. Knopf). Reprinted by permission; ©I960 John Updike. Originally in 
The New Yorker. All rights reserved. 



Updike's poem and suggests '"And do not interact a lot . . .' is better science but worse poetry." 
But it is the claim of the poem's previous line that neutrinos "have no mass" which is addressed by 
my research. Bernstein happens to agree here with the poet, claiming that "most physicists would 
be willing to give high betting odds that the mass is exactly zero." Bernstein and Updike, however, 
were writing in the 1960s, and a few decades can change the prevailing views in science. No concrete 
experimental evidence of non-zero mass has yet been found, but the existence of massive neutrinos 
is currently favored on theoretical and indirect experimental grounds. In view of the prevailing 
scientific opinion, I suggest the following replacement for the first three lines of Updike's poem: 

Neutrinos, they are very small. 

They have no charge and not much mass 
And barely interact at all. 

This dissertation will present some of the consequences of non-zero neutrino mass and will develop 
a mechanism to connect the theoretical models with the potential results of experiments. 



1.1 Some History 

The early 20th century was a time of great change for physics. Driving much of that change was 
the discovery and study of nuclear decays. Beta decay, in which an atom changes atomic number 
and emits an electron, provided some puzzling results. The decrease in the atom's energy should 
be equal to the energy carried away by the electron. But not only was the energy of the emitted 
electron in a particular reaction usually less than the expected amount, it was different for each 
occurrence of the same reaction. Energy seemed to be randomly vanishing in beta decay! In 1929, 
Pauli suggested the existence of an undetected particle which carried away the missing energy 

In retrospect, Pauli's solution may seem obvious. But physicists in the 1920s were not yet 
in the habit of postulating the existence of new particles to explain surprising results. Only the 
proton, the electron, and the photon had been discovered. The prejudice against new particles was 
so strong, Bohr preferred to discard energy conservation. He suggested that individual reactions 
need not conserve energy as long as energy was conserved on the average. Pauli hesitantly proposed 



_> 



his "neutron" as a possible solution in 1927, and Fermi renamed it the neutrino when Chadwick 
discovered what we now call the neutron [|| . 

The existence of a chargeless, weakly-interacting fermion was generally accepted by the mid- 
19308, but experimental confirmation of the neutrino's existence did not arrive until the 1950s |], 
||. Reines eventually shared the 1995 Nobel Prize in physics for this innovative experiment. Yet 
these elusive particles have continued continued to raise questions in the decades following their 
detection. Some of the questions which intrigue researchers today, such as whether neutrinos and 
antineutrinos are distinct particles, have reemerged many years after they appeared resolved. Other 
questions, such as the number of neutrino flavors in nature, and whether the neutrinos are massless, 
have never really gone away. 

According to Pauli's proposal, the neutrino and the electron share the energy lost by an atom in 
beta decay. The minimum energy given to the electron is proportional to its mass, and the minimum 
energy carried off by the neutrino should yield a measurement of the neutrino mass. But in some of 
the decays, the electron seemed to carry off all of the energy, leaving none for the neutrino. Thus 
Perrin commented in 1933 "that the neutrino has zero intrinsic mass, like the photon" §, g. Yet 
experiments cannot demonstrate that the mass of a particle (even the photon) is precisely zero. 
They may merely place upper bounds on the mass. The bounds on the mass of Pauli's "neutron" 
(known now as the electron antineutrino) from beta decay currently require this mass to be more 
than 50000 times smaller than the mass of an electron, or less than 10 eV. Fifty thousand times 
smaller than the electron mass may sound negligibly small, but the existence of even a tiny non-zero 
neutrino mass has profound effects. Nevertheless, a prevailing belief in massless neutrinos continued 
well into the 1960s, when the Standard Model of particle physics was developed. This enormously 
successful theory describing particle interactions requires neutrino masses to be exactly zero. 

Starting in 1968 with the Homestake solar neutrino experiment Q, evidence has been building 
in favor of non-zero neutrino masses. The Homestake experiment counts the neutrinos coming 
from the sun, and the number it detects is far smaller than the number predicted by solar theory. 
This "solar neutrino deficit" has since been confirmed by four other experiments 0. In addition, 
experiments measuring neutrinos produced in our atmosphere find an "anomaly" in the ratio of 
muon-type to electron-type neutrinos ^ . Both of these results could be explained by neutrino mass 
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through its quantum-mechanical implications. As for the Standard Model, it is not the "Theory of 
Everything" physicists view as their Holy Grail. Adding masses to neutrinos is a natural extension 
to the theory, requiring merely additional particle fields, without significant formal changes. These 
minimal changes to the theory, however, have profound physical results. As will be discussed below, 
non-zero neutrino mass may lead to mixing between different flavors of neutrinos and to oscillations 
of neutrino flavor. The mixings and oscillations resulting from neutrino mass could explain the 
currently puzzling solar neutrino deficit and atmospheric anomaly. Large-scale neutrino experiments 
hope to address these issues, and my work provides a model-independent path between the data 
from such experiments and the answers to the mysteries surrounding neutrinos. 



1.2 Whence We Start: A Review of the Standard Model 

The Standard Model describes the interactions of quarks, leptons, gauge fields, and a Higgs 
scalar. Quarks come in three generations of two quarks each: up and down, strange and charm, and 
top and bottom. Leptons too come in three generations: electron, muon, and tau. Each generation 
contains a negatively charged lepton and its associated neutral neutrino. Each quark and lepton has 
an associated antiparticle as well. The gauge fields include the photon, carrier of the electromagnetic 
force, the three vector bosons (W + ,W~ , Z°), carriers of the weak nuclear force, and eight colored 
gluons, carriers of the strong nuclear force. Of these particles, only the Higgs scalar and the tau 
neutrino have yet to be detected by experiments. 

Experiments show that only the left-handed particles and right-handed antiparticles are affected 



by the charged weak nuclear force (see Appendix B.3.1 for a discussion on handedness, or helicity). 
This "maximal parity violation" by the charged weak interaction was discovered in 1957 by Wu, 
et. al. ||, [fio|| , and it necessitates an asymmetry between right-handed and left-handed particles 
in the Standard Model. Left-handed particles, which participate in charged weak interactions, are 
placed in "isospin doublets." Right-handed particles, which do not experience the charged weak 
force, compose "isospin singlets." 

The Lagrangian of a theory provides information on particle interactions. Each term in the 
Lagrangian describes either a kinetic energy, a mass, or an interaction. Fcrmion terms contain 
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at least one field with Dirac conjugation -0 (defined in Appendix ^J), and one without, ip' . The 
conjugated fields represent outgoing particles or incoming antiparticles; the unconjugated fields 
represent incoming particles or outgoing antiparticles. Each term in the Lagrangian must follow 
the known conservation laws. For example, electric charge seems to be conserved in nature, so the 
total incoming charge of any term in the Lagrangian should equal the outgoing charge. Isospin and 
its component in a chosen direction (I3) must be conserved, so the field with the positive I3 in an 
incoming isospin doublet (isospin i) must couple to the field with the positive ^3 in an outgoing 
isospin doublet. Fields in an incoming triplet (isospin 1) can couple to fields of another triplet or 
to two doublet fields. Singlets have isospin and thus can couple to any term without affecting the 
isospin of that term. 

Many of the sources listed in the bibliography of this paper contain reviews of the Standard 
Model. References JlO| and treat the Model in great detail, while references p3| , and 
consider only the aspects of the Model relevant to neutrinos. We will synthesize these treatments in 
the summary which follows. 

Neutrinos do not interact by either the strong nuclear force or the electromagnetic force, so we 
will ignore the terms in the Lagrangian corresponding to those forces. The terms of the Standard 
Model Lagrangian that are of interest to neutrino physics may be represented as follows: 

a a a a a 

% ( w+ _i_ w- <rrt 

ls w 



~— Z ^NC ~ W+J£ c + W~J^ C , where 



a 

1 / I 



T^VaLl^VaL ~ ^a L j^a L + sin 2 9 W (a L j^a L + a R ^a R ) 



(1) 



^v L ^v L - I ~ - sin 2 9 W ) l L ^l L + sin 2 6 w Ir^Ir, and 



^cc = y^^aLl^^L = VL'fh- 



Here oil is the field operator spinor (expressed in full in Appendix B.2) for the charged lepton for the 



neutrino spinor v aL ^ and vl and are the vectors containing all v aL and a^, respectively. The sub- 



scripts r and l indicated left- and right-handed particles, respectively (defined in Appendix B.3.1) 



Z M , W+, and are the fields of the vector bosons, e is the dimensionless charge of the electron, 
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which equals the square root of the fine structure constant, and 6\y is the Weinberg angle: 



2 



cos^<V = -f. (2) 



Appendices [A| and B.l contain a list of conventions used in this paper and descriptions of the 
gamma- matrices. 

The Standard Model Lagrangian contains mass terms for the charged leptons, but not for the 
neutrinos. Mass terms occur when right-handed particles couple to their left-handed partners, giving 
rise to terms of the form 

V'LTO^Vfl, (3) 

where is the mass of the fermion field ip = ipi, + ipR. In such a term, the outgoing left-handed 
field has isospin ^, and the incoming right-handed term has isospin 0. We know most fermions have 
a mass, so any physical theory must include mass terms for electrons, quarks, etc., but these terms 
seem to violate isospin conservation. Higgs JL5| , drawing from the work of Goldstone, Guralnik, 
Hagen, Kibble, Lange, and many others Jig , finally came up with a solution to this apparent 
paradox. Higgs' work was incorporated into the Standard Model by Weinberg. 

In the Standard Model's "Higgs mechanism," particles obtain masses by interacting with a Higgs 
field, introduced to conserve isospin: 

= X! 9%pipaL<t>D/3R + h.c. (4) 

a,/3 

4>d is the Higgs doublet, responsible for giving charged leptons and down- type quarks their masses: 



(5) 




ip a L is the left-handed lepton doublet of flavor a: 

( 



(6) 



G 



/3r is the right-handed charged lepton of flavor (3, and the g^g are coupling constants. Up-type 
quarks obtain masses in the Standard Model by coupling with the conjugate Higgs doublet |Q 



/ 



(7) 



where 02 is the second Pauli matrix as given in Appendix B.l 



The neutral Higgs field has all the quantum numbers of the vacuum, so it may obtain a "vacuum 
expectation value," or VEV: 

($>> = W2>|0>. (8) 

The charged Higgs field cannot obtain a VEV since the vacuum is electrically neutral, so the ground 
state of the Higgs doublet is 







(9) 



The Higgs field is dynamical, so we include a term representing the deviation from the vacuum value 
and arrive at our final form for the Higgs doublet: 







\ 



Kiting/ 



(10) 



H°(x) is the massive Higgs scalar, and the VEV (4>%) gives mass to particles through interactions 
such as 

9^ol(<I>d)Pr = aTM a/3 f3 R . (11) 

Because these terms involving the VEV do not conserve isospin, but the original interactions in 
equation (^) do, the Higgs mechanism spontaneously breaks the SU(2) symmetry in the ground state. 
The extra fields present in the original Higgs doublet (^) do not simply vanish in the symmetry 
breaking; they are transformed into extra spin states of the vector bosons. Such a discussion is 
irrelevant to the focus of this paper, but the interested reader can consult the sources |l7j or Q for 
more information. 
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In the Standard Model, the charged lepton mass matrix is diagonal, M a p = m a 5 a (3, and mass 
terms for neutrinos do not arise because the Model has no right-handed singlet neutrino fields to 
participate in the coupling of equation (|4|) . We will examine the consequences of the addition of a 
neutrino mass term in Chapter ||. 

The Lagrangian in equation (frl) does not explicitly display terms involving the right-handed an- 



tineutrino fields V R = CFl T , where C is the charge-conjugation operator, discussed in Appendix B.4 . 
These fields are CP-conjugates of the neutrino fields The Lagrangian does, however, contain 
antiparticle fields implicitly. C U SM above is CP-invariant, so if you replace the particles in equation 
(|l|) with their CP-conjugates, you get the same terms as you started with. Take, for example, a 
kinetic term for antincutrinos "^r^^O^r'- 



t^^d^R = Cvl T ^d ll GT7Z T = ( 7 V) C^^Cd^. (12) 



The term in equation (12) is a number, so it is equal to its transpose: 



( 7 V) T c^V^c^ = -(^pz^Vc-SV, (13) 

where we have added the minus sign because we switched the order of two fermion fields Vi and 
vl^^ and we have used the relations (421) from Appendix B.4. The operator C anticommutes with 



7°, and CV^C- 1 = -7", so 

Vjfd„V R = -(d^Wh^L. (14) 

The last step is to bring the partial derivative <9 M to the right of the spinor T7Z . We may do this by 
remembering that Cg M is a Lagrangian density and thus always appears in an integral. We may 



The term antineutrino can be a bit confusing, since the charge-conjugate neutr ino field is different from the 



13,1 



physical CPT-conjugate antineutrino field. These subtleties are discussed in Appendix 

3 When taking the transpose of a product of fermion fields, it is necessary to account for the anticommutation of 
those fields, so (ipitl)2) T = —'4 ) \ V '4>\ if ^1 an d ip2 are fermions. When taking the hermitian conjugate, however, no 
minus sign appears: (ipi^)^ = V4"^I- The hermitian conjugate of a product is defined by the previous expression, 
and the anticommutator is never taken. Consider the number operator N = a'a from quantum mechanics, which 
must be hermitian since it is observable: JVt = N, so (a T a)t must equal a T a, and not — a T a. 



8 



therefore use the technique of integrating by parts to obtain 



- {dpVZWvL = -(dtPZ-fvL) + Wl"d^ L , (15) 

where we have moved the partial derivative through the coordinate-independent 7 M in the second 
term. The first term on the right-hand side of equation (|l^) becomes a surface term and vanishes 
for integrals over all space-time. So we are left with 

V^^B^ R = -{duVftfvi* = VL^d^VL, (16) 

which shows the kinetic term is invariant under CP-conjugation. 

The CP-invariance of the weak current terms in C U SM may be shown by a similar treatment, 
realizing that under a parity transformation, the vector field becomes its negative, so under CP, 
Z — > —Z, and — > —W T . The neutral current term, like the kinetic term, becomes itself. The 
W + term becomes the W~ term, and vice versa: 

(w+wfh + w-TE<fv L 7 = (w+jt^riR + {wffQff, a 

= -w- (-v L T &) ^cif - w+ (-iltf) rcw T 
= -w-T L c T ^ T c^ T v L - w+vzc T r T c^ T i L 

= W-JE'fvi. + W+WflL. (17) 



where we have used equation fl426|) from Appendix to substitute for P R , we have obtained a 



minus sign by switching the order of fermion fields in taking the transpose as before, and we have 



again used the properties of C given in equation (421) in the Appendix. The charged lepton mass 
terms also swap under CP, 

dt R m a dt L + b. L m a bt R = —et^C' m a C 'a r T — a^C^m a CaL T 

= orC C' T m a aL +aZCF C^ T m a aR = aR~m a aL + OLm a aR, (18) 
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Lagrangian for electron-neutrino scattering is responsible for matter-induced oscillations for massive 



neutrinos, as will be discussed in Section 3.3. 



1.3 Whither We Go: The Changes Wrought By Neutrino Oscillations 

The Standard Model contains only left-handed isospin | neutrinos and one Higgs doublet, so it 
allows no neutrino mass terms. As mentioned above, the neutrino mass was thought to be exactly 
zero in the 1960s, when Updike wrote his poem, and the Standard Model was developed during 
this same period. The Model, like the poem, reflects the prevailing views of its time and therefore 
yields a neutrino mass of zero. In the last thirty years, however, both the solar neutrino deficit and 
the atmospheric neutrino anomaly have surfaced. Experiments counting neutrinos from the sun see 
far fewer electron neutrinos than theory predicts, and experiments measuring neutrinos created in 
the atmosphere measure a muon-type-to-electron-type-neutrino ratio much smaller than expected. 
These surprising experimental results may be explained by making neutrino masses unequal and 
therefore not all zero. Non-zero neutrino mass is now preferred for symmetry reasons as well. Every 
fcrmion but the neutrino in the Standard Model has a right-handed singlet state and thus a mass, so 
the absence of right-handed neutrinos in the Standard Model seems unnatural and rather arbitrary 
to those who seek consistency in nature. Most extensions of the Standard Model, including Grand 
Unified Theories not based on an SU(5) symmetry, give the neutrinos mass. Some of the mechanisms 
by which neutrino mass may appear are discussed in Chapter The existence of mass terms in 
the neutrino Lagrangian may lead to neutrino oscillations between flavor states, as presented in 



Chapter III . Chapter IV develops a parameterization to describe neutrino oscillations which is 
independent of the model used to generate neutrino mass. This parameterization is then used to 
describe the predictions of many popular neutrino mass models (Chapter 0) The predictions are 



compared to the existing neutrino oscillation data in Chapter VI. Chapter |VII| examines some of 
the implications of neutrino mass and oscillations and addresses the exciting future for this field of 
study. 



12 



CHAPTER II 



GIVING NEUTRINOS MASS 



The Lagrangian in equation (Q) does not contain a neutrino mass term. Masses of fermions arise 
in the Standard Model through the coupling of the fcrmion with a Higgs scalar doublet, as described 
in Section 1.2. The Higgs doublet interacting with fermions such as the electron has isospin h, 



and it couples an isospin ^ state with an isospin singlet. Such interactions flip the helicity of the 
fcrmion involved, so both a right- and a left-handed state of each massive fermion are necessary in 
the Standard Model. This type of interaction is described by a Dirac mass term. Neutrinos may 
acquire a Dirac mass simply by adding a right-handed neutrino (and a left-handed antineutrino) 
to the Standard Model. Right-handed neutrinos do not participate in the charged weak interaction 
and are therefore dubbed "sterile". In this work, sterile neutrinos will be denoted by Nr as in Jl9| |, 
while the active neutrinos will be represented by the traditional vl- The associated helicity states 
of the CP-conjugate antineutrino will be given by v R and iV^.F] 

Because neutrinos are neutral, they may obtain mass by an additional method. Unlike charged 
particles, whose particle and antiparticle states differ by the sign of the charge, neutrinos may be 
their own antiparticle, like the photon is. This type of fermion is known as a Majorana particle, 
after the physicist who first examined the possibility and ramifications of a self-conjugate fermion 
p0[ . If neutrinos are Majorana particles, they could form a mass term from the left-handed vl and 
the right-handed v R . Such a combination has isospin 1, so Majorana mass terms must be generated 
by a triplet Higgs interaction, which would be an addition to the Standard Model. 



2.1 The Differences Between Majorana and Dirac Neutrinos 

Two different types of neutrino states have been experimentally observed: vl and its CP- 
conjugate v R . The Dirac mass term introduces the right-handed neutrino Nr and its CP-conjugate 

1 I prefer Nr rather than the more apparent ur to represent the right-handed component for two reasons: it clearly 
distinguishes active from sterile neutrinos, and it avoids confusion between "Ur and (u)r = FjJ. In my notation, the 
former is denoted Nr, which is clearly different from the latter. 
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N L , and has the form 



Coirac = - (WM D N R + NrM^V^ . (22) 



As discussed in Section 1.2, terms with the antiparticles V R and N L are redundant and therefore not 
explicitly written. A Lorentz boost or spin flip will transform i>l to Nr and v R to N L , as shown in 
Figure |^a, which is adapted from reference pl| . CP conjugation, however, transforms to v R and 



Nr to N L , as discussed in Appendix B.4, so four distinct neutrino states are necessary for Dirac 
mass terms. 



Boost t> j Boost „ 

v L < ► Nr v l < ► V R 



J CP ICP \QV J CP 

Boost „ Boost 

P R < — > 7Y L ^ < — > 

(a) Dirac (b) Majorana 



Figure 4: A pictorial representation of the relationships between neutrino states for a) Dirac mass 
terms, and b) Majorana mass terms. 



Majorana neutrinos, however, are self-conjugate, so the sterile states are not necessary 



ty(maj)( X ) = C "fOlp* {m aj)( x ) = e * ^P(maj)(x), (23) 



where C is the charge-conjugation operator. Acting on a Majorana field with CP only changes the 
phase of the field and takes one helicity state into the other, so it has the same effect as boosting the 
state, as illustrated in Figure [|b. Thus only two distinct neutrino states are necessary for Majorana 
mass terms. 

Majorana mass terms have the general form 

Jr'Mm^l + JlM\^ R . (24) 



Because the Majorana field is self-conjugate, its mass matrix Mm is symmetric and so may be 
diagonalized by a transformation involving only one unitary matrix, as we will show. Following the 
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treatment found in reference ]l4J , we set the Lagrangian in equation ( g4j ) equal to its transpose, 
which is allowed since the Lagrangian is a scalar: 



IJrMm^l + h.c. = (lp R M M ^ L + h.c) T = - (^lMj f ^ T ) + h.c. 



(25) 



where the minus sign comes from swapping the order of fermion fields. Using the definitions of Dirac 
conjugates found in Appendix [A], the Majorana condition (|2"3|), and the expression (425) for ipL,R 



derived in Appendix B.4, we find 



b R M M ^L + h.c. = ~ (ylMli (~^ T L tf) T ) +h.c. = ^lMl I C^ T ^ L + h.c. (26) 
= ^K-C^Ml^L + h.c.^^M^L + h.c. (27) 



For the first expression to equal the last, Mm must be symmetric. 



Still following the treatment of source 1 14 , we remember that an arbitrary square matrix M is 
diagonalized by a biunitary rotation: 



D = V ] MU, or 
M = VDU\ 



(28) 
(29) 



with D diagonal and V and U unitary. For a symmetric matrix, Mm = M M , so 



M M Ml = M T M M^. 



(30) 



Using equation (29), we find 



VDU* (VDU^ = (VDU^f (VDU^ T , so 
VT>t/ t [/L> t V t = UDV 1 = U^ T D T V T V^ T D^ T U T = J7 tT D 2 U T '. 



(31) 
(32) 
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Thus, multiplying by U T on the left and V on the right, 

U T VD 2 = D 2 U T V. (33) 

Since the matrix U T V commutes with the matrix D 2 , they must be simultaneously diagonalizable. 
D 2 is already diagonalized, so U T V must be diagonal too. A diagonal unitary matrix has the form 
of a diagonal matrix of phases, X, so 

U T V = X, with X ij =5 ij e iii . (34) 

We may define a new matrix V which incorporates the phases of X: 

V = VX\ (35) 

giving U T V =1. We may thus diagonalize Mm with the single matrix U and its transpose V't = U T : 

I) VMy.il' XY'MmI. so (36) 

U T M M U = XD, (37) 

which is diagonal, as claimed above. 

2.2 The Dirac Mass Term 

The Dirac mass term couples left-handed active neutrinos with their right-handed sterile partners 
through an interaction with the same Higgs doublet which gives mass to up-type quarks in the 
Standard Model. 

C-Dirac = - X! Qapi'otL^DNpR + h.C, (38) 

with il> a L and 4>d defined in equations (^) and (Q), respectively, and g^g the Higgs-doublet coupling 
constants for the neutrino sector. Once the symmetry is broken, the Lagrangian term in equation 
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(pq) contains a mass term of the form 



Dirac 



v a L M a pNp R + h.c. 



(39) 



The mass matrix Mp contains all of the coupling constants: 



9%{<t>V)- 



(40) 



This type of mass term is analogous to the quark mass terms and requires only the addition of 
sterile neutrinos to the Standard Model. Generating small neutrino masses by such a term would 
require the Yukawa coupling between the Higgs and the neutrinos to be many orders of magnitude 
smaller than the similar coupling between the Higgs and the quarks. While this difference of coupling 
constants may indeed exist, no ready explanation for the difference of scale has been suggested. This 
unexplained lack of symmetry is unattractive, so other mechanisms for generating neutrino masses 
have been proposed. 

2.3 The Majorana Triplet Mass Term 

If neutrinos possess a Majorana nature, no sterile neutrinos are necessary to produce a non- 
zero neutrino mass. The Majorana triplet mass term couples active left-handed neutrinos to their 
CP-conjugate right-handed partners as follows: 



a term. An SU(2) triplet operator is therefore needed to conserve isospin. Such an operator may 
be a single object such as a Higgs triplet or a combination of objects yielding a single, isospin 1 
"effective" operator. 




(41) 



Because both vl and V R have I3 of +5, the third component of isospin changes by one unit in such 
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2.3.1 The Triplet Higgs 
Gelmini and Roncadelli proposed a model using a triplet Higgs 4>t with a Lagrangian given 

by 



^Triplet — - 22 ga/3 ( fr R V a 
a,f3 ^ 

The neutrino mass matrix in this model is 



n 



h.c. (42) 



Mt Q/3 = -2y/2g%{i&), (43) 

and the charged leptons do not obtain any mass contribution from the triplet since they do not couple 
to the neutral Higgs triplet field. To conserve total lepton number before symmetry breaking, the 
Gelmini-Roncadelli model assigns a lepton charge of 2 to the Higgs triplet. When SU(2) symmetry 
is broken, lepton-number symmetry is broken as well, resulting in the formation of an additional 
Goldstone boson, the Majoron. This Majoron obtains a mass far smaller than the mass of the Z 
boson, so it provides new final states affecting the measured value of the Z-decay width. Such an 
effect contradicts experimental measurements of the width, so the Gelmini-Roncadelli triplet model 
with a lepton-number-conserving Lagrangian has been ruled out fl4|| , |]l3| . 

All is not lost, however, for the Higgs triplet model. If lepton number is not conserved by the 
mass terms, the masses of the physical Higgs particles produced have no upper bound [0. Thus 
the Higgs triplet remains a viable possibility for generation of Majorana neutrino masses without 
the addition of sterile neutrinos, but at the cost of introducing new Higgs particles and abandoning 
lepton number conservation entirely. We shall not concern oursclf further with this model, but the 



interested reader may consult sources 14 , p2|, [E3|, and JL3 



2.3.2 Radiative Mass Terms and the Zee Model 
Additional mechanisms for generating neutrino masses exist at the one-loop level. Masses gen- 
erated by such a mechanism are called "radiative masses." One of the best-known radiative models 
using additional Higgs fields was proposed by Zee ]23] (for discussions of the Zee model, see also 
p5[ , p6| |, Q, p7p , and Zee's model in its simplest form uses only two Higgs doublets, (j> a 
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and 4>b, and a charged singlet Higgs field h. Such a singlet arises naturally from the simplest SU(5) 
theories, and it interacts with lepton doublets and Higgs doublets, but not quark fields The 
scalar h may be given lepton number —2, so its interactions with the lepton fields 

Ch*l = r P (HlC^ l ) z»h + h.c. (44) 

conserves lepton number]^] Here ipaL is the left-handed lepton doublet of flavor a, containing both 
neutrinos and charged leptons. i and j represent the SU(2) index: — v ah , and ip^ L — oll 
6ij is of course the antisymmetric tensor, and C is the charge-conjugation operator. The coupling 
constant f a P must be antisymmetric because of Fermi statistics; interchanging two fermion spinors 
must change the overall sign of any fermion coupling, so 

(r aL C^ L ) eijh = -f a (V^CCl) tjih. (45) 

Obtaining the left-hand side from the right-hand side involves swapping the order of two fermion 
fields, which provides one sign change; swapping the indices of the antisymmetric tensor, thereby 
providing a second sign change; and swapping the order of the indices of the coupling f a/3 . The 
coupling must therefore be antisymmetric in its indices to arrive at the required overall sign change. 
The coupling between the charged scalar h and the doublets <f> n takes the form |Q 

Ch* = Mabdj^tfrlh + h.c. (46) 

The €ij couples the doublets to conserve I3. This interaction violates lepton number by —2. Most 
treatments of the Zee model [|l2|, pi[ |, plfl , p7| , treat only the simplest case, in which only one of 
the Higgs doublets <p a couples to leptons, and then only diagonally in isospin space: 



C^=J29a^Mu + h.c (47) 



2 Reference mm incorrectly states that individual lepton numbers L e , L M , and L T are conserved by this term. Only 
the total lepion number L = L e + L M + L T is conserved by the Zee interaction. Zee corrects this misstatement in his 
later work IBql. 
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approximate form 
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This alternative form has its virtues, but for the rest of this paper we will consider only the form 
(p9|), remembering that it is merely the simplest, rather than the only, possibility for a mass term 
generated according to the Zee model. 



2.4 The Majorana Singlet Mass Term 

The sterile neutrinos Nr, being isospin singlets, may obtain a Majorana mass without Higgs 
fields. Such a "bare mass" term has the form 

CSinglet = -\ (j^M s N R + N^M^N^ . (51) 

This term is explicitly SU(2)-invariant, so it may be included in the overall Lagrangian before 
symmetry breaking. Such a term may also be generated by interactions with a Higgs singlet, but 
the use of such an extra Higgs in not necessary. The singlet mass term does not generate any type 
of mass for the active neutrinos and thus is somewhat uninteresting on its own, but it is possible 
in models with right-handed neutrinos and so it must be included in any full treatment of neutrino 
mass possibilities. 



2.5 Putting Them All Together 

Dirac terms and Majorana terms are not mutually exclusive, and mass terms of all three types 
(p2|), (|4"l|), and ( |5l| ) may occur simultaneously. We may also add additional sterile neutrinos XvR 
which do not have active counterparts. These additional neutrinos may be Majorana or Dirac 
particles. In the most general case, we can have flavors of active neutrinos v aLl ur flavors of 
their right-handed sterile partners N@r, and n x flavors of additional sterile neutrinos Xv Most 
reasonable scenarios set til = tir, but one can imagine other scenarios, so we will not a priori set 
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them equal. The neutrino mass terms of (J22j) , (41), and ( |5l| ) may be combined with the mass terms 
involving the additional sterile neutrinos in the matrix equation 



r 



(Fl n l x l ) M v 



( , \ 

V R 

N R 

\XR J 



h.c. 



(52) 



M v is the n x n neutrino mass matrix, where n = til + tir + n x , v contains all ul active neutrino 
spinors, Nr contains the riR Dirac counterparts, and \ contains all n x sterile neutrino spinors. M v 
may be broken into block form: 



^ M T M D M vx X 



M v = 



Ml 



M s M 



N X 



M T VX Ml x 



(53) 



/ 



where Mj- is x n^, Mp is x Hr, et cetera. 

The singlet and triplet Majorana terms fall directly out of equations (|5^) and(^3|) when the other 
couplings are set to zero, but the Dirac term requires a bit of algebra. For pure Dirac neutrinos, the 
matrices Mr, Ms, and all M x are zero. We are left with the Lagrangian 



■'Dirac 



\ 



M D 


o o o 



Ml 



\ ( v \ 

V R 
Nr 



h.c 



V XR J 



[VZM D N R + N L MlP R ] + h.c. 



(54) 



Using the relations ( |425| ) and (42C) from Appendix B.4, we see that the two terms in the brackets 
are identical: 

W L Mlv R = -NlC^MlCW T = WC T M D C^ T N R = WMdNr. (55) 
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The Dirac mass Lagrangian becomes 



C-Dirac = -v L M D N R + h.c, (56) 



in agreement with the definition of a Dirac term ((22[). The extra \ in the definition of the Majorana 
terms ( |4l| ) and ( |5l| ) does not appear in Dirac terms because the particle and antiparticle Dirac mass 
terms sum when we use the general notation ([52]). This factor is also consistent with what one 
expects for the equation of motion resulting from variation of C with respect to the neutrino fields. 
For Dirac fields, ip and if>' are varied independently, whereas for Majorana fields the two are not 
independent and their variation consequently returns a factor of 2. 



2.6 The See-Saw Model 

A popular class of neutrino mass models which includes both Majorana and Dirac mass terms 
is that of seesaw models. Seesaw models have achieved their popularity because they naturally give 
rise to small neutrino masses with minimum additions to the Standard Model. In such models, 
active neutrinos compose Dirac mass terms with their sterile counterparts, and the sterile neutrinos 
also form a bare Majorana mass term. Equation (|5^) becomes 



£mass — — ^ ( V L N^) 



M D 



M L D M s 




(57) 



In see-saw models, the elements of Mo, which arise from coupling to the Higgs doublet and the 
standard SU(2)xU(l) symmetry breaking, are on the order of lepton masses, while the SU(2)xU(l)- 
invariant elements of Ms are on the order of a larger scale, perhaps the scale of grand unified theories. 
Thus we may expand to leading order in Mg 1 . Roughly following the treatment in |Q, we find that 
U v has the form 

1 A MdM^L 

D 



-m~ 11 <mIa 



(58) 
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where A and B diagonalize the "light" and "heavy" mass matrices: 



MnMg 1 MpA T ^ = D Hght , and B^M S B T ^ = D h , 



(59) 



The mass eigenstates are 



N r - 



At {v L - MnMfNj) 

fit (m s 1 Hi^ d u l + n l 



(60) 



These are Majorana states, so the right-handed component is also the charge-conjugate of the left- 
handed state: 



1 v m \ 


- 








c 

■I 


K N m j 


R. 




\ N m j 


L- 





A 2 



'R. D 1V± S 



(61) 



B T {Ms l M£p R + N R) 



The seesaw model is so- named because it naturally gives small masses to the states v m , which are 
composed mainly of active neutrinos, while giving large masses to the states N M , which are largely 
made up of sterile neutrinos. 

In the above treatment, Mt has been set to zero. The seesaw model still works for non-zero Mt, 
provided the elements of Mt are still much smaller than the elements of M$- Such an inclusion 
of Mt has negligible effects on the heavy mass states, and changes the light mass matrix from 



-M D Mg L Mjj to M T - M D Mg X M^ |K 



2.7 Mass States vs. Flavor States 

Any mass term, regardless of the origin of the term, has the form tJjlM^r, where ip is a n- 
dimensional vector and M is an n X n matrix, ip can contain the active v ai the sterile N a , the 
additional Xv, or all three. As is the case for quarks, the mass matrix M for neutrinos is not 
necessarily diagonal. Neutrino states of definite mass correspond to eigenvectors of the mass matrix, 
which are generally distinct from the states of definite flavor. Adding neutrino masses, the lepton 
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Lagrangian of equation ([]]) in matrix form is 



C = ^id^ + Ud^l-c 1 Z^(wYv L -{l-2s\n 2 e w )l Ll ^l L + 2&in 2 e w l Rl n R ) (62) 
- c 2 (W+VZ~fl L + W-l^u L ) + JEMiIr + Tr~m}i l + JZM y ^ R + lh.Mt^ L . 

The subscripts on the mass matrices v and I label a matrix by the type of particle for which it 
provides mass; they do not denote flavor indices. The lepton spinors in equation (^2|) correspond to 
states of definite flavor. States of definite mass (to be denoted by the superscript m ) are achieved 
by diagonalizing the mass matrix, thereby rotating from the flavor basis to the mass basis. The 
diagonal mass matrices, Di and -D„, are given by 

D^Ur^MiU^, and D v = U~^M v U Nr . (63) 



As discussed in Section 2.1, when tp represents a Majorana spinor, M v is symmetric and we may 



take Un r = U^. At this point, however, we will consider a three-generation case without Majorana 



terms. The following treatment is extended to allow higher generations in Section 4.3 . 

The rotation from flavor states to mass states leaves the kinetic and neutral current terms 
unchanged, but the charged current is not invariant under the transformation: 

+ Wi L U{-Hd^u lL u^i L + 'Qu lR u l - 1 id^u lR Uf R 1 i R 

- Cl Z^ (VZU VL U-^U-^U VL v L - (1 - 2 sin 2 6 W ) ~l L U lh U^U 1l UfH L 
+ 2 sin 2 ewhU^U^U^U^h) 

- c 2 {w+veu vl u-^Ui l ui-H l + w-hUi L ui- L x ru VL u-^ L ) + hu lL u h 1 M l u lR ur R 1 i R 
+ TrU 1r u^m}u 1l u^H l + vzu vl u-^m v u Nr u^ 

+ F£U£U lL i8,sfV£ + VRUrn U i* id ^ l R 

- Cl Z^ {WU£U VL -fvZ ~ (1 - 2 sin 2 6 W ) TZU^UtrfVE + 2 sin 2 OwIrU^U^Ir) 
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+ l^D„N% + N%Dlv'l l 

- Cl Z M (v^^L - (1 - 2 sin 2 6 W ) X^-fVS + 2 sin 2 9 w T^Ir) 

- C2 {w+1%u£u 1l t»12 + w-lfu^u.^T) 

+ ifDd'R + VRD\V£ + v?D u N% + N%Dtv£ . (64) 



The charged current term in equation (|64|) is not simultaneously diagonalizable in both the flavor 
basis and mass basis; upon rotation to the mass basis, it acquires a mixing matrix, V = U VL . 
For Dirac neutrinos, the lepton mixing matrix V has the same properties as the CKM matrix in the 
quark sector. If the Lagrangian does not contain a neutrino mass term, as in the Standard Model, 
or if the neutrino mass eigenvalues are degenerate, the neutrino states may undergo an arbitrary 
rotation, so we may set U UL = Ui L , which yields V A In such models, flavor states and mass states 
are indistinguishable, so we may assume a diagonal charged lepton mass matrix for the Standard 



Model, as we did in Section 1.2, with no loss of generality. 

Even when neutrino mass terms are present, the matrices U VL and Ui L are not unique. Multi- 
plying each by any unitary matrix R leaves 

V = (RU lL ) 1 (RU„ L ) = U} L U VL (65) 

invariant. In the quark sector, an implicit choice R — U upL is made so that quark mixing is 
constrained to the down- type quarks, with the up- type flavor states assumed identical to mass 
states. Mixing in the lepton sector is typically similarly restricted to neutrinos. Under this implicit 
assumption, the neutrino flavor states are connected to the neutrino mass states by the mixing 
matrix V which appears in the weak charged current term: 

n 

v aL = ^ v ai v *L, with (66) 
1=1 
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Jcc = ^T7"iaL + h.c. = WlV^^ljL + h.c. (67) 

a— flavors — 1 

The mixing matrix occurring in the charged current term has two mass indices, since it couples the 
neutrino mass states to the lepton mass states. Under the assumption that charged leptons have 
a diagonal mass matrix, the lepton mass state Ul is equivalent to the lepton flavor state l a L 1 and 
the mixing matrix may be given a mass index and a flavor index and identified with the neutrino 
rotation matrix U VL . In the following Chapters, we will assume the lepton matrix is diagonal and 
give the mixing matrix one mass index and one flavor index. 
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CHAPTER III 



OSCILLATIONS 



If the eigenvalues of the neutrino mass matrix are non-degenerate, then neutrinos may change 
flavors as they propagate. We will show that the flavor of a neutrino changes in a cyclic fashion, 
so the neutrino flavor oscillates with distance traveled. To see how these oscillations arise, let us 
consider the neutrino Lagrangian and derive the appropriate equations of motion. 

3.1 The Equations of Motion 

For simplicity, we will consider a Dirac mass term as in equation (|22]). The neutrino Lagrangian 
in vacuum is then given by 



a— flavors 



v a Lid^ v aL + N aR id^N aR - ^ (v aL M Da pNp R + N aR M\ )ap Vi3 1 



f3— flavors 



(68) 



Varying (p8|) with respect to v ah and then with respect to N aR yields the following equations of 
motion: 



-— = id^v aL - V M Da gN pR = 0, and (69) 

(3 

spy 

— = id^N aR -J2M j Daf3 u 0L =O. (70) 

oJy aR *y 

Solving for N aR in the first equation and then substituting into the second, we obtain an equation 
involving only the active neutrino v-y L '- 

^^^N aR = id^ MBi, {iduY^ L ) = E M L^l- (71) 

7 7 
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Following the treatment of Kuo and Pantaleone in p0[ , we make the following simplification using 
the properties of the gamma matrices: 



drfdvY = - + dudrf-f] = dpdvgr = d^. (72) 

With this substitution, equation ( O ) may be written as 

d„d^j3 L + J2(MnMl)^ 1L = 0. (73) 

7 

This has the same form as a Klein-Gordon equation. Following the traditional treatment of a 
Klein-Gordon equation, we assume a plane wave solution of the form 

^^) = E e " lPV ^. (74) 

Vfi is a /3-flavor neutrino spinor, with no time or spatial dependence. The 4-momentum is repre- 
sented by a matrix because it is not diagonal in flavor space. Eventually, we will diagonalize the 
4-momentum and move into the mass basis, but first a word about the plane wave assumption. 

As mentioned above, a given flavor of neutrino created at a source does not necessarily have 
a definite mass, and therefore the 4-momentum is not diagonal in flavor space. Most discussions 
of neutrino oscillations, such as those found in JT3] and jl3], treat the neutrino of definite flavor 
as a particle with definite momentum, composed of a superposition of energy states corresponding 
to the different mass states. Reference |H] shows that this treatment yields the same results as 
a treatment using a definite energy and a superposition of momentum states, and that both yield 
results consistent with a more rigorous treatment of the problem using wave packets, provided the 



neutrinos are relativistic. Appendix C.l summarizes the wave-packet approach and the validity 
of its approximation by the plane wave solution (|74|). The energy-momentum assumptions of the 
standard oscillation treatment are examined in reference p2|, and we present those arguments in 



Appendix C.2. In what follows, we will use the plane wave solution (74) rather than wave packets 



for simplicity, but we will use the relativistic results of Appendix C.2 to conserve both energy and 
momentum. 
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Examining the plane wave v aL in equation (174) , we find 



d d°is aL = - ^ ( E afi) 2 "PL' and 

0= flavor s 



(75) 



a— 1 0— flavors (3 



(76) 



Note that the index a is not summed over in equations ( [76] ) and ([75]), but a is. As usual, a and /3 
are flavor indices, and a is a spatial Lorentz index. The conventions we use for summing and indices 
are described in Appendix |X| Using the plane wave solution in (73), we obtain an equation for the 
matrix p^p a /3 M . 

d^v aL = J2 {E 2 afj - \p\l p ) v 0L = Y,( M ° M dU v 0l- (77) 

(3 

Neutrinos propagate in mass states, not flavor states, so we must diagonalize the mass matrix in ( |77j ) 
to find an equation of motion for the mass states. The same mixing matrix V diagonalizes both the 
mass matrix and the four-momentum matrix, leaving us with an expression for the four-momentum 
of the ith mass state Vi L : 



a, /3,7 



(3 V PL 



£ Vl(M D M^ D )^{VV^) ip v^ so 

EE 

a, 7 i 



Vl { M D Ml , , , 

Q7 



v; 



E^Wi 



(78) 



with 



i L 



(79) 



and 



^2vl(M D MlUV yj = m 2 i S i: 



07 



(80) 



The matrix on the right-hand side of equation (|78|) is diagonal, so the matrix on the left-hand side 
must be diagonal too. Consider first the rest frame of the neutrino. The three momentum is zero, 
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so the left-hand side of equation (y8[) is 



E 

a, 7 



\vt (E 2 ) 

[ JQ \ cry J 



ViL = EiSi 



(81) 



Making a Lorentz boost back into the moving frame only changes Lorentz indices, not flavor indices, 
so equation ( |8l| ) remains valid in the frames with non-zero momentum. Since the left-hand side of 
equation ( |7q ) must be diagonal and the first term is diagonal by equation (|l]), the second term 
must be diagonal too: 

— r ' 1 (82) 



Comparison of equations (|7q) through (p2|) yield the familiar energy-momentum relationship 



{E 2 -\p l \ 2 ) VlL =m 2 u 



iL- 



(83) 



Having rotated into mass states of definite energy and momentum, we may now describe the 
propagation of the neutrino mass eigenstates through space: 

u iL {t h x i ) = e i ^-*^v i . (84) 



We have placed mass indices on the time and position, since the propagation of the state in space 
and time depends on the mass. (Different mass states will travel the same distance in different times 
or different distances in the same time). An equation similar to (|4|) may be derived for the sterile 
neutrino Nr, and the two may be combined in vacuum to describe the behavior of a Dirac neutrino 
v. 

3.2 Vacuum Oscillations 

Weak interactions at a neutrino source produce neutrinos of a definite flavor, v a , so the wave- 
function of a neutrino when it is produced is 

v L (t = 0,x = 0) = v aL = ^2 VaiV iL , (85) 
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The wavefunction a distance x from the source, and a time t after emission, (taking the x-axis along 
the momentum direction) is 

n n 

v L (t,x) =Y,y^ie- i{Eiti - piXi) =J2v<*me-^, (86) 

i=l i=l 

where <pi = Eiti — piX% is the kinematic phase of the propagating mass state V{. Although neutrinos 
travel in mass states, they are detected in flavor states, so we must rotate back to the flavor basis 
to compare with experiment: 



(3— flavor s Li— 1 



vp = y^^apgVp- (87) 





ap a is the amplitude for the transformation v a — > vg, and the probability of the transformation is 



P (x) = \a 0a \ 2 =J2T,( V ™ V k K Vf>j)e<+*-*A 

Va ~* V& i=l j=l 



We call the phase difference the relative phase 



= 5 - . (89) 



where the factor of i is included for future convenience. 
Because the mixing matrix is unitary, 



]T (v*tVpiV2jVf>j) = <W> (90) 
1=1 3=1 



and (pq) may be written as 



p (*) = E E Wi^i) ( e ~ l2 * ,j - !) + <w- ( 91 ) 



"a — V P 

i=\ j = l 



Terms antisymmetric under the interchange i <-» j will cancel in the double sum, leaving only 
symmetric terms. The product of matrix elements in equation ( |9l| ) becomes its complex conjugate 
under this interchange, so the real part is symmetric and the imaginary part antisymmetric. The 
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exponential may also be broken into its symmetric cosine and antisymmetric sine contributions. 
Keeping only the symmetric products, equation ( |9l|) becomes 



P (x) = yV[Re(^ m ^K^)(cos(2$ y )-l) 
+ Im (VctVjiiKjVpj) sin (2$ y )] + S aP 



E [ Re {VoiVpiKjVpj) 2 sin 2 (92) 
i=l j=i 

- Im (VaiVpiKjVpj) sin (2$y)] + <5 Q/3 . 



The terms for which i = j do not contribute since $jj = 0, so we remove them to arrive at our final 
form for the neutrino oscillation probability: 

p (x) = -2j2Y,M v ^ v 0^J v ^ sin2 ^ ( 93 ) 

+ E E Im (K»V£ V£ sin (2$ y ) + <5 Q/3 



For relativistic neutrinos, is derived in Appendix C.2 to be 



2 2 

m~ — m~ 



- — ( 94 ) 



which is linear in x, the distance traveled. As a neutrino propagates, its transition probability thus 
oscillates in distance as the square of a sine function with an oscillation length 

L-- = — = 4lTp = 4?rp (95) 
13 ~ $y m? - m 2 ~ Am?- ' 

Am?- = mf — m 2 is called the mass-squared difference for mass states i and j. 



We see upon examination of equations (|93j) and (94) that the oscillation probability when all mi = 
rrij is 5 a p. Thus if the neutrino masses are all zero, the probability for a neutrino to change flavors is 
and the probability for the neutrino to stay the same flavor is 1. Detection of neutrino oscillations 
would therefore be proof of non-zero neutrino masses and potentially the first contradiction of the 
Standard Model. 
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3.2.1 Two- Way Oscillations 
Consider the simplified case of two neutrino flavors, n = 2. An arbitrary 2x2 unitary matrix 
will have one angle parameterizing the real degree of freedom, and three phases corresponding to 
the three imaginary degrees of freedom. But all three phases may be absorbed into the definitions of 
Dirac fermion fields (or cancel in oscillation probabilities for Majorana neutrinos as described later 



in Section 4.2.2 ), so a 2 x 2 mixing matrix has only one (real) degree of freedom, which we will 



describe by the angle 9. The mixing matrix V then has the form 



cos 6 —sin ( 

V = I I • (96) 



sm f cos ( 



Since the matrix is explicitly real, the term in equation (93) involving imaginary matrix elements 



disappears, and the oscillation probability in the two flavor case is simply 



P (x) = +4 cos 2 9 sin 2 9 sin 2 ^-x + 5 a0 



sin 2 26» sin 2 ( ) <V .■ (!I7) 



3.2.2 Three- Way Oscillations 
The three-generation case provides a bit of freedom in parameterizing the mixing matrix. An 
arbitrary 3x3 unitary matrix has three real degrees of freedom, and six phases. As will be shown in 



Section 4.2.2 , 2n — 1 = 5 relative phases are not observable, so the three-generation mixing matrix 
needs three mixing angles and one phase to describe it. The quark sector contains such a matrix, 
called the CKM matrix after Cabbibo, Kobayashi, and Maskawa, who developed the theory of quark 
mixing. The original choice of these four parameters, by Kobayashi and Maskawa, is perhaps the 
best known parameterization. Their choice of V is JllJ 



ci 



sic 3 



S1S3 



\ 



-8\C2 C1C2C3 - S2Sze lS C1C2S3 + S2C3e lS 



-S1S2 C1S2C3 + C2Sze lS C1S2S3 - C2Cse lS 



(98) 
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where c a = cos# a , and s a = sin^a. We will refer to this particular choice of angles as the "standard" 
or "KM" parameterization, while the name "CKM matrix" refers to any three-generation mixing 
matrix parameterization. 

In the standard KM parameterization, the phase only appears in the lower right-hand "corner" 
of the matrix. The placement of the phase is somewhat arbitrary, since we absorb five relative 
phases into the field definitions. Choosing to absorb a different set of relative phases repositions 
the remaining phase in the mixing matrix. Clearly the location of the phase cannot be measurable! 
Indeed, as discussed in the next chapter, individual mixing matrix elements cannot be measured 
independently, and the measurable probabilities include a single function of the phase, called J , not 
the phase itself. In the standard KM parameterization, the invariant J has the form 

J = Cjsf C2S2C3S3 sin S. (99) 



This function has a maximum value of ^= |35| ] : 



l 3 ) 

I max V 2 



J max = (cosfli - coBfli'l - sin2tf 2 - sin20 3 (^§) max = ^ (1) = 



1 . \ , . « 2 1 1 , s 1 



9 / ^ 'max q fo 9 (i fo' 

max \ Z /max J V O * * DVo 

(100) 

Another popular parameterization chooses slightly different mixing angles and places the phases 
elsewhere in the mixing matrix 0], p3| : 

C12C13 S12C13 s 13 e 



-C23S12 - Ci 2 S 2 3Sl3e 1 ' 5 C12C23 - S 12 S23Sl 3 e lS C13S23 

\ Si 2 S 23 - C 12 C 2 3Sl 3 e l5 -C12S23 - C 2 3Sl2S 13 e lS C13C23 J 



(101) 



Both of these parameterizations (|98|) and (|10l| ) are equally valid, but they are not equivalent. In 
addition, we will demonstrate in Chapter |y| how complicated the observables are in terms of the 
standard KM parameterization. Our development of a model-independent parameterization has 
been motivated by the arbitrariness and complexity of the traditional approach. 
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3.3 Matter Effects 

When neutrinos travel through matter, they may interact with the electrons, protons, and neu- 
trons contained in the matter. Once these interactions are averaged, the resulting terms in the 
effective Lagrangian resemble mass terms, as shown below. All flavors of neutrinos will scatter off 
of electrons, protons, and neutrons by the neutral current interaction, so the mass matrix receives a 
contribution proportional to the identity matrix which does not affect the mixing matrix. (Rotating 
this identity term by the mixing matrix leaves it diagonal while diagonalizing the mass terms. Thus 
the same mixing matrix that diagonalizes the vacuum Lagrangian will diagonalize the Lagrangian 
with neutral-current terms added.) However, only electron-type neutrinos and electrons may scatter 
via the charged current. This type of interaction produces a non-trivial contribution to the mass 
matrix, affecting the mixing matrix and oscillation parameters. 



As discussed in Section 1.2, neutrino-electron charged current scattering may be represented by 
the effective Lagrangian of equation (pOl). Substituting for J cc , we obtain 



—AG 

Ceff = [WLi^l] [eZj^eL] ■ (102) 



The order of the spinors may be changed by a Fierz transformation (see, for example, Appendix F 
of reference Jill or Exercise 2.12 of reference Il(| for details on these transformations), resulting in 



4G 

C c e f f = --j=[V^^v eL ][eZ llx e L }. (103) 



The electron factor may be reduced by averaging over the electron field and assuming non-relativistic 
electrons Q to give 

{eZYe L ) « (eZ 7 °e L ) = {e\e L ) = n eL « -n e , (104) 

where n e is the electron number density in the medium interacting with the neutrinos, and we assume 
that the medium contains equal numbers of left- and right-handed electrons. Thus the contribution 
to the Lagrangian from the charged current neutrino-electron scattering is 

C°f s = -V2Gn e u^7oVeL- (105) 
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A similar treatment may be applied to the neutral current contribution to the effective Lagrangian, 
resulting in 

(106) 



V2G 
~2 



V2 



where the factor of ^ comes from the difference between charged and neutral currents, and we have 
set n e = n p for electrically neutral matter. As we shall see, this neutral current term does not affect 
neutrino oscillations since it is uniform in flavor. 

Because these effective Lagrangian terms contain only the spinor products Vv, they are similar 
to mass terms and affect the diagonalization of the mass matrix. When these interaction terms are 
added to the vacuum neutrino Lagrangian (pq), equation (69) becomes 



6H 1 

- — = id^faL - X! M D a p N i3R. + -j=Gn n j is aL - V2Gn e S ae ^ v aL = 



8v a 



(107) 



while equation (frOl) remains unchanged: 



5C 



The equation of motion for the active neutrino (|71|) becomes 



id^N aR = M Da-y id ^ (idrf + -^=Gn n j - ^/2Gn e 5 iel ^j v. 
= V M f n u 

/ j JJaj 



(108) 



If the densities n e and n n do not fluctuate quickly in space or time, we may approximate 



(109) 



This approximation is valid if % < < 1 
From the Dirac equation we know 



d^Va L ~ m v v a L 0. 



(110) 
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The flavor states don't have definite mass, but we can assume m„, the effective mass of a flavor 
state, to be on the order of the individual masses mi, and therefore small. From this approximation, 
we find 

dol°VaL ~ d a j a V QL , SO (111) 



d^YloVaL = 7o(3o7° + d al a )v aL ~ 7o(25o7°)^ qL = 2d v al 



(112) 



With these substitutions and the simplification of equation (72), equation (|108|) becomes 



apVt3 L + iGy2 (2n e 5 ae — n n ) d v aL = 0. 



(113) 



This looks like the Klein-Gordon equation ( [73] ) with an amount G\/2<9o(2n e <5 Qe — n n ) added to the 
diagonal elements of the mass matrix. Again we try a plane- wave solution of the form ( ]74| ) . The 
energy-momentum relationship, however, is affected by the matter, becoming 



^(E*«- \P\1 ) Vf3 r = X 







'L / j 





(m d mI) + iGV2 (2n e S ae - n n ) {~iE aP ) 

\ / a/3 



v pi 



(114) 



The neutral current term will automatically be diagonalized by the same mixing matrix which 
diagonalizes E a p, but the charged current term will affect the diagonalization. We may define a 
new mass matrix M M which includes this charged current term and is diagonalized by the matter 
mixing matrix, V M : 



V, 



(MnMjjU + 2E aP V2Gn e S ae ] V$ = (M M M M ^) = (mf) 2 (115) 



where mf 1 are the effective neutrino masses in matter. The mixing matrix is independent of the 
neutral current term, as predicted before. As for the vacuum case, the momentum and energy are 
also diagonalized by the mixing matrix: 



v£*\P\lpVj% = and 



T/Aff p2 T/M _ P 2r 



(116) 
(117) 
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The oscillation equations in matter are the same as they were for neutrinos traveling through a 
vacuum, provided the mixing matrix V is replaced with V M , and the masses mf are also replaced 
with their corresponding (mf 1 ) 2 . For example, equation ( |93| ) would become 

P ™M = - 2 EE Re (^f*^T O sin2 ( — )2 4 " — ) ( 118 ) 
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CHAPTER IV 



A NEW PARAMETERIZATION 

4.1 The Unmeasurables 

The individual elements of the mixing matrix V which appear in the Lagrangian are unmeasur- 
able. Consider the two different parameterizations of the three-generation mixing matrix presented 



in equations ( |98|) and (101) of Section 3.2.2 . The freedom to absorb phases into the definitions of 
fermion fields produces a freedom in the placement of the remaining phase. For a Dirac mass term, 
we may absorb 2n phases into the definitions of the n left-handed neutrino fields and n left-handed 
fields of the charged leptons. Such an absorption removes In — 1 relative phases from the mixing 
matrix without changing any other terms in the Lagrangian, provided we absorb conjugate phases 
in the right-handed fields: 



C = VZM„N R + l L Mil R + c 2 W+Wy tl lL + h.c. 

= ^D I/ N^+T^M l l^ + c 2 W+l^Vj^ + h.c. (119) 
= ^X V X- X D V X N X^N^ +TfX lL X[ L 1 M l X ln X[ R 1 l% 
+ czW+vfXvX^VX^X^ri? + h.c, 



with 



; = Sye-'^i (120) 



Setting Xn ~ X v and Xi R = Xi L , which leaves the diagonalized matrices D v and M\ unchanged, 
and redefining ip' = X^ lr ip m , we get 

1^JD V N R ' + T L Mil' R + c 2 W+^J [A,; 1 VX lL ] ^I'l + h.c. (121) 

So we may absorb n phases through X v and n phases through Xi L in the field definitions, thereby 
removing 2n — 1 relative phases from the mixing matrix, without changing the rest of the Dirac 
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neutrino Lagrangian. Note that it is the phase differences between X~ x and Xi that enter into the 
new mixing matrix, and so 2n — 1, rather than 2n, phases are removed from V. 

If the neutrinos are Majorana particles, however, we do not have the freedom to set the phase 
matrix Xn independently of X v [jl4|. The Majorana mass term has the form 

vZM v V R , with V R = CvZ T - (122) 

Making the transformation — ► VjJ = X~ 1 V[ j m results in the transformation i> R — > V R = 
X~ 1 *V R m — X v P R m ; the phase matrices in the mass term do not cancel, but combine, and the 
Lagrangian is not invariant under such a transformation: 

VEM V P R ^ U?D V V R '. (123) 

Therefore phases may be absorbed only in the charged lepton fields, and just n; (which may be 
different from n = n v in models with sterile neutrinos) phases may be removed from a Majorana 
mixing matrix. 

These In — 1 — n; extra phases in a mixing matrix for Majorana neutrinos are not, however, 
observable in oscillation experiments JLlj . Each row and each column of V may be multiplied by an 
arbitrary phase; id est, any mixing matrix element V a i may be replaced with V^e - ^ without 
changing the product of four matrix elements found in the oscillation equations: 

VaiV^VfrVn = e^V^e-^e-^V^e-^V^e^V^e-^ . (124) 

Thus whether the mass matrix is Majorana or Dirac, In — 1 relative phases will not appear in 
the oscillation probabilities even if the number of neutrinos is greater than the number of charged 
leptons. Processes sensitive to the "extra" Majorana phases include double-beta decay and the 
radiative decay of unstable neutrinos contained in some models; in these processes the same vertex 
occurs twice (as opposed to a vertex and its hermitian conjugate) in the relevant Feynman diagrams. 
A discussion of neutrinoless double-beta decay and its relation to Majorana masses is contained in 
Appendix 
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4.2 The Measurables 

The immeasurability of the mixing matrix elements in the quark sector has been addressed 
by numerous authors, such as those of references |35|| , p6f , |37f| , and ]3S|j . Measurable quantities 
include only the magnitudes of mixing matrix elements, the products of four mixing matrix elements 
appearing in the oscillation equation (B3J), and particular higher-order functions of mixing matrix 



elements 36 . Neutrino oscillation probabilities depend linearly on the fourth-order objects, 

= V ai y^V 0j Vl = VaV^V^Vpj, (125) 

which we will call "boxes" since each contains as factors the corners of a submatrix, or "box," of the 
mixing matrix. For example, 

"□22 =VuV^V^V22. (126) 

These boxes are the neutrino equivalent of the "plaques" used by Bjorken and Dunietz for another 
purpose in the quark sector in reference |37j . In this work, we will examine the applications of using 
boxes to describe neutrino oscillations. 

We have developed a graphical representation (discussed in detail in Appendix ^) which proves 
useful when considering relationships between boxes. Figure || illustrates the representation of 
11 ^22 • This representation will be used throughout this Chapter to illustrate various relationships. 



I t 



Figure 6: The graphical representation for the box 11 D 2 2 in three generations. Vertical arrows 
point from the matrix elements which are not complex conjugated in the box to the complex- 
conjugated elements. 



In general, each of the box indices i, j, a, and (3 may be any number between 1 and n, the number 



of neutrino flavors. We therefore initially have n 4 possible boxes. Examination of equation (125) 
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however, reveals a few symmetries in the indexing: 



= = ^a* aj = . (127) 



If the order of either set of indices is reversed [id est, j «-> i or /3 <-> a), the box turns into its complex 



conjugate; if both sets of indices are reversed, the box returns to its original value. Equation (127) 



demonstrates that boxes with a = (3 or i — j, are real. Indeed, these are given from equation (J125|) 



as 



ai a aj = \V ai \ 2 \V aj \ 2 , and (128) 
m u pi = \V ai \ 2 \V 0i \ 2 . (129) 

Those boxes with both sets of indices equal are 

ai U m = \V m \\ (130) 

We call boxes with one or two repeated indices "singly-degenerate" and "doubly-degenerate," 
respectively. As can be seen from equation (^3|) , singly-degenerate boxes with repeated flavor indices 
enter into the formulae for flavor-conserving survival probabilities, but not for flavor-changing tran- 
sition probabilities. Degenerate boxes with repeated mass indices (including the doubly-degenerate 
boxes) do not appear in any oscillation formula. Degenerate boxes may be expressed in terms of 
the nondegenerate boxes with a ^ j3 and i ^ j, as will be shown shortly. This possibility and the 



symmetries expressed in equation (127) allow us to express combinations of boxes in terms of only 



the nondegenerate "ordered" boxes for which a < (3 and i < j. 



n 



The number of flavor-index pairs satisfying a < ordering is N = | I = N mass- 

2 

index pairs similarly satisfy i < j ordering. Thus, the number of ordered nondegenerate boxes is 
TV 2 . The number of flavor-degenerate boxes is n 3 , and that of mass-degenerate boxes is also n 3 . The 
n 2 doubly-degenerate boxes appear in both these counts, so the total number of degenerate boxes 
is 2n 3 — n 2 = n 2 (2n — 1), implying n 2 (n — l) 2 nondegenerate boxes. From the nondegenerate boxes, 
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we want to select only the ordered ones for which j > i and (3 > a, which comprise one fourth of 
the nondegenerate boxes. The number of nondegenerate ordered boxes is therefore ^ "("^ ; as 
argued above. Applying ordering to the singly-degenerate boxes gives a count of \v? (n — 1) = nN 
for flavor-degenerate boxes, and the same for mass-degenerate boxes, yielding 2nN total singly- 
degenerate ordered boxes. For n = 3, N = 3 too, and we have nine ordered nondegenerate boxes, 
nine ordered singly-degenerate boxes, and nine doubly-degenerate boxes. These counts are recapped 
later in Table [I], along with other box counts. 



Using the symmetries expressed in equation (127), equation (B3J) becomes 



p W = -EE[ 2Re ( MD ft-) sin2$ -r Im ( tt ' D ft> ill2$ . 3 ] 



% = 1 j > i 



- E E t 2Rc ( ) sin2 - Im ( Q<D « ) sin 2 *«] ( 131 ) 

i=l j<i 



= -EE [ 2Re ( + mu h ) sin2 - Im ( mn « - mu h ) sin 2c M + * 



a/3 



2—1 j>i 



= -2 E t 2Re ( m D « ) sin2 $ « - Im ( ai a P'j ) sin 2$ «] + s <*e > 

with $>ij defined in equation (p4T). The survival probabilities P (x) may be found from the transition 



probabilities P (x) by 

Va, — *■ >/0 



V a — >■ V a 1 * V a — * Vq 



Interchanging a «-> (3 in equation (131) gives the time-reversed reactions P (x) 



P (^) - -2^ E [ 2Re ( ) sin2 ®v + Im ( ) sin 2^-] + S af3 , (133) 



vp ->Vt. 

2—1 ,?>2 



so a measure of T-violation, or equivalcntly CP-violation, in the neutrino sector is 



p (x) -p (x) = 4 E Im ( ma pj ) sin 2$ ^ • ( 134 ) 

v„ va va — > l><v * — ' ' — ' 



"/9 "0 

2—1 J>2 
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Equation (131) may be expressed in matrix form. For three flavors, we find 



7>(n = 3) 



f P (x) ^ 

v e — > 
P (x) 

— > Vt 

\ P (*) / 

\ !/ e — > Kt / 



-4 Re(B) S 2 ($) + 2 Im(B) 5(2$), 



(135) 



where 



B : 



12 n li n 

u 23 u 23 



2 1 1 — i 22 i — i 2 1 1 — i 

u 32 u 33 u 33 



ll n 12n Hn 

u 32 u 33 u 33 



' sin* 3 ^ 



, and S*($) 



J 



sin k $ 



23 



sin k $ 



13 



/ 



(136) 



Here the operation Re(B) is defined so that [Re(B)]jy = Re(Bki); Im(B) is defined in a similar 
manner. B and 5 fc ($) are defined for three flavors in equation ( |136| ), and these definitions may be 
extended to any number of flavors. This extension, however, rapidly becomes difficult to manage. 
Adding a fourth, perhaps sterile, neutrino flavor increases the number of boxes from nine to thirty- 
six, since there are now six independent transition probabilities and six mass differences. For the 
six flavors of neutrinos in the see-saw mechanism, B is a fifteen-by-fifteen matrix, and we have two 
hundred twenty-five boxes! 

If CP is conserved, the boxes are real. CP conservation also implies time-reversal symmetry 
if CPT is a good symmetry, so P (x) = P (x) . These symmetries leave us with just three 

V<x — * "p vp — > v a 

independent probabilities (assuming three flavors): 



— » f M 

P (x) 

Vfi — > V T 

\ P (*) / 

\ V e — * V r / 



= -4 



u 22 u 23 u 23 



21 ^32 22|=l 33 21|=l 33 



V 



n n 12 n 1:L n 

u 32 u 33 u 33 



sin $12 
sin 2 $23 



y sin 2 $13 J 



(137) 



Our three probabilities are, however, still functions of E and x, so different experiments may not 
measure the same values for a given oscillation probability, regardless of whether CP is conserved. 
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4.2.1 Mixing Matrix Elements and Boxes 
Neutrino oscillation experiments will measure the boxes in equation ( |131| ), not the individual 
mixing matrix elements, V a %- But one would like to obtain the fundamental V a i from the measured 
boxes. Some tautologous relationships between the degenerate and nondegenerate boxes are easily 



confirmed using equation (125); they hold for any number of generations: 



\V m \ 2 \V a A 2 



Ql|-|* Oil 



'□ 



n 



"J 



a), 



(138) 



W, 



l2 \V 8 . 12 



ai\ \vf3i\ 



04q* Oil 



/3i 



fix 



□ 



/3y 



/3y 



(x^y^i), and 



(139) 



|K»i| 2 
l%l 2 



px rjj 



(V ^ a ^ P, and x ^ i ^ j). 



(140) 



In these equations and what follows, a, f3, 7, i, j, and k will usually be reserved for indices that are 
chosen at the start of a calculation, while other indices such as x, y, 77, and A primarily represent 
"dummy" indices which are chosen arbitrarily except to respect the inequality constraints following 



equations such as equation (140) 



The tautologies (138) to (140) become evident not only by algebraically using the definitions 
of the boxes, but also by considering the graphical representation described in Appendix [e]. For 
example, the degenerate box 12 Di3 is "created" when one uncanceled vertical arrow enters and 
another leaves each of the matrix elements V12 and V13, while arrows at all other points cancel. This 
is illustrated in Figure |a. Those arrows are the result of the combination of ordered boxes given in 
(138), as shown in Figure 0b. 



Equations (138) and (Jl 39j) are themselves special cases of the more general 



(141) 



and the analogous relation 



Cxi 1 — 1 ak 1— 1 



aii — 1 akr— i 



131 ^fij ■ 



(142) 



4G 



13 
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(a) 



i i 



(b) 



Figure 7: The graphical representation in three generations for (a) the degenerate box 12 Di3 , 
and (b) the combination of nondegenerate boxes which "create" it. 



These relations hold for both degenerate boxes and nondegenerate boxes. Including disordered 
boxes adds no new information, so we will as usual consider only ordered boxes on one side of the 
equation. (Swapping indices to produce the other side of the equation may lead to disordered boxes 



on that side, but that can quickly be remedied by using equations (127)). If a = 7 or j3 — 5 for 



equation (141), or if % = k or j = I for equation (142), the relationships are trivial. Ordering of these 



pairs, however, should not be imposed. 

Consider the second relation, (|l42| ). Each of the \n[n + 1) ordered plus degenerate pairings of 
(a, 0) contains j^n(n + l)(3n — 2)(n — 1) nontrivial relations with j 7^ I, i =/= k, i < j, and k < In 
all, equation ( |142| ) implies ^n 2 (n + \) 2 (n — l)(3n — 2) nontrivial relations, as does equation ( |141[ ). 
On the other hand, if we restrict ourselves to considering only the N boxes nondegenerate in (a,/3), 
and (k, I), the number of constraints drops to ^71(71 — l)( rl ~ 2)(3n — 5) j 7^ I, i 7^ k, i < j, 
k < I combinations for each of the N a < j3 pairings ,^ leading to a total number of nondegenerate 
constraints given by in 2 (n — l) 2 (n — 2)(3n — 5). Degenerate boxes participate in the majority of 



the relations implied by equations (141) and(|142[). For n — 3, the number of i, j, k 7 I combinations 



1 We arrive at this number by considering all (in(n + 1)) (i < j, k < I) pairings, then subtracting the ^n(n + 1) 
(I = j, k = i, i < j) possibilities, followed by the ^ ™ J = |n(n-l)(n-2) (1 ^ j, k = i, i < min(j, !)) combinations 

and the same number of (j = I, k 7^ i, j > max(i, ft)) combinations. We must also subtract the 2 ^ ™ ^ = n ( n — 1) 
(I 7^ j, k = i = min(j',Z)) possibilities and the same number of (ft 7^ i, j = I = max(i, k)) possibilities. So the number 
of terms in the sum X/j^i i^k i< ' k<l ls (,\ n ( n ■O) ~ \n(n + l) — 2^n(n — l)(n — 2) — 2n(n — 1) = j^n(n + l)(3n — 
2)(n-l). 

2 To obtain this number, we remove the i = j and k = I possibilities. So we start with only fin(n— 1)) 

(i < j, k < /) pairings and subtract |n(n — 1) (J = j, k = i, i < j) combinations followed by ^n{n — l)(n — 2) twice 
as above. The 2n(n — 1) pairings for [I 7^ j, k = i = min(j, I)) and (k 7^ i, j = I = max(i, k)) do not occur in our 
original sum since % 5^ j and k 7^ I in that sum. The number of terms in the sum . , , is therefore 

(|n(n - l)) 2 - in(n - 1) - 2±n(n - l)(n - 2) = ^n(n - l)(n - 2)(3n - 5). 
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in equation (142) including degenerate boxes is 14; without degenerate boxes this number decreases 
to 2. These numbers are doubled once the relationships implied by equation (|l4l| ) are considered as 
well. 



This index rearrangement of equations (141) and (142) may straightforwardly be extended to 
obtain higher-order relationships (trilinear, quadrilinear, et cetera) in the boxes. These higher-order 
relationships do not enter in our treatment of oscillation probabilities, so we will not examine them 
here. 

i 

We may find \V a i\ — ( m O ai ) 4 in terms of three singly-degenerate boxes by setting a = (3 



in equation (|139|). Using equation (138) to substitute for the singly-degenerate boxes yields an 



expression for the doubly-degenerate box in terms of nine nondegenerate boxes: 



\ V at\ — LJ (ij, — 



"■n m n 

'—'ax ay 



'□ri 



L - '<TX l — 1 pi 



I [—1 UX [—1 



ay 



ifl Pifl * at/n axri 
^ax F u Ck u ui u /jy 



t ^ a ^ a 

pi 7^ lo ^ a 
x^y^i 



(143) 



As noted earlier, flavor-degenerate boxes give the probabilities for flavor-conserving oscillations, 
while nondegenerate boxes give the probabilities for flavor-changing oscillations. Thus individual 
\V a i\ may be deduced from a set of measurements of cither kind. One can obtain a relationship 



similar to equation ( |143| ) by setting i = j in equation (13£) and then using equation (139) to 
substitute for the singly-degenerate boxes: 



\v -I 4 — ai a ■ — 

I v ai | — L - 'ai — 



-'Tjl 



Xl D r 



1 -^pi — i 

J \n u ap 'r]r '—'as '—'Xu 

i«n v r n A»n 



n =/= p =/= i 
r ^ s ^= i 
t u =/= i 



(144) 



In the three-generation case, these equations are unique and equivalent to each other, since all 
nine boxes are used in both equations. This equivalence may be seen by choosing s = n = u = x, 
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r=p = t = y,(T = p = uj~X, and r = C = A* = Wi which yields for either expression (143) or (144) 



l^iil 4 



"□a, 21 n 13 11 □ 



33 



"□ 22 "nja "Daa 22 n 33 



22 D 13 33 Dl2 21 D 33 31^ 



12 D * 3 12 D33 21 D33 21 □ 



2l n * 

32 



(145) 



for the case a — i = 1. The second equality merely represents a substitution of ordered boxes for 
disordered boxes. Two other examples are 



1^21 1 



2l n llr-i 2l n n n 12 n 
U32 U22 U33 U23 U33 



"□ 32 "Dsa 22 n.33 12 □ 



J 33 



33 



J 23 



and 



(146) 



|V 3 i| 4 



"□ 32 21 a h aiQga i 2 n 2 



"□22 n ^ 3 12 n 33 22 n 



22 n * 
33 



(147) 



The equalities (145), (146), and (147) are illustrated in Figure §1 



(a) 



(b) 



(c) 



Figure 8: The graphical representation for (a) |Vn| 4 , (b) |V2i| 4 , and (c) | V31 1 4 . Notice that 
each matrix point except the one being represented has an equal number of vertical arrows as 
horizontal arrows entering and leaving. The element being represented has two vertical arrows 
leaving (representing V^) and two entering (representing V*f) to produce |K«| 2 . Appendix [|] 
describes how these were produced, using | V21 1 4 as an example. 



Before continuing, we should note that all of the relationships in this Section follow from the 



definitions of the boxes in equation (125) and so are valid for any matrix, unitary or otherwise. The 
unitarity of V will impose further relations among the boxes which may not hold for models which 



only approximate unitarity of the mixing matrix. Equations (138) to (|142| ), however, will always 
hold. 
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4.2.2 Using Unitarity to Reduce the Number of Independent Boxes 
An n x n arbitrary matrix has n 2 elements, which may be complex, leading to 2n 2 parameters. 
Our mixing matrix, however, is unitary, and has n 2 unitarity conditions as constraints. Of the n 2 
remaining parameters, 1) (the number of parameters in a real unitary, or orthogonal, matrix) 

we may 



are real, and the remaining jn(n + 1) are imaginary phases. As discussed in Section 4.1, 
remove 2n—l relative phases from the mixing matrix, leaving ^n(n + 1) — (2n— 1) = ^(n — l)(n— 2) 
phases measurable in oscillations. If the phases arc all zero, all of the boxes are real, so CP violation, 



as expressed in equation (134) disappears. Similarly, if any of the phases are not zero, CP is violated 
in the neutrino sector. 

The total number of CKM mixing matrix parameters is n 2 — (2n — 1) = (n — l) 2 . The number 
of real and imaginary independent box parameters should be the same as the number of real and 



imaginary CKM parameters 1 39 , but the reduction to a basis is complicated for boxes. The number 
of ordered boxes A 2 grows much faster than the number of mixing matrix elements n 2 and exceeds 
that number once n > 3, so we need more box constraints than mixing matrix constraints for the 
higher numbers of generations. The parameter counts for the mixing matrix and the boxes are 
summarized in Table [7j. 

Table 1: Parameter counting for the mixing matrix and boxes. N 2 is the number of ordered, 
nondegenerate boxes, with N = ^n(n — 1). 



Number of generations 


n 


2 


3 


4 


6 


Params. for arbitrary matrix 


2n 2 


8 


18 


32 


72 


Unitarity constraints 




4 


9 


16 


36 


Relative phases 


2n-l 


3 


5 


7 


11 


Real params for unitary V 


N 


1 


3 


6 


15 


Remaining phases in V 


i(„-l)(„_2) 





1 


3 


10 


Initial boxes 


n 4 


16 


81 


256 


1296 


Doubly-degenerate boxes 


■> 
n z 


4 


9 


16 


36 


Ordered singly-degenerate boxes 


2nN 


4 


18 


48 


180 


Ordered nondegenerate boxes 


N 2 


1 


9 


36 


225 


Independent Im( m D,3j )s 


i(„-l)(n-2) 





1 


3 


10 


Independent Re( m ^pj )s 


A 


1 


3 


6 


15 



Unitarity requires that 



Ew 



6 ij , and 



(148) 
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(149) 



These equations are not independent sets of constraints at this point.^ Equation (148) states V^V =4 



and equation (149) states VV* =1, which is implied by the previous expression and the associativity 
of matrix multiplication: 



VV 1 =4 =S> (VV*)V = V(V^V) = V =S> (t/V) A 



(150) 



We can, however, use these equivalent equations to obtain two separate sets of constraints on boxes 



by multiplying equation ( |148[) by V^V\j and equation ( |149j ) by V* x Vp x : 



n=i 



V&VxjVM^VZVvSv, so 



(151) 



^ "'Daj = ^ Xi a Xl Sij, and 



7J=1 



XI V ™ V <*y V Pv = V «x V Px S <x0, OT 



I/=l 



(152) 
(153) 



^px — V u ai y ap- 

»=1 



(154) 



Separating the manifestly degenerate boxes and the nondegenerate boxes, equation (152) becomes 



77^A 



(155) 



and equation (154) becomes 



y^x 



(156) 



The riiiTri 



ber of unitarity constraints n becomes clear from examination of either equation (148) or (149) 



Equation (148) gives a real constraint for each of the n i = j choices, TV constraints from the real parts of the sum 
when i ^ j, and TV constraints from the imaginary par ts of the i ^ j sum. Thus we have n + 2TV = n 2 unitarity 
constraints from either one of the equations ( p.4q ) and (149). As argued next, the two equations are redundant, so 
the total number of constraints on the mixing matrix is just n 2 . 
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Summing equation (152) over A in the i ^ j case, we find 



o = EES = E 

A=lr;=l A=l 



'□ 



\? + E " iD ^ 

7,#A 



A=l ?;<A 



(157) 



A=l 



Comparison with equations (136) and (13?) shows that the sum of a column of Re (23) which has a 



fixed i and j should equal — ^ 53a=i I^A*I I^Ajl ■ We mav i n a similar manner show that the sum of 
a row of Re (23) with fixed a and 29 equals — \ J22=i \ Vax\ 2 \Vp x \ 2 . 

An alternative way to obtain constraints from unitarity is to start with the definition of the boxes 



(125) and use the unitarity of the mixing matrix: 



rj=£a 



hi - E v * v > 

A^/3 



Xi 



(158) 



After a bit of algebra, this becomes 



EE " iD Ai = mn « -M-i + i^f + iiW 



(159) 



This relation is equivalent to using our result (|152j), as shown below 



E E 

A#/3 njta 



A#/3 Vij=1 

^■ei^i 2 -(e mq 



03 



\X=1 



Si j (l-\V m \ 2 )-\V ai \ 2 5 ij + m n M 



(160) 



Equations (152) and (154) each present \v?{n + 1) relations among the N 2 complex ordered 
nondegenerate boxes, the 2niV real ordered singly-degenerate boxes, and the n 2 doubly-degenerate 
boxes. Not all of these constraints are independent. We can identify one source of redundancy by 



summing equation (152) over A: 



E E ^ = E = J2 \vxiWi = % 



(161) 



A=l ti=1 



A=l 



A=l 
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which is pure real, so 



EE^ = °> 

A=l TJ=1 



(162) 



where we have introduced the notation tfpj = Im ( "'D^- ). We will similarly define 



Since the imaginary sum in equation (162) vanishes, at most n — 1 of the values for A give 
independent imaginary constraints. The real constraints exhibit no such obvious dependency, so 



all n values of A appear at this point to offer independent constraints. Equations (152) and (154) 



therefore each offer at most ^n(n — l) 2 imaginary constraints and ■^n 2 (n — 1) real constraints for 
i 7^ j and a ^ j3, respectively. 



The number of unitarity constraints from equations (|152|) and (154), both real and imaginary, 
grows as n 3 , while the number of ordered boxes grows as rt 4 , so additional relationships between 
ordered boxes must be used to identify the (n — l) 2 independent Js and Rs. These additional 



identities will come from the definitions of the boxes, and the resulting tautologies (141) and (142), 
each of which grows as n 6 , implying a high degree of redundancy. Some of these relationships are 



explored below in equations (J201J) to (205) 



The constraints (152) and ( |154| ) hold independently for the real and imaginary parts of each sum. 
Because the right-hand sides of the equations are manifestly real, as are terms on the left-hand side 



involving degenerate boxes, the sums of nondegenerate boxes in equations (154) and ( 152 ) must be 
real. This leads to imaginary constraints of the formQ 



J2 % = 0. and 



= 0. 



4 The constraints (165) and (166) may be written exclusively in terms of ordered boxes as 

X/ ^ ~ X/ Aj, ' j = °' and 



r?<A 



E 



!7>A 

E 



(165) 
(166) 

(163) 
(164) 



While expressing e very thing m te rms of ordered boxes is our goal, we will find it more convenient to use the complete 
sums of equations (165) and (166) in mathematica ma nipula tion s and switch to ordered boxes at the end rather than 
to deal with the two separate terms in equations (163) and (1164). 
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In three generations, each sum in equation ( |165| ) or (166) contains only two terms, leading to the 
equality (up to a sign) of two Js. For example, choosing A = 1 in three generations yields 



$2i — — ^3i- 



(167) 



Choosing A — 2 relates the imaginary parts of -'□31 and ■'□^j • As discussed above, at least 



one value of A leads to redundant constraints from equation (165); choosing A = 3, which relates 
the imaginary parts of 25 and ^'d^ , is not independent of the first two constraints. For 
each of the n — 1 choices of A not manifestly redundant, we may choose \n(n — 1) combinations for 



i < j. Equation (16E) has hn(n — l) 2 possible index combinations as well, giving n(n — l) 2 total 
imaginary constraints. In three generations, this results in 12 constraints on 9 ordered Js; clearly 



more redundancies among relations (16?:) and (fU36J) must exist! 



The first set of imaginary constraints (165) relates boxes that are part of a column of £>, which 



was defined in equation (|135|); the second set (166) relates boxes that are part of a row of B. Given 



this separation, it becomes clear that the constraints derived from equation (152) and those derived 



from equation (154) are both necessary. Once one combines the two sets of constraints, however, 
more redundancies appear. 



For example, we may relate ¥22 to ¥33 by using equation (165) for i = l,j — 2 with A = 2 and 



then A = 3, followed by equation (166) for a = 1, /3 = 3 with x = 2 and then x = 3: 



1 J l 22 — T32 



12 _ n 

J33 — J33- 



(168) 



This path is indicated in Figure g by the solid arrows. Note that the matrix points in this Figure 
represent boxes in the matrix Im(£>), not mixing matrix elements. We may alternatively use equa- 



tion (166) for a = l,/3 = 2 with x = 2 and then x = 3, followed by equation (165) for i = l,j = 3 
with A = 2 and then A = 3: 



J22 — %3 — — J23 — — T33 — X; 



J 33- 



(169) 



This path is represented by the dashed arrows in Figure 
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Figure 9: An illustration of two different paths relating ¥22 with ¥33. An equality (up to a 
minus sign) between two Js is r epre sented by an arrow connecting those two Js. The solid a rrow s 
represent the path of equation ( 168 ), and the dashed arrows represent the path of equation ( 169 ). 
The matrix represented here is Im(S), not V; this illustration is not an application of the graphical 
representation developed in Appendix ^. 



The number of independent constraints contained in equations (165) and (166) may be deter- 
mined in three generations by counting the number of non-redundant connecting lines in a graph 
like Figure |j| Figure |lO| a contains all 18 connecting lines implied by equations (165) and (166) for 
n = 3. Figure |Io| b illustrates the minimum number (8) of independent lines required to connect 
all nine imaginary parts. Thus eight of our original eighteen constraints are independent, leaving 
us with only one imaginary part as expected, which we'll call J = $11. Using these remaining 
constraints, we find 

J J -J 
Im(B)= J J -J 

\-J -J J ) 

for three generations. 



\ 



(a) 



(b) 



Figure 10: An illustration of a) all imaginary constraints, and b) eight independent imaginary 
constraints. 



For n > 3 generations, the sums in equations ( |165| ) and (166) contain more than two terms, 
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so individual Js become equal to sums of other Js rather than just equal to another individual J. 
Our illustrative method is therefore not applicable to the n > 3 case. In addition, the reduction 



of imaginary parameters using equations (165) and (166) is no longer sufficient to eliminate all the 
extra degrees of freedom. In four generations, the total number of sums relating the 36 Js is 48, but 
we can demonstrate the redundancy of 21 constraints, leaving 9 apparently independent imaginary 



parameters. But we know from consideration of the mixing matrix and from the work of 39 that 
only 3 imaginary parameters are truly independent. The additional 6 constraints must come from 
the expressions relating real and imaginary parameters discussed below. Such a reduction to a basis 
for the four-generation case, while possible, is fodder for another's dissertation and therefore will 
not be presented here. 

The equality of imaginary boxes for n = 3 has an interesting implication for CP invariance. The 



CP-violating contribution to oscillation probabilities (134) becomes 



± 4J (sin 2$i2 + sin 2$ 23 - sin 2$ i3 ) , (171) 

where the plus sign applies for P (x) and P (x) , and the minus sign for P (x) . CP is violated 
if two conditions are met [^3|: i) J ^ 0, and ii) mf ^ mj, for all The first condition should 



be obvious, since a vanishing J clearly makes the expression ( 171 ) vanish. The second condition 



follows rather simply also. If any two masses are degenerate, mf = m?, the corresponding <I>i 



vanishes, and the remaining two terms in equation (171) will be equal but opposite. For example, 
if m\ = m 2 2 , then $ i2 = 0, $23 = j| (m\ - mf), $13 = ^ (mf - mf) = ^ (m\ - mf) = $23, and 
the CP-violating term disappears as predicted. 

Even if the conditions i) and ii) are met, CP- violation may be difficult to directly observe in 
experiments. The CP-violating term contains the sum 

• „^ • • • xAm^ ? . xAmh . x (Am?, + AmL) , „„ N 

sin 2<J> 12 + sin 2$ 2 3 - sin 2<J> 13 = sin — + sin — - sin — i . (172) 

2p 2p 2p 
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Using trigonometric identities, we may write 



x (Amf 2 + Am 23 ) . xAml 2 xAm% 3 xAm\ 2 . xAm\ 



2p 



sm 



2p 



2p 



cos ■ 



2p 



sm ■ 



2:1 



2p 



(173) 



so the CP-violating term is 



sm 



xAm\ 2 
2p 



cAr 



sm ■ 



2p 



xAmin x Am 
sm — cos 



2p 



xAm? 9 

sm — 

2p 



1 — cos ■ 



2p 

xAm 2 , 



cos ■ 



'12 



2p 



sm ■ 



' 23 



2.'! 



2p 



xAm\ 3 



2p 



2p 



1 — cos ■ 



= 2sin2$i2 sin 2 <5> 23 + 2sin2<J> 2 3 sin 2 $ i2 . 



If $i 2 is small, CP violation is of the order 



xAm 2 2 



(174) 



4$i 2 sin 2 $ 23 . 



(175) 



If all phases are small, sin 2$^ ~ 2$^ — f^ij 3 , and the CP- violating term becomes 



(m 2 -m 2 ) 3 + (m 2 -m 2 ) 3 -(m 2 -m 2 ) 3 + • • •) 
= ±3% (( m ? ~ TO l)( m l - ™§)(™3 - ™i)) + ' ' ' • 



(176) 



In particular, the first order term in has vanished; even though the kinematic phases scale with 
the square of neutrino mass, the CP-violating linear combination of phases scales with the sixth 
power of neutrino mass. 

This phenomenon holds even for higher numbers of generations: the CP-violating term (|134[) is 



4 EE^ sin2 ^ = 



i j>i 



2 2\ 

[ m, ; — rrij J 



(177) 



t j>i 



03 



3>i 
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The double sum in equation (177) X/j>i contains all i < j combinations, as does the double sum 



J2j Si<ji so we ma y substitute the latter for the former in the second term of equation (177). We 
may then exchange the dummy indices i,j in the second term to obtain an overall sum which is 
equivalent to the sum in equation (|165|) : 



5>?£*«-£I>?*« - E^E^-EE*^- 

i j>i i j>i i j>i j i<j 

= E™?£^-££™^ ( 178 ) 

i j>i i j<i 

= £™?£%+££^ 

i j>i i j<i 

Thus, for any number of generations, the CP-violating contribution to neutrino oscillations when 
the phases are small becomes on the order of ($y) 3 , which is the order of (Amy) and therefore 
negligible for small m. 5 



In Chapter VI we show that reasonable values for the mass-squared differences yield phases 
too small to produce a CP-violating effect observable at current accelerator experiments. Longer 
baseline experiments would, of course, allow larger phases, so CP violation too small to be detected 
at accelerator sites could possibly be detected by long-baseline experiments unless matter effects 
conspire to reduce its effect [[l0| . We will discuss below the possibility of detecting the effects of CP 
violations, even when CP violation is not directly observable. 



The real parts of the constraints (152) and (154) also lead to interesting conclusions. Consider 



first the homogeneous constraints for which the Kronecker delta is zero. These relationships give 
the singly-degenerate boxes as sums of ordered boxes: 



\V Xi \ 2 \V Xj \ 2 = Xi a Xj = - E * D Ai = " £ Mi, and (180) 



5 The manipulation of the double sum performed here can also be used on the exact CP- violating term (134). The 
result is 

4 £ £ sin 2 ^ = 8 £ sin £ cos ^ • (179) 

i— 1 j>i i= 1 
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\V ax \ 2 \Vp x \ 



J 0X 



/3x 



(181) 



For three generations, each of the sums contains two terms, so we now have expressions for our 
singly-degenerate boxes in terms of two nondegenerate boxes measurable in neutrino oscillations. 

For fixed in equation ( |l80| ), A can take n possible values, implying n constraint equations. 
N ordered nondegenerate boxes appear in these n equations. Thus, for N < n, which is true for 



< 3, the unitarity constraint (180) and (181) may be inverted to find a nondegenerate box in 



terms of singly-degenerate boxes. Adding the sums in equation (181) for x — i and x — j, then 
subtracting the x = k sum yields the expression: 



%i = ~\ {\Vai\ 2 \Vpi\ 2 + \v aj \ 2 \v Pj \ 2 - \v ak \ 2 \v^ 



(182) 



where we have used the property that Rs are unchanged under mass (or flavor) index interchanges. 

Many sources, such as reference fl3|l , use these sums as the coefficients of the oscillatory terms 
in oscillation probabilities. In terms of the |Vai| 2 |V/3i| 2 for three generations, 



P (x) = 2(sin 2 $ 12 +sin 2 $ 13 -sin 2 $ 23 ) |Ki| 2 |Vsi 



2 (sin 2 $ 12 - sin 2 $ 13 + sin 2 $ 23 ) \V a2 \ 2 \Vp 2 \ 2 
2 (- sin 2 $ 12 + sin 2 $ 13 + sin 2 $ 23 ) | V^ 3 1 2 1 Vs 3 | S 



(183) 



Such an expression does not include the possibility of CP- violation. Our boxes are clearly preferable 
for their compactness and versatility. 



Similar manipulation of equation (181) gives an alternate expression in term of the cyclic triad 
(a,/?, 7): 

%i = -\ {\V ai ?\V aj \ 2 + |^| 2 |%| 2 " \V 7i \ 2 \V 7 j\ 2 ) . (184) 



The unitarity constraints (180) and (181) greatly simplify our expressions for a doubly-degenerate 
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box: 



1 — 1 c 



1 — 1 f 



ax n 
1 — 'rvfl* 



<— >r 



xj^y j^i, and (185) 



(186) 



For three generations, doubly-degenerate boxes are expressible in terms of the real parts of six 



ordered boxes, rather than the nine complex boxes used in equations (143) and (144). For example, 
in three generations 



I V r _ 



|4 _ 11 



□i 



($22 + $3 2 )($23 + $. 



:s:s 



t23 + ^33 



(187) 



Summing equation (180) over j ^ i yields another expression for |V>u| 2 in terms of nondegenerate 
boxes: 

ivxii 2 X)i^i a = i^i 2 (i- i^i 2 ) = -EE ( 188 ) 



The same result is obtained by summing equation (152) over all j. The explicit solution, valid for 
any number of generations, of the above equation has a two-fold ambiguity: 



l± j 1+4 EE^- 



(189) 



This solution gives the doubly-degenerate box in terms of the (n — l) 2 nondegenerate boxes of the 



double sum. When the double sum is small, an approximation to the exact equation (189) is 



i-EE^,i+EE^ 



(190) 



For three generations, this yields | V>i | 2 as a linear equation of four nondegenerate boxes. For 
example, 



l^l 2 = 5 



1± v/l +4(^22+ ^32 + ^23+ H33) 
( $22 + $32 + $23 + $33) , 1 + ( $22 + $32 + $23 + $3s)) 
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We may use the homogeneous unitarity conditions (18C) and (181), along with the tautologies 



( |138| ) and ( |139| ) to obtain constraints between ordered boxes, thereby reducing the number of real 



degrees of freedom. Recognizing that the tautologies (138) and (|139| ) give 



= Re 



°^vj j j + °J 7)3 °J Xj ^A j — °$y j ^A j ^Xj + A 



and 



'□ 



Pi 



= Re 



(° R |3») + ("*/3w) 



(191) 



(192) 



our unitarity constraints fll80D and (181) become 



m m , ad iu m cwi m rij . cm, oti rjd _ 

ttrjj K>\j ^Aj T" •ir/j J A j J^Aj ^Xj J Aj ~T~ -^77 j J A j JAj — 

-((^) 2 + (^) 2 )^^, 



and 



(193) 



Ttfla: Ttflu Rfiv + Tfe V)/ 



/3a; tl/3y tLfSy + J fix J fly tlf3y ~ °Sf}x rif3y Jf3y + °&f3x °$f3y Jf3y 



(194) 



with the usual inequalities r\ 7^ A 7^ a, i 7^ j, in equations (|191|) and (|193|), and a ^ [3, and x =/= y =/= i 



in equations (|192j) and (191) satisfied 



We may rename some indices and solve both equations (|193|) and (194) for one 



%i %Aj ~ 'RAj »A,) + ((^A,) 2 + (»A,) 2 ) gr 



%y %y ~ %y %y) + {(%y) + ( jfr) j Ez^J ^ 

^ ( Wv + %y) + °&Py ( ^Py + fyy) 



(195) 
(196) 



Similar expressions could be found for °fpj if so desired. 

These constraints interrelate imaginary and real parts of boxes through unitarity for any number 
of generations. If CP is a good symmetry, all of the vanish, and the above expressions for 
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^i/Sj will be different than if CP were not conserved. For three generations, t must equal A, and 



fpj = — f\j = W\j by equation ( |165|) . Similarly, z = y, and — — fp y — "fpy by equation ( |166[ ) 



Equations (195) and (|19q) for n — 3 become 



" ^ + ^ (197) 

Rearranging the equations gives a measure of CP-violation in three generations: 

( y,,f = « w & Xj + 1< 'K v . + H Xj & Xj (199) 

= °&pj "fifty + fifty + °&[3y < R/3j / . (200) 

These relationships between Rs and Js in three generations exhibit a simple parameter symmetry. 



Cfpj) 2 in equation ( |199| ) equals l^- M/3j + terms cyclic in (a, /3, A); in equation (200), it equals 
"Jfysj °^/3i+ terms cyclic in (i,j,y). 

If CP is violated, 'fpj ^ 0, so the combination of real parts on the right-hand sides of equa- 



tions (199) and (EOOl), measurable with CP-conserving averaged neutrino oscillations (discussed in 



detail in Chapter VI), cannot be zero. Thus, even if CP violation is not directly observable in an 
experiment because of a short baseline or small mass differences, the effects of CP violation may be 
seen through the relationships among the real parts of different boxes! 

We may infer similar constraints by setting the imaginary part of a singly-degenerate box to 
zero: 



Cf% ati 7&i m Qf3 Cf% ogi rki cm cha tvl 

n v j o\j tX\j — Jr/j i^Xj — *x v j tx\j J \j — J v j J Xj J Xj _ q 



and 



'"' ( "' n : - ) = '"' ( n llu^ V ) = 
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Qti Ota 



■fix J By K/)y ~ ™$Bx °&Py "H-By ~ ~H-8x ~K-8y "Jf3y ~ J0x ~->~3y "J By 
py fly 



ofi aft asp 



■ JBv Jfl 



SO 



ad a*i m ad a4 m ad ad ni ad ad ni r\ 

tirjj J\j tX\j J. n j tL\j H\j ttfjj J\j J v j J\j 3\j — U, 



(203) 



and 



ad ax ad ad a& ad 

Jfiv n-Bn — JBx rtrt,, Tlfl,, — Ti 



■fix Jf3y Kffy — Jf3 X iXjjy iXjjy — K,/3 X ttfiy Jfiy ~ J @ X J fjy J fly 



ax ad ad ax 

JBv — J Br. JB-ii JB 



(204) 



With a bit of index switching, these lead to constraints of the form 



3y °R/33/ - ^2/ ' 



(205) 



For three generations, unitarity enforces the simple relations (|165|) arid (|166|) among the Js, and 



we find that equation (205) reproduces the same expressions (|l97|) and (|l98[) which were derived from 



the unitarity constraints (195) and (U9q) above. For n > 3, however, equation (205) is different from 



the prior equations, since the sums Ylr^a p an( ^ ^2z^ij m equations ( 195 ) and (196) ensure that 



boxes with all n flavor or mass indices enter into equations (197) or (J198|), respectively. Equation 



(205) involves only three such indices, so it will be different from the unitarity constraints in higher 
generations, providing additional identities among the real parts of boxes. 



We have not yet used the inhomogeneous unitarity constraints (152) and (154) that occur when 
the Kronecker delta is not zero. The homogeneous unitarity constraints cannot provide the absolute 
normalization of the V a i or the boxes, so we must need the inhomogeneous unitarity constraints. 
These inhomogeneous constraints are functions of degenerate boxes: 6 



1 — 1 f 



■ = V mQ « . and 



(208) 



In terms of mixing-matrix elements, equations ( |20S[ ) and (209) are trivial. For example, equation (208) is 

2 



\V ai \ A + ^|V ai | 2 |K 7i | 2 = iv^l 2 

n 

Jj^| 2 = i. 



(206) 



(207) 
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■ 4- ai n — \l ai n - 

u qi i / J '—'az ■ — V u m 
z^i 



(209) 



Parenthetically, we note by comparing equations ( |208| ) and ( |209| ) that a sum over mass-degenerate 
boxes equals a sum over flavor-degenerate boxes: 



E °' u 



E ai n 
— OLZ 



(210) 



Equations (206) and (209) can be rewritten strictly in terms of nondegenerate boxes by using the 



homogeneous unitarity constraints (|180|) and ( 181 ) . (If these expressions were not available, we would 



need to use the tautologies ( |143| ) and ( |144|) for the doubly-degenerate boxes, and equations (13? 



and (13E) for the singly-degenerate boxes, resulting in an incredibly ugly mess (as compared to the 
merely unincredibly ugly mess derived below)). We find 



Off 



ars 



E 



E ^ l 

z^i 



(211) 



with x =/= y =/= i, and 



E -£**) =°> 



(212) 



with A ^ tj ^ a. 

In three generations, each of the sums has only two terms. Eliminating the square root, and 
multiplying through by the denominator, we find constraints of the form 



(213) 
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and 



( °&\x ^Ay) ( ^r}x + ^^y) ( ^77X + *"^73/) [l + 2 ^ IIax + ^Ay + ^t^x + ^r/y)] + 
( °&\x + ^Ay) ( °&r)x + *^)7y) H~ ( ^Ax ~h ^Ay + ^772; + ^Tjy) ( ^?7^ + ^73/ ) — 0, 

where x ^ y ^ i, and A 7^ 7} ^= a. So, for a = 2 and z = 1, we find 



(214) 



&22 + I32) ( &23 + %33) ( $23 + ^33) [l + 2 ( & 2 2 + %32 + &23 + %33 J 
( &22 + %32) 2 ( &23 + ^33) 2 + ( ^22 + %2 + %3 + %3 3 ) 2 ( %23 + ^33) ^ = 0, 



(215) 



and 



( &22 + % 3 ) ( %2 + I33) ( $32 + ^33) [1 + 2 ( %2 + &23 + %2 + %a)] + 
( % 2 + ^23) 2 ( ^32 + % 33 ) 2 + (%2+ %3 + %2 + %s) 2 ( ^32 + &3 3 ) 2 = 0. 



(216) 



As with the constraints on the imaginary parts of boxes, many of the constraints ( 195 ) to ( 214 ) 
are redundant, but an independent set may be used to reduce the number of parameters significantly. 
The homogeneous constraints are much simpler than the inhomogeneous constraints, so we want to 
use as many of those as are possible. For the three-generation case, we may eliminate ¥£23 by either 
equation ( |l97|) or equation fll9S| ) 

10 _ $33 $33 ~ J 2 _ $22 $23 - J 2 / 01 „x 

^ 3 - -S33-H33 " - $23 - $22 ' ^ 

We may similarly eliminate ^,32 and $33 



u $22 $32 — J 2 _ $33 $33 — J 2 /„, x 

1132 ~ — S h — - — w h — ' y lV6 > 

— Tt32 — T122 — K33 — ^33 



and 



24 $23 $33 — J 2 _ $32 $33 — J 2 /II OA 

^ 3 " - £33 " »23 " -^33-%32' (219) 



S33 may be eliminated from equation (217), and $33 may be eliminated from equation (|219|) 



« _ $33 $22 $23 + J 2 ( $23 + $22 " $ 33 ) , 

^ " $33 $23 + $33 $22 = $22 $23 + ^ 2 ' ^ ^ 
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ft; 



32 



V33 



ft 



23 J 



ftj 3 %32 + ^33 &23 - %32 ^33 + J 2 



(221) 



Equation (218) will not provide an additional constraint; it is redundant to the other two. We must 



turn to equations (213) and (214) to eliminate the last degree of freedom. This necessity is expected, 
since without this set of inhomogeneous constraints, the rows/columns of V will not be normalized 



and our set of boxes will not yield a unitary V. We will here choose the constraint ( 215 ), with a = 1 
and i = l since its expression in terms of our four remaining boxes is the least complicated. We may 



substitute the second equality from equation ( |219D for ft^ and the second equality of equation (217) 
for %23, leaving a constraint which contains the four parameters ft 22 , ft23, ^32, and ^33: 



( ft 22 + ft 32 ) 
+ I32 + ft 23 



ft 23 + 



^32 ^33 — J 2 \ ( ^-22 ft23 — J 



"^33- 
%32 ^33 — J 2 



ft 



22 



^33" 
■ ^32 



'32 



ft 



23 



^32 
+ (ft 22 
^32 ^33 



&23 ~ &22 
> 2 



% 
J 



32) 
2N 2 



ft 



I33 J [1 + 2 ( ft 22 

ft 3 2«33 ^ 2X2 



J 2 



23 



^33 



% 



32 



-«33- 
ft22 ft-23 — J 2 



32 



ft 



23 



ft 2 



33 



(222) 



Multiplying through to place all of the terms in the numerator, we are left with 



= ( ft 22 + ft 23 ) 2 ( ft-22 + ft.32) 2 (- ft.32 % 3 + ft.23 ( %2 + I33) + J 2 f 
%2) 2 + (%3 + ft 22 ) (%2 + I33) + J')" 

x(%3%3+%2 (-%3+% 3 ) J')" 
( ft 2 2 + 1 32 ) ( ft 22 + ft 23 ) ( ft 2 2 ( ft 23 - %3) ~ %3 133 - J 2 ) 

x ( % 2 + 2 ft 22 I32 + 2 ( ft 32 ) 2 + 2 132 ft 23 

+ I33 + 2 ft 22 %3 + 2 ft-23 ^33 + 2 J 2 ) . 



(223) 



This equation is quartic in II33; we may eliminate this box by solving the equation 



A + B(% 3 )+C(% 3 ) 2 + D(% 3 f + E(% 3 ) 4 = 0, 



(224) 
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with 



A = 



( %2) 6 (1 32 ) 2 - 2 ( % 2 ) 5 % 2 (2l 3 2 ^23 + J") 
( ^ 22 ) 4 [9 ( 1 32 ) 2 ( % 3 ) 2 + 10 1 3 2 % 3 J 2 + J' 
( %2) 3 [6 (1 32 ) 3 (^2 3 ) 2 + (1 32 ) 2 (^2 3 ) 2 (1 + 10^23) 

+ %2 ( %s + 20 ( ^ 23 ) 2 - 2 J 2 ) J 2 + 6 % 3 J* 
+ %2 { ( 1 32 ) 2 (1 + 10 % 3 ) J A + 4 % 3 J* (- ( % 3 ) 2 + J 2 ) 
-2(% 2 ) i [(% 3 ) 3 -2% 3 J 2 

+ %2 % 3 J 2 (-2 ( £ 23 ) 3 + J 2 + 6 % 3 J 2 ) 

+ ( I32) 3 £23 (- ( % 3 f - 2 ( £ 23 ) 3 + J 2 + 8 % 3 J 2 )' 

- ( ft 22 ) 2 [3 ( 1 32 ) 4 ( % 3 f + ( 1 32 ) 3 % 3 ( ^23 + 8 ( tt 23 ) 2 - 6 J 2 ] 

+ (9 ( £ 23 ) 2 - 2 J 2 ) J* + % 2 J 2 [( ft 23 ) 2 + 14 ( ^ 23 ) 3 - J 2 8 % 3 J< 
+ ( I32) 2 [( % 23 ) 3 + 4 ( tt 23 ) 4 + 2 ( % 3 ) 2 J 2 3 J 4 ] } 
+ J 2 [( 1 32 ) 3 ( % 3 f (1 + 2 £ 23 ) + ( 1 32 ) 4 (2 ( % 3 f J 2 ) 
+ ( I32) 2 ( I23 + 3 ( 1 23 ) 2 - 2J 2 ) J 2 J 2 (( % 3 ) 4 + J 4 )" 



(225) 



D 



-{%2+ % 3 ) {2 ( £ 22 ) 5 1 32 - ( 1 32 ) 4 (( % 3 f 2J 2 ) 

(l32) 3 (l + 2tt 23 )((^ 2 3) 2 -J 2 )+4^ 6 

- 1 32 ^ 2 (2 ( £23) 3 + J 2 ) + ( £22) 4 (-2 ( I32) 2 + 6 1 32 £23 + 2J- 
+ (%2) 2 J 2 (-3% 3 -6(% 3 ) 2 +8J 2 ) 

- 4 ( ^ 22 ) 3 [( 1 32 ) 2 % 3 - 2A 23 J 2 + % 2 (-3 % 3 + J 2 ) 



22 



6 ( 1 32 ) 4 % 3 + 2 ( I32) 3 ( % 3 + 5 ( % 23 ) 2 - 3J : 



+ 2 1 32 ( ^23 + 9 ( ^ 23 ) 2 - 2 J 2 ) J 2 + J 2 (-2 ( 1 23 ) 3 + J 2 + 12 % 23 J 2 ) 

+ (%32) 2 tt23(-tt23 + 2(tt 23 ) 2 + 18 t 7 2 )" 

- (1 22 ) 2 [l2(l 32 ) 3 ^ 23 + 2(l 32 ) 2 ^ 2 3(l + 4^ 2 3) 



+ J 2 



% 3 ~12(% 3 ) 2 +4J 2 



(226) 
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%2 (-2 ( 1 23 ) 2 - 8 ( % 3 f + J 2 + 2A % 3 J 2 )] } 



C = - 



V23J 



( %2 + %s) { ( ^ 22 ) 5 + 3 ( % 2 y % 3 + ( % 2 f % 3 (1 + 2% 
( % 2 ) 4 (-2 % 2 + % 3 ) - 2 % 2 % 3 J 2 + 6 % 3 J* 
( % 2 ) 3 [4 ( % 2 ) 2 - 16 ! 32 %3 + 6 ( I32) 2 - 2 J 2 
( 1 32 ) 2 [-2 ( 1 23 ) 2 - 2 ( 1 23 ) 3 + J 2 + 10 % 3 J 2 

%2 [3(! 32 ) 4 + (! 32 ) 3 (i + 8£ 23 ) 

2 ( 1 32 ) 2 ( % 3 + 4 ( 1 32 ) 2 - 5 J 2 ) - 2 ( % 3 + 6 ( % 3 ) 2 - 3 J 2 ) J' 



6(1 



.">2 J 



; 1 22 ) 2 (1 - 2 Baa) + ( 1 23 ) 2 + 4 ( 1 23 ) 3 - J 2 



14 l 23l 7 2 + %2 (-3 £ 23 - 22 ( 1 23 ) 2 + 8.J 2 )] } 



(227) 



D = -(1 22 +1 23 ) 2 {4(1 22 ) 2 ( 1 32 -1 23 ) 

- %2 %3 + ( 1 32 ) 2 (1 + 4 % 3 ) + 4 % 3 J 2 



%2 + 4 (1 32 ) 2 - % 3 + 4l 32 1 23 - 4 ( 1 23 ) 2 + 4J 2 ] } , and 



(228) 



E 



tt22 + TK-23 I 



(229) 



Eliminating with this constraint leaves us with three real parameters ¥l 22 , : ft 23 , and 2 tl 32 , and 
the one imaginary component J as our basis. This final constraint relating the four boxes, while 



not as ugly as it would be without the use of equations (18C) and (181), is quite complicated. The 
algebraic solutions (which are obtainable) are too messy to be particularly illuminating, so we will 
choose to eliminate our final spurious parameter by solving equation (|215|) numerically rather than 



algebraically when we compare to experiment in Chapter VI 
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4.3 More Unmeasurables 

The previous Sections have developed relationships between the lepton mixing matrix and the 
boxes. Now we turn our attention to relationships between the neutrino mass matrix and the boxes. 



We showed in Section 2.7 that we cannot distinguish experimentally between theories containing 
a diagonal Mi with U VL — V and theories containing mixing in both the neutrino sector and the 
charged lepton sector. We thus lose no physical information by restricting the mixing to the neutrino 
sector. Most models for fundamental lepton mass matrices, however, do not make this ansatz of 
a diagonal mass matrix for the charged leptons. On the contrary, many models postulate as a 
symmetry principle the same initial form for the charged lepton mass matrix as for the neutrino 
mass matrix. Such models will yield predictions equivalent to those of many other models, including 
a model that contains all of the lepton mixing in the neutrino sector. 

Consider a model for lepton mass matrices that contains a neutrino mass matrix M„ and a 
lepton mass matrix Mi . We will first assume a Dirac mass term for the neutrinos for simplicity, then 
extend our treatment to the general case. The Lagrangian in the mass basis (Q) which we derived 



in Section 2.7 contains the terms 



U^B v N R m + l" L l D t l% + c 2 W+UL m V- 1 1 n r £ + h.c. (230) 



We illustrated in that Section how the Lagrangian is invariant under a rotation of the mixing matrices 
Uv L ,i L — ► RU Ul j l by any unitary matrix R. Such a transformation is equivalent to rotating the 
mass matrices M v and Mi by the same arbitrary unitary matrix. In particular, we may choose this 
arbitrary rotation such that one of the mass matrices is diagonal from the beginning. This implies 
the mixing matrices U VL = V and Ui L which we assumed when we derived our box relationships. 
The charged lepton mass matrix is diagonal in this basis, which we will call the standard basis. The 
neutrino mass matrix M v in this basis is related to the original mass matrix M v by M v — Uj~^~M v , 
which can be seen by equating the diagonal mass matrices in the two bases: 

D„ = U^M v U Nr = V- x M v U Nr = U-^U 1l MU Nr , so (231) 
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M v = VU- 1 M U = U£U VL U£M V = U^M V . 



(232) 



When more generations of neutrinos are introduced, the left-handed neutrino mixing matrix U VL 
no longer has the same dimensions as the left-handed charged-lepton mixing matrix. We will use 



the general block form of the neutrino mass matrix from Section 2.5 



M v = 



M 7 



Ml 



M T 



M D M vx 
M s M Nx 

M x ) 



The mass term (B2I) containing this matrix is of the Majorana form, so the matrix M v must be 



symmetric and may be diagonalized by U v X M V (U^) , as described in Section 2.1 . We may 
combine the spinors of all three types of neutrinos into one n-dimensional vector i\> v : 



N 

V x J 



(233) 



In addition, we will sometimes break the neutrino rotation matrix U v into three submatrices, U\, 
U2, and U3 such that 



u 2 



(234) 



U\ is «i x n, U2 is tir x to, and U3 is n x x n, with n = til + n R + n x- Such a decomposition is 
necessary because the charged current term contains only vl, which is related to the n mass states 
contained in i/j™ by vl = U^i/j^. These Ui are not square and so are not unitarity, but UiU\ =U( 
here has dimensions Ul x Ul), so W is the right inverse of Ui. 



Using these definitions, the mass-basis Lagrangian (230) is 



il> V LU v U- x M v Tf*- 1 U*il) U R + IlU^U^M^UI-Hr + VZUxU^Ui^Hl + h.c. 
= WzPvWr + W d Jr + WlV'^l + h - c - 



(235) 
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The matrix in the charged current term V — Uf^Ui is not square, but has dimension x n. 

We still have the freedom to incorporate all of the lepton mixing in the neutrino sector, effectively 
replacing U± with V, so the neutrino mixing matrix in this standard basis is 



U u = 



u 2 

\U 3 J 



( rr-l 



ur L o o 



\ 



1 
1 



(236) 



This replacement implies a different non-diagonalized neutrino mass matrix M„: 



( ur 1 0^ 



1 
1 



1 uj- 1 ^ 



UuDUt 



\ 



\ 




1 
1 



^ f uT- 1 o o 



1 
1 

\ 



V 



1 
1 



(237) 



This transformation includes both a right and left rotation because it involves Majorana neutrinos. 
The transformation for til generations developed in equations (|23l| ) and ( [232| ) assumes the right- 
handed fields Nr are independent of the left-handed fields vl , so the right-handed rotation matrix 



Un r appearing in equation (231) is unaffected by our new choice of U VL . For the Majorana mass 
term indicated in equation (|55|), this is no longer true. The right-handed fields are related to the 
left-handed fields, and the right-handed rotation matrix must equal (C^J) . So changing to the 
standard basis rotates the mass matrix on the right as well on the left. Such an additional rotation 
could be performed in the three-generation Dirac case for consistency's sake (since we cannot measure 
Un r in the absence of a right-handed current), but it is not necessary on mathematical grounds. 



We note that under the transformation in equation (237) the Dirac subblocks Md and Mp, in the 



matrix (53) do transform to Ui L Mu and M^Uj ', in agreement with equations ( |23lD and (232) 



The issue here will become moot in the next Section when we restrict ourselves to the hermitian 
mass-squared matrix MMK The extra rotation cancels in this product. 
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Even the diagonal mass matrix D is not completely measurable. The mass matrix M does not 
have to be hermitian, so the rrii in D are not necessarily real. We have the freedom to perform a 
chiral rotation e 1 * 75 on the fermion field definitions; under -0 — > e iei5 tp, mass terms become 

ipEMipn + h.c. -> ^Me 2, ^ 5 fe + /i.e. = (238) 
ipi,M (cos 20 + 175 sin 29) ^ R + ifi R M^ (cos 29 + ijs sin 29) tpL, (239) 

where we have used the properties 75 =4 and 75 2fc+1 = 75 for integer k in the Taylor expansions of 
e 21015 , sin0, and cos#. Since j^Pl = —Pl, an d 75 -Pr = Pr, we have 

~ L Mi 1 ^ R + ^Aft*76Vi = t^MiipR - i^MUip L , (240) 

so our whole mass term has become 

JlMe 2i6 il)R + V^MV^Vl = lpZe 2ie Mi}R + h.c. (241) 

A chiral rotation by 9 = ? has the same effect as the transformation M — > — M; a rotation by 
6 = j effects M — > iM. We can separately perform chiral rotations on particular flavor fields, so the 
phases of individual mass eigenvalues are arbitrary. These rotations lead to corresponding effects 
on the diagonal D; the phase of the mass matrix, and the phases of the fermion masses contained 
in the diagonal mass matrix D, cannot be measurable since we may absorb these phases into the 
fermion field definitions. 

In the next Section, we will derive equations expressing the elements of the matrix M„M' in 
terms of the observable boxes. To use these relationships to place constraints on the elements of the 



mass matrices given by particular models, the transformation (232) must be invoked. Once extra 
neutrinos are included, the CKM matrix V — U± appearing in the charged-current term is no longer 
equivalent to the full neutrino mixing matrix U v . The matrix V knows about all n of the neutrino 
mass states, since it connects vl with ip™, but it knows nothing of the sterile flavor neutrinos and 
therefore cannot be used to describe mixing between active and sterile flavors. The entire matrix 
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U v is necessary to completely describe neutrino mixing. In what follows, however, we will use the 
symbol V to denote the left-handed neutrino mixing matrix, as we did in the previous Sections of 
this Chapter. For three generations, this assignment presents no ambiguity; for n > 3, the neutrino 
mixing matrix V will have different dimensions than the CKM matrix V , so again there should be 
no confusion. 



4.4 More Measurables 



We showed in Section 4.1 that 2n — 1 relative phases can not be observed. For Dirac neutrinos, 
we also have the freedom to rotate Un r ad nauseum since it is not observable in the absence of 
right-handed currents. The effects of this last degree of freedom may be ignored if we consider 
only the hermitian product MM^ rather than the generally non-hermitian mass matrix M. Should 
right-handed currents be found, one may explore Um r via the hermitian matrix M^M. Note that 
in the absence of the symmetry imposed by Majorana neutrinos, M'M = M^M, since the rotation 
by U VL in M — U~^MUn r is canceled; this matrix is diagonalized by Un r alone into D 2 : 

U\,M^MU Nr = D 2 . (242) 



The phase ambiguity of the mixing matrix requires that any observables which are functions of 
the left-handed neutrino mixing matrix V and/or the product MM' must be invariant under the 
simultaneous transformations 

V -> YVX, and MM^ -> X^MM^X, (243) 

where X and Y are diagonal matrices of phases. These invariances are identical to those possesed by 
traces of products of MM' and diagonal basis matrices, or of traces of products of V and diagonal 
basis matrices. Such traces offer an elegant method of obtaining the measurable functions of the 
mixing matrix and the mass matrix. 
We introduce the notation 

H = MM\ (244) 
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and the diagonal basis matrices 



(245) 



We may express mass-matrix invariants in terms of these newly-defined quantities as 



Tr 



(246) 



Thus, 



Tr [HE a ] = H aa 



(247) 



is an invariant, as is 



Tr [HE a HEp] = \H a p\ 



(248) 



Higher-order traces such as Tr \HE a HEpHE^\ are also invariants. 

Turning to the mixing matrix, we see that equal numbers of Vs and V's provide invariants. V 
placed into a trace by itself is not invariant, since the phase matrices X and Y occuring on cither 
side of the rotated mixing matrix do not cancel in the trace: 



Tr [VEt] Tr [YVEiX] . (249) 

Here we have used the commutativity of the diagonal matrices X and Ei. The product V'V, 
however, is rotated by the single transformation X, which will cancel in the trace, so 

Tr [V*E a VEi] = Tr [X- 1 V*E a Y- 1 YVE i X] = \V m \ 2 (250) 

is invariant |35|, as is the fourth-order trace 

Tr [VE^EiVEpV^Ej] = V ai V^V Pj v} a) (251) 

which is a box. Again, traces of products higher-order in even powers of V are also measurable. 
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The mass-matrix invariant traces are related to the mixing-matrix invariant traces by 



H = VD 2 V^ = V 



E^ m i 



(252) 



So wc find 



tf QQ - Tr[ff£j - J2 m P T [VEjV^E a ] = ^2m 2 j \V aj \ 2 = E m V " J ' D 



(253) 



and 



|ff Q/3 | 2 = Tr [/f £ Q /f ^ 



1=1 3=1 



i=l 3=1 



The next-order expression becomes 



n n n 



H aP H Pl H ia = Tr [HEpHEyHE a ] = E E E rf m Wk v * v V V M V ft V r'° v L- ( 255 ) 

i=i j=i fc=i 

These products of six mixing matrix elements and even the product of eight mixing matrix elements 
could be measured in processes such as the decay fi^ — ► e ± e + e~ which proceeds through penguin 
and box graphs with internal neutrino lines. 



The right-hand side of equation (254) may be written as 



Pi 



EE 

»=1 j^i 



WTjTTl 2 m U 

1 J 



03 



1=1 



mf|y m | 2 |^| 2 + 2^^ 

i— 1 j>i 



TO 2 TO 2 ' 



(256) 



From this we see that not only is information about relative phases lost, but that the measurable Js 
do not enter either. The missing information is carried in the higher-order observables, such as those 



found in equation (255). An alternative approach is suggested by the unitarity constraints (|195|) 



and (196) which may be rearranged to express in terms of Rs. Unitarity therefore provides a 
mechanism for implicitly obtaining the Js from real trace relations. 



Equations (|253j) and (254) involve the neutrino mass matrix in the standard basis, with all of 
the mixing in the neutrino sector. If we want to compare the predictions of mass models specified 
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in a non-standard basis to the measurable quantities in the right-hand side of equations ([2 5 3D and 



( |254| ) , we need to rotate M back to the M v of the specific model as described in equation (|232j ) 
Performing this rotation on H a p, we find 



a/3 1 



(257) 



and M„ are given by the model, while the right-hand side of equation (254) (which equals (257)) 
contains measurable quantities. Finding individual elements of M„M^ in terms of the observables 
might be messy, but it is possible using the above equations. 

Consider next the converse issue of constructing the boxes from a given M (and therefore a given 
H). Equations (|253| ) and ( |254|) present equations for each of the n 2 possible a, (3. These equations 



involve ./V 2 ~ n 4 ordered boxes. Equations (|253|) and (254) therefore cannot be simply inverted to 



obtain the boxes from the mass matrix in the general case. |41 . The unitarity constraints presented 
in the earlier Sections of this Chapter reduce the number of independent boxes so the inversion 
could in principle be done. But the relations would be highly non-linear and intractable. To obtain 
the boxes from a given mass matrix, we will diagonalize the specified matrix and solve for the 
eigenvectors. These eigenvectors may then be used to construct the mixing matrix V, which gives 
the boxes. This procedure is illustrated in Chapter |v| for a few selected mass matrix textures. 
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CHAPTER V 

APPLICATIONS OF THE BOXES TO SPECIFIC PHENOMENOLOGICAL MODELS 

Having developed a parameterization for neutrino oscillation which is independent of neutrino 
mass and mixing models, we will now examine how different models may be expressed in terms of 
the boxes. 



5.1 Models of the Mixing Matrix 

5.1.1 The Standard KM Parameterization 
For Dirac neutrinos, the mixing matrix V becomes the familiar CKM matrix from quark mixing 
(although the mixing angles need not have the same values as for the quark case), and V may be 
characterized by three mixing angles 9 a and a phase S. A popular choice of V is the Kobayashi- 
Maskawa matrix of equation (Bq): 



/ 



Ci S1C3 S1S3 

-sic 2 C1C2C3 - s 2 s 3 e iS C1C2S3 + s 2 c 3 e lS 



-sis 2 C1S2C3 + c 2 s 3 e lS C1S2S3 - c 2 c 3 e lS 



\ 



In terms of these angles, we find that J , the imaginary part of the boxes is 



J = cisfc 2 s 2 c 3 s 3 sin 5 

= - sin 6\ sin 26\ sin 26 2 sin 29 3 sin 5 = J sin 5, 



(258) 



as predicted in Chapter Ill's equation (99). Using this definition of J, we may express our boxes as 



22 



12 



□ 



2:1 



2:5 



s i c 3 s 3( c i c 2 — s 2) + J cos6(c 3 — s 3 ) + iJsinS 

2 2 2 2 _ j AS 
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21 1 — i 2 2 2/ 2 2 2\ , t r/ 2 2\ , • t • r 

n 32 = SiC 2 s 2 (c 1 c 3 — s 3 ) + J cos o(c 2 — s 2 ) + smo 

22 D 33 = c§s§ [slcl(st + 6c? + 2c? cos 25) -c?] 
J 



s 

21 m 2 2 2/22 



+ 3(1 + c 2 )(c^ - s^s 2 - cj) cosd + iJsind (259) 



n 33 = s 1 c 2 s 2 (ciS 3 — c 3 ) + J cosS(s 2 — c 2 ) — iJ sinS 



"□32 = -c 2 s 2 S 2 c 2 _ 



12 D 33 = s\c\s\{c\s\ - c 2 ) + Jcos(S(s 3 - c 3 ) - iJsinS 



"□33 = -c 2 s 2 s 2 s 2 + Je w 



In terms of the standard parameterization the observable probabilities, presented here as boxes, 
are big ugly messes. It is no wonder that approximations have been commonly used when comparing 
experimental results to the standard mixing angles and phase! 



In principle, one can solve equations (259) with S = for Oi, (9 2 , and 63 in terms of a basis of the 



real parts of any three independent nondegenerate boxes. Then substitutions into equations (259) 
give the remaining six nondegenerate real boxes in terms of these chosen three. 5 is then given in 
terms of these three boxes and the imaginary part ±^7 of any box from 



5 = sin" 1 ( J n ^ ) . (260) 



J {0i, 8 2 , 3 ) 



In practice, we have been unable to analytically invert equation (25£) with only three basis boxes.Q 



Thus we have another reason to eschew the arbitrary unitary parameterization of V and work directly 
with the box algebra. 



1 Allowing four nondegenerate boxes as a dependent basis does provide an analytic inversion of equation (259). 
We find 

sin 2 26»i = -4( U rj 2 2 + U n23 + "^32 + "ass ) , (261) 

tan 2 6» 2 = „ — — n — , and 262) 
"□ 22 + n n 2 3 

tan 2 6> 3 = „ 3 33 ■ 263 

We see, however, no way to re duce this basis to three independent boxes. We choose the four basis boxes by 



examination of the relationships (259); the boxes which are most simply expressed in terms of the mixing angles are 
the boxes we choose for our basis. Thus the choice of which four boxes to isolate as the basis depends on the unitary 
parameterization used for V. V is constructed by three successive rotations and a phase matrix; one may choose the 
rotation axes (i,j,k) in V(8 = 0) = R~ i (9i)Rj(02)Rj.{O3) such that each of the desired four basis boxes has only two 
terms, being composed from three of the simplest Vs and one binomial V. 
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5.1.2 The Wolfenstein Parameterization 
Wolfenstein suggested [Q a now-popular parameterization of the CKM matrix based on data 
from the quark sector. The parameter A is defined by A = V12, and the rest of the matrix is 
determined by comparing matrix elements with powers of this small angle ^2|. (For example, the 
element V23 is on the order of V 2 2 in the quark sector, so it is expressed in terms of A 2 . Keeping 
the matrix unitary to the desired order of A yields the higher-order terms in V\\ and V 22 . The 
Wolfenstein mixing matrix to 0(A 3 ) is 



V 



1 



A X 3 A(p - ir,) 



, X 3 A(1 - p - ir,) -\ : 



A 



X 2 A 
1 



(264) 



This matrix is only unitary to order A 3 , so any relationships between our boxes based on the unitarity 
of the mixing matrix, such as the relationships proving the equality of the imaginary parts of boxes, 
will hold only to order A 3 . Wolfenstein describes the mixing matrix in terms of 4 parameters: 
A, A, p, and r,. We choose the three boxes n D22 , 21|=l 32 , and n D23 as our basis. Combining the 



definitions of our boxes (125) with Wolfcnstein's mixing matrix (264), we find 



"□ 22 _ A 2 (1-±A 2 ) 2 



-A 2 , 



(265) 



21 



□ 32 = A 2 A 6 (1 - -A 2 )(l -p-irj)~ A 2 \ & {\ -p- ir]), and (266) 



ll D23 = -A 2 A 6 (1 - -A 2 )(l - p - ir,) ~ -A 2 X 6 (p - ir,). 



(267) 



11 D 2 2 is real, and the imaginary component of 21 D 32 has equal magnitude but opposite sign as the 
imaginary component of 11 ^23 , so the three boxes act as four independent parameters. (Again, the 
boxes do not appear to have the same imaginary parts because Wolfenstein has thrown out terms 
higher order in A. If we were to keep the leading-order imaginary term for all of the boxes, all tfpj 
would be ±A 2 \ 6 r,.) We can now find the rest of the boxes and express them in terms of our chosen 
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three boxes. 



33 ^A'X b (p-ir])(l- P -ir])^ — 



11 n 21 n 

u 23 u 32 

n* 21 n 
U 23 — U 32 



4 D 



23 



33 



A 2 A 6 (p — irj) 



23 



32 



L n 33 ~-v4 2 A 6 (l-p-*77)~- 21 n 



32 



(268) 
(269) 
(270) 



B, our matrix of boxes in (1135]) becomes the following for the Wolfenstein parameterization: 



I ii 



B 



□ 



22 



V 



LJ32 

21 Q 



11 n 

— u 23 



J 32 



32 



23 



J 32 



23 



ll n 21 n 
u 23 u 32 

LL n* - ^32 



(271) 



/ 



The dominant box in this approximation is n n22 • Thus the v e <-> i/p transition would be most 
probable, with its primary contribution coming from the Am 2 2 term. The transition <-> ^ T is 
next-probable, down by the order of A 2 from the v e <-> z/^ transistion, and it receives its dominant 
contribution from the Am^ term. The least likely transition according to the Wolfenstein model, 
down the order of A 4 from the v e <-> transition, is the v e «-> v T transition, and it receives equal 
contributions from all three mass differences. 



5.1.3 The One- Angle Approximation 
A related near-unitary parameterization fitting the quark data uses only one real parameter, 8, 
and a phase 6. Experimental determinations of the quark mixing matrix elements are consistent 
with the following form for V [p4j : 



/ 1 



V 



3 „iS 



-6 6 e 



9 3 e~ lS ^ 



1 J 



(272) 
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This approximation is a special case of the Wolfenstein parameterization, setting A = \p — irf\ = 1 
and results in the following matrix of boxes: 



with 



B 



( 11 



□22 

u 32 
. 21 Q 



:!2 



□ 



:52 



("□ 22 )' 
21 n » 

I 1 c 



J :*2 



.'52 



32 



( n ° 22 ) 3 ; 



(273) 



n n 22 = -6 2 , and 21 D 3 2 = ~9 6 e iS . 



(274) 



The mixing matrix in (272) contains only two independent parameters, but our matrix B seems 
to contain three: EI22 > which is real, ( ^32), and ( ¥32). But the three are related by 



21 n 32 I 2 



: 2 A 3 2) 2 + ( 2 r32) 2 = ( 11 n2 2 ) 6 



(275) 



So we are indeed left with two independent parameters. 

Because this one-angle approximation is a special case of the Wolfenstein parameterization, it 
produces the same predictions for the relative significance of transition probabilities: v e <-> is 
the most probable, and v e «-> v r is the least probable. Discerning between the general Wolfenstein 
parameterization and this extension would require measurements of probabilities at different lengths 
so boxes within a row may be measured individually. 



5.1.4 The Dominant Mass Scale Approximation 
A popular approximation for neutrino masses and mixing assumes to 3 >> rai,m2 



so 



23 ~ ^13 



4j> 



Under this approximation, 



Re (B) S 2 {<$>) = 



^ ^22 ¥^23 + ^23 ^ 
%32 ^33 + %33 
¥^32 ^33 + %33 



/ 9 

' sin $12 



(276) 



sin $ 



13 
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and 



Im(S) 5(2$) 



J22 J23 + J23 



2J 22 , 2J 

J32 J33 + J33 



^32 %3 + "^3; 



33 



sin(2$i 2 ) 
sin(2$i 3 ) 



( H \ 

J22 



'32 



V ^ J 



sin(2$ 12 ) = Jsin(2$ 12 ) 



I 1 \ 

1 



(277) 



The second column of Im(y6) vanishes because the sums of boxes are real by equations (165) and 



am- 

For small x, 



sin 2 $12 << sin 2 $ 



13 



(278) 



due to the mass hierarchy. Thus experiments with a small x will measure the factors multiplying 
the bigger phase $13, such as the sum 

^23+^23- (279) 



In Chapter VI, we apply a dominant-mass-scale approach to the oscillation data currently available; 



in that treatment, the LSND experiment measures the sum ( 279). This sum is equal to the degenerate 



box K23 by equation (181) 



Sources producing neutrinos with uncertainties in p or x large enough that all the oscillatory 
terms average will provide measurements of a sum of all three boxes in a row of B: 



^ ^22 + $23 + $23 ^ 



ReB(S 2 ($)) = - 



^32 + $33 + $33 
$32 + $33 + $ 



33 



1 ELil^l 2 l^l 2 ^ 

ELil^| 2 |^l 2 
V ELil^l 2 !^! 2 J 



(280) 



where the last equality comes from equation ( 157 ). This type of asymptotic measurement could 
not detect CP violation directly since the average of the sin2$y function is zero. Our analysis in 



Chapter VI applies this averaging to both solar neutrinos and atmospheric neutrinos. 
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CP violation is a very small effect in this dominant mass-scale approximation, even when the 



oscillatory terms do not vanish. As seen in equation (277), the CP-violating term depends on 

2 _ 2 

sin(2<I>i2), so the CP-violating effects for local neutrino sources are of order mi 2p " 2 x. Oscillatory 
behavior would eventually develop over longer baselines. 



5.2 Models of the Mass Matrix 

5.2.1 Using a Fritzsch Mass Matrix 
In the late 1970s, Fritzsch jisf proposed a form of the quark mass matrices which has proven to 
be consistent with experiment. Extending this model to the lepton sector means the mass matrices 
M v and Mi for neutrinos and charged leptons will both be Dirac matrices and have the Fritzsch 
form HH]: 






ae I?1 










be 1 ^ 





6e^ 4 





\ 



J 



Mi 



Ae^ 1 
Ae iE2 Be 1 - 3 
Be lS4 Ce Ss f 



and 



(281) 



(282) 



The parameters a, b, c, A, £?, and C may be taken to be real positive and will be determined by 
the lepton masses. The postulated zero entries are today known as "texture zeros." They can be 
motivated by invoking symmetries in the Higgs sector; such symmetries may be natural in grand 
unified theories. Each of the mass matrices originally has five independent phases, but we may 
rotate the mass matrix and absorb the phases in the lepton fields: 



WM V N R + IlMiIr + h.c. = VZX Vl X-^M u X Nr X^ 1 r N r + l L Xi L X^MiXi R X^l R + h.c 

= 1^£M U N R + f L Mil~L + h.c, with (283) 
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v~ h = X~}v u N R = XzINr, f L = Xr% and l R = X,H R . 



(284) 



As before, vl,Il : Nr, and Ir are the 1x3 vectors containing all three generations of spinors. 
Following Fritzsch's treatment in p5|, we choose 



X v 



( 1 

g-' i(?l-€3-S4+&) Q 



X 



n r 



-i(Si+fo-S 3 -£4+S 5 ) o 



V 



, and 






e -*«l~«4+«5) 

with Xi L and Xi R found by substituting E for £. Under this transformation 



M v -> M u = X~}M v X Nr = 



( a ^ 



a b 
b c 



, and 



' A » 



AOS 

s c 



These real symmetric matrices may be diagonalized by the orthogonal matrices R v and Rf. 

m e 

, and R^MiRi = 



R v 1 M V R V = 



mi 
-m 2 
m 3 J 



V 



m e 






















m 



/ 



(285) 



(286) 



(287) 



(288) 



(289) 



The nii are positive, and the minus signs in front of the second mass eigenstates occur because the 
determinants of the M are negative. Such a sign factor may be absorbed by a chiral redefinition 
of the appropriate lepton fields. Under ip — > e l ^ 15 ip = *7sV'i rnipip — » —rmpij}, as discussed in 



Section 4.3. This chiral rotation, however, would put the mass matrix into a non-Fritzsch form. The 
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issue of the sign of the fermion mass is obviated if we work with MM 1 ' as in the previous Chapter, 
but we have here chosen instead to parallel Fritzsch's development. 

These transformations R and X take us to the mass basis: i/L m = R^ 1 X~^vl 1 Nr" 1 = 



R v 1 X N 1 Nr, et cetera. As discussed in Section 2.7, the mixing matrix arises in the charged-current 
interaction: 



Jcc = WflL = WX Vh R v (R^X^X^Rt) R^Xf^h = v^V^l 



L • 



(290) 



Once R v and Ri are found from (D89|) ,the mixing matrix may be calculated from 



^10 ^ 



V = R v 1 X v ^Xi L Ri — R v 1 



e ia 
e ir > 



Ri. 



(291) 



The two phases a and rj contain all of the phases in X UL and Xi L : a — (^i — Si) — (£3 — S3) — (£4 — 
S4) + (^5 — S5), and 7/ = (£1 — Si) — (£4 — S4). One of these phases may be determined by measurement 
of the standard CKM phase S in ( |98| ) . The other only appears in real terms (see below) and therefore 
may be absorbed in the definition of the mixing angles in the standard parameterization, as shown 
in @. 

This extra phase occurs because Fritzsch postulates a set form for the mass matrices, thereby 
depriving us of some extra degrees of freedom. As discussed above, any measurable quantities may 
be described by the masses and the mixing matrix, V. This matrix is unitary and therefore may have 
at most six phases. Five relative phases may be absorbed in the definitions of the left-handed lcptons 
and right-handed neutrinos, leaving one phase in this observable matrix. Fritzsch, however, assumes 
a form for the mass matrices which yields ten phases. Only eight relative phases may be absorbed 
in the lepton fields (the left-handed neutrinos and charged leptons must be rotated together, so nine 
rotations may be performed on the left-handed lcptons, the right-handed neutrinos, and the right- 
handed charged leptons), leaving two phases in the mass matrix. Under this method, the rotation 
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matrices have the form 



X VL - X L 



gi(?l+«2-?3+e4-?5) o 

e i{i2+a) 

e *(?2-«3+2C4-«5+r)) 

\ 



1 

e j(€2-53+C4-?5) o 











and 



e i(6-H 2 ) o 

e i(«l--YU+«2-?3+?4-« 5 ) 











»(£a— S s ) 



(292) 



(293) 



(294) 



These rotations yield the matrices M; and X M^, with X equal to the net phase matrix above: 



X 



110 

e la 
V e l?) 



(295) 



Once we have simplified the mass matrices to M v and M; which have the same form, we need 



only solve for R v by diagonalizing M v in fl289| ) and then obtain i?; by substitution. The first step 
in diagonalizing a matrix is to solve the eigenvalue equation, a task which is simplified because we 
know the eigenvalues, mi, —m2, and 7713, and want only the unknown a, b, and c in terms of the 
masses: 



Det(A/ l , — ml) — —m + cm + (a +6 )m — a c = (mi — m)(—m2 — m){m^, — m). (296) 



Solving for the unknown a, b, and c by matching coefficients of m, we find 



C = mi — TO2 + TO3 
6 2 + a 2 = mi?7l2 + J712TO3 — TOaJTli 



(297) 
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a 2 c = m 1 m 2 m 3 . 



The first of these relations gives c explicitly, and we can solve for a and b: 



2 m\m,2mz 
a = 

mi — 777,2 + W3 
6 2 = 77717712 + 7712 7773 — 7773 7771 — a 2 . 



(298) 



Because a and & were defined to be positive real numbers, we must always take the positive square 



root of the right sides of (298). 

We may construct the diagonalizing matrix R u out of the eigenvectors of M„ . These eigenvec- 
tors satisfy the relationship 

M„ifj = m0, (299) 



which provides three constraints on the elements tp\, 02, and ip3 of 0: 



aip2 = mtpi, 
q/01 + 6-03 = 77702 , and 

602 + C0 3 = 77703 ■ 



(300) 



The first and the second of the equations in (BOOJ) describe the form of 0. The third, combined with 



the conclusions of the first two, reiterates the eigenvalue equation (296). We thus find that has 
the form 



= 77 



777 & 



(301) 



where 77 is a normalization constant, chosen to make the length of equal to 1: 



2i2 2 l2 1 / 2 2\2 

— — a + 777 +(777 — a ) . 



(302) 



These relations hold for m equal to each of the eigenvalues mi, — 777, 2 , or 7773, giving the apparently 
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simple form for R v 



R v 



n±ab 
nimib 



n2ab 



n 3 ab 
n 3 m 3 b 



, nx(m\ — a 2 ) ^(m^ — a 2 ) n 3 (m\ — a 2 ) , 



(303) 



By analogy, Ri is given by 



Ri = 



N e AB 
N e m e B 

,2 _ A2 



N^AB 



-N^m^B 



N T AB 
N T m r B 



y N e (m 2 e - A 2 ) N M (ml - A 2 ) N T (m 2 T - A 2 ) J 



with 



(304) 



C 

A 2 

B 2 
1 

Nf 



m e m u m T 



m e — rrin + m T ' 



— m e m^ + ni^m-r — m T m e — A 2 , and 



A 2 B 2 + m 2 B 2 + (m 2 - A 2 ) 2 , {I = e, (i, r). 



(305) 



(306) 



Using R v and 7?; from (304) and ( p03[ )in (291), we find the following form for the elements of the 
mixing matrix: 



V a i = n a Ni (abAB + {-lf" 2+s '^m a mibBe' a + (m 2 a - a 2 )(m 2 - A 2 )e^ 



(307) 



for a — 1, 2, 3 and I — e, fi, t.[] 

Obviously the boxes from the Fritzsch mixing matrix will be rather complicated functions of 
the masses. In a parameterization based on the mixing matrix elements, we would have to invert 



equations quartic in the Vs given by equation (307) to find the Fritzsch parameters a, b, et cetera 
from the observable probabilities. The boxes, however, are the probabilities (more or less), and the 



2 Our assignment of numbers to the neutrino index and lepton flavors to the charged lepton index is slightly 
different from our normal convention of a = e,/i, r and I = 1,2,3. Because the charged lepton flavor states are not 
equivalent to mass states in the Fritzsch model, we have chosen to use flavor indices to represent the charged lepton 
flavor. To avoid complicated indexing of the neutrino states, we have labeled them by numbers. In retrospect, this 
appears a bit inconsistent, but it seemed like a good idea at the time. 
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Fritzsch parameters may be found by using equation (E54I) . Unfortunately, we must first rotate the 



Fritzsch M v to M using the charged lepton rotation matrix (304) and equation (257). This again 
gives daunting equations, so we will content ourselves with considering only two limiting cases. 

We know that the charged lepton masses obey a strong hierarchy: m T ^> 3> m e . With 

these constraints, Ri becomes 



Ri = 



^ 1 Ai A X A| ^ 
Ai -f A 2 
-A X A 2 A 2 f J 



where 



(308) 



Ai 



and Ao 



(309) 



If neutrinos also obey a strong mass hierarchy, 7713 3> m 2 3> mi, R v will have the same form 
(|308|), with A substituted for A and 



Ai = 



mi /m 2 
and A 2 = 



(310) 



The diligent reader may notice that the matrix in equation (308) does not exactly diagonalize Mi. 
But the mass eigenvalues do appear on the diagonal, and the off-diagonal elements are negligible 
under the mass hierarchy approximation. The mixing matrix for a Fritzsch model with dual mass 
hierarchies takes the form 



V = 



1 + AiAie" 7 
Ai - Aie i<T 



^ A 2 A ie iCT - AiA 2 e 



Ai - Aie iCT AiA 2 e iCT - AiA 2 e' 

A1A1 + e ia + A 2 A 2 e l " -A 2 e J(J + X 2 e^ 
-A 2 e lCT + A 2 e l " 



J 



(311) 



This result disagrees with the CKM matrix given in |45| and [g9[ by the sign of V\%, V21, V23, and 
V32 . One possible explanation of this discrepancy is that the authors of ji^l and |^] have absorbed 
an extra phase e l7r into the fields of the muon and muon-neutrino. Such an absorption will not affect 
the measurable boxes. Indeed, every box has an even number of the terms with the disputed sign, 
so the extra minus signs cancel. 
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The boxes formed from the CKM matrix in ( pll| ) are rather complicated functions of mass ratios. 
For example, 



11 



□ 



22 



33 



= e la (A1A1 - h\e la - \\e w + \ x k x e 2la ) , and 
= e lT > [-\ 1 \ 2 K 1 K 2 e 2ia - (AiAlAi + AiA x A^) e 4(CT+J)) - (A^A 2 + A 2 A 
+ X 2 K\K 2 e l{ - a+2 ^ - AiAaAiAae 24 " - A 2 A 2 A 2 e 3lff 



2\ J(2°+V) 



(312) 
(313) 



(314) 



Under this double-mass hierarchy, the dominant boxes are 22 a 33 , which are on the order 
of a mass ratio (two powers of A or A). The other boxes are on the order of a mass ratio squared. 
11 D 22 appears in the expression ( |135| ) for the probability of a v e «-> transition, multiplied by 
sin 2 $12, where $12 is given by (|3) to be 



$12 



4p 



(315) 



which is on the order of m| £ for small x, so its contribution is lessened by X 2 compared to boxes 
multiplied by $23 or $13. This reduction more than cancels out the effect from the box's dominance. 
22 D 33 , on the other hand, appears in the probability of a «-> v T transition, multiplied by sin 2 $23, 
making the <-* v T transition the most probable transition in the double-hierarchy Fritzsch model. 

Another popular scenario for neutrino mass is the reverse hierarchy: m\ S> m 2 3> Using 
this approximation in the Fritaxch model yields a rotation matrix given by 



^ AfjAj 3 A' 2 



\ 



A', 



An 



with 



A'i — A 2 A! J 



(316) 



A'i = an d A 2 = . I — -. 

V mi V m i 



(317) 



Using a straight hierarchy for the charged leptons and this reverse hierarchy for neutrinos, we arrive 
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at the following form for the mixing matrix: 



V 



' X' 1 A 1 e la - AiA 2 e iv 
X' 2 - Aie" 7 
1 + X' 2 A 1 e ia 



-X[e t,T + A 2 e lr > X\A 2 e lG + e"' 

A 2 Ai + e ia + AiA 2 e i " -A 2 e l,J + X' x e i7 i 



Ai - X' 2 e l 



X' 2 A 2 e ia — A' 1 A' 2 e i ' ) , 



(318) 



Under this approximation, the dominant boxes are 21 CI32 and 12 d 23 , which contribute to 
v t mixing scaled by sin 2 $12 and to v e <-* i/^ mixing scaled by sin 2 $23, respectively. In 
this case, sin 2 4> 12 is large at small distances, so once again the v„ «-> v T mixing predominates for 
neutrinos from local sources. These predictions of the Fritzsch model will be compared against 



experimental data in Chapter VI 



If neutrinos are not pure Dirac particles, then the matrix in ( |29l| ) does not necessarily apply. 
But if a light-neutrino 3x3 matrix M v is obtained by the see-saw mechanism M v = M]jMg Md, as 



in Section 2.6, then the mixing matrix may have the form in (291), as long as Ms, the mass matrix 
for the right-handed neutrino singlet term, is proportional to the unit matrix in the basis where Md 
is diagonal [B9j. 



5.2.2 Using a Democratic Mass Matrix 
Another texture of mass matrix popularized by the disparities in quark masses is the "democratic" 
mass matrix [fl6l: 



(1 1 A 



M = 



1 1 1 
1 1 1 



(319) 



This matrix has eigenvalues 0, 0, and 3. The mass matrix must be normalized to set the non-zero 
eigenvalue equal to the heaviest mass. We know from experiments the mass of the top quark is much 
heavier than those of the up and charm quarks, and the bottom quark mass is much heavier than 
the other down-type quark masses, so the democratic mass matrix is a reasonable approximation for 
the quark sector. The charged leptons also have one dominant mass scale, and it is not unreasonable 
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to assume the neutrino masses might follow the same pattern. So we have 



fx 1 1^ 



M„ = 



m 3 



1 1 1 
1 1 1 



and Mi = 



1 1 1 
1 1 1 



(320) 



Because of the degeneracy of two of the masses in each matrix, the eigenvectors depend on an 
undetermined parameter which we call a. They have the forms 



Ip 1 = Til 



\ 



\ -(o+l) / 



( a + 2 \ 

-(2a +1) 
a- 1 



and V 3 = 3 3 



1 

V 1 / 



m and n 2 are fixed by normalizing the eigenvectors: 



(321) 



m = 2(a 2 + a+l); 
n 2 = 2(3a 2 + 4a + 3). 



(322) 
(323) 



As above, we will choose lowercase a and for the neutrino matrices, and uppercase A and Ni for 
the charged lepton matrices. These eigenvectors may be used to form the R v and Ri that diagonalize 
M v and M\. Constructing the CKM matrix from the R matrices, we find a rather simple result: 



^ niN e (2 + 2aA + a + A) mN^a - 3A - 2) 0^ 



1/ = 



V 



n 2 ^V e (3^-3a-2) 3n 2 N fl (2 + 2aA + a + A) 
1 



(324) 



a and A are completely arbitrary; the eigenvectors corresponding to the lighter particles are degen- 
erate, so one can rotate freely between the two states with no observable effect. We can represent 
this freedom by parameterizing the mixing matrix with rotation angles 6 for the neutrino sector and 
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for the charged-lepton sector. Orthonormality of the eigenvectors is met by setting 



N X A 



cosw 


sm 9 




v/6' 


COS0 


sin 6 


V2 


f 



and 



(325) 
(326) 



The mixing matrix has a quite obvious form with this choice of parameterization: 



1 cos(6 - 9) sin(0 -9) ^ 



V 



- sin(0 - 6) cos(6» - 6) 
1 J 



(327) 



The lighter states are rotated by a combination 9 — 9. If the charged lepton sector is described by 
the same angle as the neutrino sector, no mixing occurs at all. Also, if the lighter-mass degeneracy 
is exact, there is complete freedom to choose 9 — = 0. One expects, however, that "accidental" 
symmetries will be broken by higher-order effects. Parameterizing this symmetry breaking of the 
2x2 light mass matrix as 

/ 

cti e 

one finds 



(328) 



tan 29 = 



2|e| 



Cl — 02 



(329) 



and 



Ami 



(ax + o 2 ) ^(o x - ( r 2 ) 2 +4|e| 2 . (330) 
Whether one chooses a and A or 9 and 0, the mixing matrix as given has only one non-zero box: 



□ 22 = t4t4N(N^(2 + 2aA + a + A) 2 (3a -3A- 2) (3 A - 3a - 2) = -- sin 2 (2(6» - 6)). (331) 



This box appears in the expression (135) for the probability of a v e — > f M transition, multiplied by 



93 



sin 2 $12, where $12 is given by (|^) to be 



$ 12 = ~ L ^x, (332) 



which is very small for a nearly exact democratic mass matrix except in the case of extremely large 
x. Symmetry breaking may also induce small nonzero boxes mixing the third mass with the first 
two. Here $13 and $23 would be large, but the amplitude of oscillation would be small. Oscillations 
may therefore not be measurable in this scenario for terrestrial experiments. Experiments looking 
for oscillations from cosmologically distant objects might, however, be sensitive to the v\ <-> V2 
transition p7j . 

The democratic mass matrix is only an approximation. The masses of the electron and muon are 
not exactly zero, and we reasonably expect non-zero masses for the two lighter neutrinos too. The 
masses of the muon and electron are, however, much smaller than the tau mass, so the democratic 
mass matrix may be a good first approximation, with the smaller masses arising from corrections to 
it. 



5.2.3 Using a Zee Mass Matrix 
As discussed in Section [2.3.2 , the Zee mechanism with diagonal coupling between a Higgs doublet 
and charged leptons gives rise to a neutrino mass matrix of the form (^9|). Such a matrix may be 
parameterized in terms of three parameters a, a, and m pa], p6[, (l2[: 



( 



\ 






a 


cos a 


a 





sin a 


cos a 


sin a 






(333) 



These new parameters are related to the charged lepton masses and the coupling constants of ( |49| ) 

by 



fe 



cos a = m 

m 



fe 



f er f lit 
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sin a 



— vK 

m 



f er f fir 



(334) 



cos a, and 



mo 



fer-ml 
m T \J fer + ffiT ■ 



These relationships ignore m\ and with respect to m 2 T , which is a fairly good approximation 
given the mass hierarchy of the charged leptons. We will also assume a <C 1, which holds unless 



/e M ~10 4 /er 1 12 



Upon diagonalizing the matrix in equation (|333|) , we find the neutrino masses |1J 



-TOocrsin2o; 



1 



to 2 = toq ( 1 — -cr sin 2a 



(335) 



m 3 = m [ 1 + — u sin 2a ] , 



where an overall minus sign on has been absorbed into the definition of the neutrino field through 
the chiral rotation ip — > i"/§ip. The mass differences for this Zee model are 



ml 



ml 



-Toq (1 — 2cr cos a sin a) , 

-Toq (l + 2cr cos a sin a) , and 



(336) 



— 4to i7 cos a sin a. 



TO3 and TO2 are nearly degenerate and much larger than m%. 

Because the lepton mass matrix has been assumed diagonal for this model, the weak-interaction 
mixing matrix V is equivalent to the neutrino mixing matrix U VL , and is found by diagonalizing the 



mass matrix in equation ( 333 ) 



^ \/2 cos a sin a + a cos a sin a + a cos a 



V = 



V2 



sin a cos a 

-V2a 1 



cos a 
-1 



(337) 
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The boxes from this mixing matrix have the forms 



11 D 2 2 = — — sin a cos 2 a (sin a + a cos a) 
11 ^23 = — — sin a cos 2 a (sin a — a cos a) 



u 32 



1 

■ —a cos a (sin a + a cos a) 



11 CI33 = — -cr cos a (sin a — cr cos a) 

12 D 23 = \ cos 2 a (sin 2 a - a 2 cos 2 a) (338) 



4 

4 



12 D 3 3 = - - (sin 2 a - a 2 cos 2 a) 



21 1 

□32 = — cr cos asm a 

21 1 

□ 33 = ——cr cos asm a 

22 D 33 = -icos 2 a 



If fer ~ fun then sin a ^> cr cos a and the boxes have the following relative magnitudes: 

"□22 ~ n n 2 3 12 °23 ~ -sin 2 acos 2 a, 

n n 3 2 - n n 3 3 ~- 21 n 32 ~ 21 n 33 ~ -<r sin or cos or, (339) 

12 D 33 ~ —sin 2 a, and 

22 D 33 ~ -cos 2 a. (340) 
If, on the other hand, f eT 3> fy. T we may neglect the first terms of the sums in the box relationships 



(338). In this scenario we find 



11 a 2 2 ~ — ^23 ~ — crsinacos 3 a, 

"□32 ~ - n D 33 ~ 12 D 33 ~ -cr 2 cos 2 a, (341) 

12 n 23 ~ -a 2 cos 4 a (342) 

21 D 32 ~ — 21 n 33 ~ cr sin a cos a, and 

22 D 33 - -cos 2 a. (343) 
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The predictions of both of these scenarios are compared with experimental results in Chapt 
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CHAPTER VI 



THE CONNECTION TO EXPERIMENT 



The advantages of the box notation become clear when drawing conclusions from experimental 
results. To that purpose, we will examine the currently existing data. No direct evidence of neutrino 
oscillation has yet been confirmed. The LSND experiment at Los Alamos has preliminary data 
which may be explained by oscillations of muon antineutrinos into electron antineutrinos [jl8| . This 
oscillation interpretation has not yet been confirmed by another experiment, but we will consider 
LSND's result as a possible measurement of an oscillation probability. Other accelerator experiments 
put bounds on different oscillation probabilities, and we will use these bounds as constraints in what 
follows. Indirect evidence for neutrino oscillations comes from the solar neutrino deficit §,[| and the 
atmospheric neutrino anomaly Q , and we will use the oscillation interpretation for these phenomena 
also. The probabilities derived from such indefinite results are not meant to be definitive, but they 
will allow us to demonstrate the usefulness of the box notation. Once more conclusive data becomes 
available, the numbers used below may be adjusted, but the procedure will remain the same. 



6.1 The Data 

Solar neutrino experiments measure the flux of neutrinos from the sun and compare it to the 
flux predicted by solar models such as Pinisonneault and Bahcall's [Q. With the exception of 
Kamiokande, the older solar neutrino detectors are sensitive only to electron type neutrinos, so the 
ratio of observed flux to predicted flux gives the survival probability for electron neutrinos, P ^ . 
The Kamiokande experiment measures neutral current events too, so it may detect muon and tau 
neutrinos. The neutral-current cross section in water is 0.17 (about one-sixth) times the charged- 
current cross section for the neutrino energies measured by Kamiokande [ j50| . Thus the probability 
measured by Kamiokande (and now by Super K) is 

P^ ue + 0.17(P e ^ + + P^J, (344) 



()<s 



where P ^ is the probability as measured by charged-current-only detectors. The analysis of solar 
neutrinos used below Jsj]] takes into account this additional contribution to the Kamiokande data. 
The quoted solar deficits contain hidden uncertainties, because the expected neutrino spectrum given 
by solar models is itself somewhat uncertain. 

Atmospheric neutrino experiments measure the flux of muon and electron neutrinos produced in 
the atmosphere. The chain of events producing these neutrinos should create roughly two muon-type 
(i/p or ) neutrinos for every electron-type (y e or V e ) neutrino. Experiments, however, have mea- 
sured nearly equal amounts of electron- type and muon-type neutrinos. These results are expressed 
as a ratio of ratios, R a tm = ^expected/ f measured, where the ratios r are of electron-type to muon-type. 
For three generations, this anomaly could be due to v a <-> v e mixing, <-> v r mixing, or three-way 
mixing v e v r , and these three possibilities are described in detail in the paper which first 

suggested oscillations as a solution to the atmospheric neutrino anomaly, ]5^ ], We will abbreviate 
r expected as simply r, and take its value to be 0.48 |p3[ , the result of sophisticated Monte Carlo 
simulations assuming no oscillations. This number is somewhat dependent on detector efficiencies, 
but the differences between detectors will not affect our conclusions significantly. 

If the atmospheric anomaly is caused solely by v a <-*■ v e mixing, the ratio depends only on P _^ 
and its time-reversal conjugate P ^ : 

1 - P +rP _ 
R ^m = x _ p ^'\ r ^ ■ (345) 

When CP is effectively conserved, such as when sin2$y averages to zero, the two probabilities are 
equal, and the measured ratio is simply 

1-P (1-r) 
Rat m = l _p^ {l _ r -ry (346) 

The transition probability may be found from data by inverting this expression: 



-, (347) 

x 
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where x and y are denned for our use as 

x = 1 - r - Ratm + r^ 1 R atm , and 

y = I -Ratm- (348) 
For the value of r = 0.48 used here, 

a; = 0.52 + l.liW (349) 
If the atmospheric anomaly is caused by <-> v T mixing, the ratio is simply expressed: 

Ratra = 1 - P . (350) 

This possibility is similar to one with «-> v x mixing, or oscillations to both v T and u x , since v T is 
not measured by atmospheric neutrino experiments and therefore acts like a sterile neutrino. The 
transition probabilities in these scenarios are related to the measured ratio by 

1 - P , and (351) 
l-P -P , (352) 

respectively. Probabilities may be found from 



Rat 



P _ - 1 - Ratm - V, (353) 



using the definition of y given in equation (348) 



The final possibility for the atmospheric anomaly involves all neutrino flavors. The measured 
ratio in this case, including a possible sterile neutrino, is 

l-P - P - P +rP 

Ratm = x _ p ^'_ p ^_p^\ r ^ ■ (354) 
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Once some of the probabilities in equation (354) are measured directly by experiment, this four- way 
possibility could become a preferred solution. In what follows, we will confine our consideration of 
this possibility to cases for which P _^ and P ^ are negligible, and for which CP is conserved by 
atmospheric neutrinos. The expression for the measured ratio in this three-way mixing scheme is 

1 - P - P + rP 
Ratm = i _ P " +r -ip T _ p ~ ' ( 355 ) 

which may be rearranged to give 

P +XP -R atm P =V- (356) 

The accelerator experiments considered here, LSND, CHORUS, KARMEN, and CHARM II, 
are, with the exception of KARMEN's v e — ► v x "disappearance" (of an v e to any other flavor 
x = fj,,T,x) measurement, "appearance" experiments, so they measure (or constrain) individual 
transition probabilities. Analysis of these experiments' results is thus simpler than the analysis of 
the "disappearance" solar neutrino and atmospheric neutrino experiments. Many other experiments 
(Bugey in France, Gosgen in Switzerland, NOMAD and BEBC at Cern, CCFR and E531 at Fermilab, 
E776 at Brookhaven, and Krasnoyarsk in Russia) have searched (and/or are searching) for neutrino 
oscillations as well, but none of these have yet found any evidence of oscillations. 

The data from the experiments which we will use is summarized in Table ||. Only the data from 
Kamiokande, SuperK, and 1MB have been listed under atmospheric neutrinos; other experiments 
(Frejus, Soudan2, Baksan, and Nusex) have also measured the fluxes of atmospheric neutrinos, but 
their uncertainties are considerably larger than those of the listed experiments, which makes them 
less useful. 

If the uncertainty in a phase of a neutrino is ~ one-fourth an oscillation length, the oscillatory 
terms in the probability equation average, with sin 2 $y averaging to one-half, and sin2$jj averaging 
to zero. This assumption that the oscillations average thus becomes an assumption about the 
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Table 2: Experimental data pertaining to neutrino oscillations. Sources for the data are listed un- 
der reference When both systematic and statistical uncertainties were given for an experiment, 
they were added in quadrature to obtain the total uncertainty listed below. The measured ratio for 
solar experiments is the measured flux divided by the flux predicted by Bahcall's solar model Jf9| . 
Kamiokande and SuperK do not measure P ^ but the combination given in equation (344). The 
addition to their measured ratios is noted by "+iVC" in the Table. For atmospheric experiments, 
the measured ratio is the ratio of ratios, R a tm- The probabilities for the atmospheric experiments 
are found from R a tm by assuming CP invariance. CP invariance is also assumed for the combined 
probability from the CHARM II experiment. The limits on measured probabilities from the accel- 
erator experiments are reported for the 90 % confidence level. Energies marked with a ~ represent 
peak energies in the energy spectrum for the reaction specified. Accelerator neutrino energies with 
no such demarcation correspond to monoenergetic neutrinos. LSND-1 refers to the decay-at-rest 
results; LSND-2 refers to the newly-presented decay-in-flight results, for which values of L and E 
are not presently available. 



Type 


Experiment 


Length 
(mj 


Energy 

1 1V1C V J 


Measured 


v a 


vp 


P (x) 




Homestake 


1.5 x 10 11 


> 0.814 


0.28 ± 0.04 


Ve 


V e 


0.3 




Kamiok. 


1.5 x 10 11 


> 7.5 


0.50 ± 0.06 


Ve 


Ve+NC 


0.4 


Solar 


SAGE 


1.5 x 10 11 


> 0.233 


0.53 ± 0.10 


Ve 


Ve 


0.5 




Callex 


1.5 x 10 11 


> 0.233 


0.51 ± 0.06 


V e 


V e 


0.5 




SuperK 


1.5 x 10 11 


> 7.5 


C\ A A I C\ f\ A 

0.44 ± 0.04 


V e 


v e + NC 


0.4 












v M 


v T , v x 


0.46 




1MB 


10 4 to 10 7 


^10 3 


0.54 ±0.09 






0.37 












3-way 


x = l.l, 
y = 0.46 












V M 


v T , v x 


0.40 




Kamiok. 


10 4 to 10 7 


^10 3 


0.60 ±0.06 






0.30 


Atmos. 










3-way 


x = 1.2, 
y = 0.40 












v M 


v T , v x 


0.33 




SuperK 


10 4 to 10 7 


^10 3 


0.67 ±0.08 




Ve 


0.27 












3-way 


x = 1.2, 
y = 0.33 




LSND-1 


30 


- 50 






K 


3.l± 1.3 x 10~ 3 




LSND-2 










K 


2.7± 1.4 x 10~ 3 








30 




v» 


Ve 


< 0.026 




KARMEN 


17.6 


- 35 




Ve 


Vr, Vf,, V x 


< 0.197 


Accel. 






- 50 






Ve 


< .0038 




CHORUS 


590 


- 27000 






V T 


< 0.004 














Ve 


< 0.0047 




CHARM II 


650 


- 20000 




K 


V e 


< 0.0026 












both combined 


< 0.0028 



magnitude of the mass-squared differences: 



a * > 77 9 ? > 

A$i, ~ — =>■ mj — m - 



2tt 



3 A(f) 



(357) 



where the symbol A represents the uncertainty of a quantity. In our data analysis below, we will 
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consider only mass-squared differences for which this type of averaging occurs for both solar neutrinos 
and atmospheric neutrinos. 

According to stellar evolution theory, solar neutrinos are produced in the inner part of the sun 
below 0.2i?Q |^0), so the distance traveled varies between 1 a.u. — (0.2)R@ and 1 a.u. + (0.2)i2©. 
(1 a.u. — 1 astronomical unit = 1.5 x 10 8 km, and Rq = solar radius = 7.0 x 10 5 km.) This 
leads to an uncertainty in distance of less than one percent. The energy of solar neutrinos, however, 
ranges from 0.2 MeV to 25 MeV Q. Thus f varies from 3.7 x 10 12 eV~ 2 to 3 x 10 10 eV~ 2 , and 
we find 

TO 2 _ m 2 £ 2tt (2.7 x 10~ 13 ) eV 2 w 1.7 x 10~ 12 eV 2 (358) 



as a rough criterion for the solar neutrino oscillations to average. A more refined calculation would 
include an energy integration over the neutrino spectrum and the detection cross section starting 
from the detector energy threshold, but the above estimate is sufficient for our purposes. 

Making this same assumption for atmospheric neutrinos results in a stronger constraint on mass- 
squared differences. The detected neutrinos separate by energy into three distinct groups: contained, 
stopped, and through-going . Neutrinos in contained events produce charged leptons in the detector 
that are absorbed by, and thereby contained in, the detector. These events are caused by lower- 



energy neutrinos, with energies ranging from 0.2 GeV to 1.2 GeV |53|, ]56|], |52 . Because the lepton 
is formed in the detector, down-going background events from cosmic rays are easily identified. 
Down-going signal events dominate the contained data, due to the closeness to the Earth's surface. 

Stopped events are caused by upward-going muon neutrinos which create a muon outside the 
detector which travels into the detector and is stopped. Electron neutrinos do not contribute to this 
type of event since electrons will not travel far through the material surrounding the detector, and 
down-going events are vetoed because the down-going background is significantly larger than the 
signal. The muon neutrinos producing the stopped muons are more energetic than those causing 
contained events (since the muons from stopped events may travel further than the width of the 
detector), having energies between 1 and 100 GeV f36|] . Still higher-energy up-going muon neutrinos 
may produce muons outside the detector which travel through the detector and out the other side. 
These through-going events result from neutrinos with energies between 1.1 GeV and 10 4 GeV [p6[. 
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In our data analysis, we follow the lead of reference [ pi] and consider only the contained events, 
which are identified roughly ten times as often as events of the other two types J5(| . These events 
have an energy uncertainty of about a GeV. Because down-going events dominate, the majority 
of contained events' neutrinos originate in a cone of about 90° about the vertical J5l| . Thus the 
distance traveled by a contained-event neutrino varies from 30 km, the depth of the atmosphere, 
to 6.5 x 10 3 km, the radius of the Earth. Contained events therefore have an uncertainty in % of 
^OcV 1 ~~ l 2°GcV ~ 1-7 X 106 eV ~ 2, This unccrtaint y causes the oscillations to smear out if 

m, 2 - > (p g x 1Q -7-J e y2 _ 3 y x 10 -6 e y2 

When this average is made, the probabilities become simple sums of boxes: 



(K a ^) = -2 ( %2 + %3 + %a) ■ (360) 



As discussed following equation (157) in Chapter IV, this row sum equals the sum \Vax\ 2 \Vp 12 



We cannot average oscillations for accelerator experiments. Upon observation of Table we 
see that the oscillatory behavior must occur if the LSND results are to be consistent with the limit 
placed by CHARM II. 



6.2 Interpreting the Data 

The atmospheric and solar data may be analyzed independent of neutrino masses, provided the 
mass-squared differences are large enough to yield the averaging described above. The probabilities 
measured at accelerators, however, are functions of these mass-squared differences. We will follow 
the authors of reference [j5lj and choose the mass-squared difference Am| 3 which puts the LSND 
experiment at an oscillation maximum, and the difference &m\ 2 to be much smaller: 

m l ~ m l ~ m 3 ~ m \ ~ 2 ev2 > and m \ ~ m i ~ 1(r2 ey2 - ( 361 ) 
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The smaller difference Am] 2 is small enough to be negligible in accelerator experiments but large 



enough to be above the limits set in equations (358) and ( 359 ). With this choice, oscillations average 
and CP is effectively conserved for both solar and atmospheric neutrinos. 

The oscillatory factors for accelerator experiments are listed in Tabic || for our choice of mass- 
squared differences. We have used the typical E quoted in Table |^ to obtain the first column of 
numbers. For each of the ^ given in the Table, sin 2 $12 and sin2$i2 are much smaller than the 
sin functions involving the other mass-squared differences, which are roughly equal. Under these 



conditions, the sum of sin functions in equation (171) is nearly zero. (The biggest value for this sum 
comes from the LSND experiment and is ±0.015J\) CP is therefore roughly conserved, and P ^ 
measured by an accelerator experiment will be approximately 



p acce 



1 -2 ( % 2 sin 2 $ 12 + % 3 sin 2 $ 23 + ^3 sin 2 $ 13 ) 

i -2 (~ 1(T 4 % 2 + %3 + %a) sin 2 $ 23 . 



(362) 



Unless the real part of al Op 2 is much larger than the other real boxes, its term will not contribute 
significantly to the transition probability. 

Table 3: Values of oscillatory terms for accelerator experiments. The numbers were calculated 



using the mass-squared differences given in equation (361). The typical values of energy indicated 
in Table || were used to calculate the second column. A more complete treatment would integrate 
over the entire energy range, but the approximation we have chosen is sufficient for our purposes. 



Experiment 


M^ 2 ) 


sin 2 $12 


sin 2 $23 


sin 2 $13 


sin2$i2 


sin2$23 


sin2<I>i3 


LSND 


0.76 


5.8 x 10- b 


1.0 


1.0 


0.015 


0.10 


0.10 


KARMEN 


0.64 


4.1 x 10" 6 


0.92 


0.92 


0.013 


0.55 


0.55 


CHORUS 


0.028 


7.8 x 10~ s 


3.1 x 10~ 3 


3.1 x 10~ 3 


5.6 x 10~ 4 


0.11 


0.11 


CHARM II 


0.041 


3.7 x 10~ 2 


6.7 x 10~ 3 


6.7 x 10~ 3 


8.2 x 10~ 4 


0.16 


0.16 



Acker and Pakvasa | pl[ have analyzed the solar neutrino data, including the possible contributions 
from muon and tau neutrino neutral current interactions to the Kamiokande data. They find that all 
the solar neutrino data is consistent with (P _^ ) between 0.4 and 0.55 for reasonable values of the 
initial 8 B neutrino flux from the sun. These values were attained by assuming the oscillatory part 
of the probability sin 2 (3>y) averages to one-half for all three mass-squared differences, so (P ^ ) = 
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1 — (P ^ ) — (P )• The range of probabilities allowed by the solar neutrino deficit is therefore 



0.45 ~ (P ) + (P _ ) ~ 0.60. (363) 



Using Kamiokande's results (which are comparable to IMB's and SuperK's numbers but with 
smaller error bars) with a range of 1.5a, the atmospheric neutrinos in this CP-conserving scenario 
should obey 

< 1 -(P^)-(P „)+r(P„^) < 
0.51 ~ , ° 7 ,7 ; ; - 0.69. (364) 

Up to this point we have not used our choices for neutrino masses, except to argue that oscillations 
average for atmospheric and solar neutrinos. 

The LSND experiment measures P _^ % , which equals P _^ by CPT conservation. Since only the 
decay-at-rest (LSND-1) result was available when this analysis was carried out, we restrict ourselves 
to consideration of that data. Using that measured probability and our choices of mass-squared 
differences, we find the constraint (again using 1.5 standard deviations as our limit) 

1.2 x 10~ 3 ~ -2 (10~ 4 % 2 + %3 + ^23) sin 2 $ 23 ~ 5.1 x 10~ 3 (365) 



from equation (362). Since 10 4 ^22 < 10 3 for reasonable values of ^22 and sin 2 $23 = 1 from 



Table [|, this constraint requires 

2.5 x 10" 4 < - ( % 3 + % 3 ) < 2.8 x 10~ 3 . (366) 

KARMEN's 90% confidence level places a tighter upper bound than the 1.5cx bound from LSND: 

- (^23 + ^23) < 1-0 x 10- 3 (367) 

By a similar argument, the constraints from CHORUS and KARMEN may each be turned into 
a constraint on the sum of two boxes for the other oscillation possibilities. From CHORUS, we find 
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a limit on v^-v.,- mixing: 

- (£33+ £33) < 0.32. (368) 
KARMEN provides a limit on v e -v r mixing: 

- (IS33+ ^33) < 5.4 x 1(T 2 . (369) 



We may combine the constraints on v e -v^ mixing with equation (217) to eliminate S23, and solve 



for If we call the sum — ( S23 + ^23) as measured by LSND or other accelerator experiments 



%3 = -,±V,<-4^s-43< pn) 

If TB.23 is to be real, which it is by definition, 

AJ 2 < s 2 -4tt 2 2S. (372) 

Assuming — M.22 has its maximum value of half (P (x) ) mm , or 0.3, and s equals 10~ 3 , the maximum 
value allowed by KARMEN, we find the following minimal constraint on J 2 : 

J 2 < 3 x 10" 4 . (373) 

Even though CP is essentially undetectable due to our choice of masses, the imaginary part of 
the boxes, J7, may still play a role in determining the oscillation probabilities. Using the above 



maximum values of — ^22 and s along with this maximum value of J 2 in equations (370) and (371) 
yields 

%3 = %3 = ~ = -5 x 10- 4 . (374) 



107 



If, on the other hand, J 2 = 0, 



(&23, $23) = (1.7 x 1(T 2 ,1.8 x 1CT 2 ) or (1.8 x 10~ 2 ,1.7x 1CT 2 ) (375) 

A set of measurements of probabilities at different distances which pinpointed the values of individual 
boxes would thus probe the value of J' 2 , even if CP were hidden from observation in the oscillations 
being observed! The above ambiguity between the values of ^23 and ¥£23 arises because sin 2 $23 ~ 



sin 2 $13. We have chosen one mass scale to dominate the other, a situation examined in Section 5.1.4 



Consideration of equation (276) shows that only the sums "lips + °!kpa are measurable. Unitarity 
constraints of the form (217) for a known J 2 and °kp2 allow us to determine the choices for 



and °Kf33 but not to say which has which value. 



6.3 Finding Solutions 



As promised in Section 1.2.2, we numerically find boxes which match the data. Different values 
of the five parameters ^22, ^23, ^32, ^33, and J 2 are stepped through, the other five Rs found by 



equations (217) to ( |22l| ), and the resulting oscillation probabilities compared to the solar neutrino, 
atmospheric neutrino, and LSND results discussed above. The boxes must also produce a mixing 
matrix which is unitary. 

We find many solution sets which meet all of these conditions. For J 2 = 0, the two best fits are 
found with 



B l = 



which is symmetric, and 



Bo = 



-18 


-1.0 


0.95 


-1.0 


-0.67 


-2.0 


0.95 


-2.0 


-1.8 


-10 


0.91 


-1.0 


-13 


-1.0 


-0.93 


5.6 


-10 


-13 



x 10" 



(376) 



x 10" 



(377) 



/ 



The first solution gives a solar ratio (P ) of 0.42, an atmospheric ratio R a tm of 0.55, and an 
LSND probability of 2.1 x 10~ 3 . The second gives a solar ratio of 0.55, an atmospheric ratio of 0.69, 
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and an LSND probability of 3.6 x 10 



We may obtain the magnitudes \V a i\ from the boxes from equation (|143| ). Using the notation 
\V\ for the matrix of these magnitudes, the two solutions yield 



/ 



and 



\ 



0.88 


0.47 


0.20 


0.50 


0.88 


0.12 


0.12 


0.20 


0.99 


0.50 


0.41 


0.74 


0.65 


0.73 


0.04 


0.53 


0.51 


0.65 



\ 



(378) 



(379) 



Inspection reveals appropriate places to insert minus signs to make these matrices unitary to within 
the tolerance required by our numerical technique. Such an assignment, however, is not unique due 



to the phase ambiguities discussed in Section 4.1 



Comparing the solutions ( |376] ) and (377), we see that by the first solution, v e -v^ mixing comprises 
most of the solar neutrino deficit, while the second solution uses v e -v T mixing to account for most 
of the solar deficit. 

Acker and Pakvasa in [[sif expressed their solution in terms of the mixing matrix: 



V A p = 



I 0.70 0.70 0.14 ^ 



-0.71 0.69 0.12 
-0.010 -0.19 0.98 



(380) 



We find the corresponding box matrix to be 



B A i 



i -24 0.84 -0.87 ^ 



-0.92 -1.6 0.87 
-0.92 -1.8 -0.96 



x 10" 



(381) 



This solution was not selected by our technique because the resulting atmospheric ratio, R a tm — 0.50, 
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is not within the 1.5a range of Kamiokande's data. Acker and Pakvasa's solution yields (P^ _^ ) = 
0.48 for the solar neutrino deficit, and an LSND probability of 1.2 x 10~ 3 . 

The above solutions all contain no CP violation. We found solution sets for CP-violating cases 
too, with J 2 going as high as 2.7 x 10 -4 . The only solution for that maximum J 2 is 



^ -22 -0.12 5 x 10~ 4 ^ 



B 3 = 



-1.0 -1.0 
0.84 -2.3 



0.85 
-3.2 



x 10" 



(382) 



yielding 



^ 0.73 0.63 0.25 ^ 



W\ a = 



0.65 0.74 0.14 
0.21 0.19 0.95 



(383) 



This CP-non-conserving solution gives a solar ratio of 0.53, an atmospheric ratio of 0.51, and an 
LSND probability of 4.9 x 10~ 3 . None of the solutions which conserved CP found this combination 
of high solar ratio and low atmospheric ratio. The CP-conserving solution set which most closely 
reproduces these CP-violating probabilities is 



B 4 



( 



V 



-10 -1.1 1.0 
5.0 -16 7.3 
-10 1.0 -1.1 



\ 



x 10" 



(384) 



which produces a solar ratio of 0.40, an atmospheric ratio of 0.52, and an LSND probability of 
4.5 x 10~ 3 . One can see that CP violation may have a measurable effect, even if the mass-squared 
differences conspire to hide direct evidence of it. 
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6.4 How Do Specific Models Stack Up? 

None of the models presented in Chapter |v| are compatible with the solutions we have presented 
here. Wolfenstein's parameterization requires 

#23 « - #23 « #33, (38 5) 



as seen by examination of equation (269). The first approximate equality is met by all three solutions, 



(376), (p77[) , and (381), but the second is not met by any of them. The small-angle approximation 
is merely a special case of the Wolfenstein parameterization, so it is not consistent with the current 
data either. 



Equations ( |339D and (341) of Chapter |y| present two possibilities for Zee-model- induced oscilla- 



tions. Both possibilities predict 11 D 2 2 ~ 11 ^23 > a condition which is not met by any of the above 
solutions. This equality is very hard to reconcile with both the LSND results and the solar neutrino 
measurements, but it could become relevant if matter effects enhance the effect of an otherwise small 



11 D22 • And as discussed in Section 2.3.2 , we have only examined the Zee mechanism for two limits 
of coupling constants; our work does not address the applicability of other coupling-constant choices. 
The Fritzsch mass matrix model of Section |5.2.1 predicted that for a mass hierarchy, the diagonal 



elements n D22 and 22 D33 should be much larger than the other elements. These two boxes are of 



the same order in solution 2 (377), but they are not significantly larger than the off-diagonal CI32 
and 12 CI33 . The reverse mass hierarchy considered in that Section predicts the same two dominant 
boxes, once we swap the 3 and 1 mass indices to account for the new labeling of mass states, so 
it is inconsistent with our solutions too. Again, our work does not comment on the more general 
Fritzsch model, but only on its limits when one mass scale dominates. 

Before concluding this Chapter we must remind the Gentle Reader that the solutions derived 
above assume an energy-independence of the solar neutrino data and the atmospheric neutrino data. 
Choosing this interpretation allows the fitting of all current data with only three neutrino flavors, 
but we do not mean to present this choice as the only one possible. Many other excellent analyses 
of the data, such as those by Fogli and Lisi |57j] and those by Cardall and Fuller |58j, take different 
approaches which are equally valid. 
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The neutrino-oscillation interpretations of all three data sets considered here (solar, atmospheric, 
and LSND) are still unconfirmed and controversial. If an energy dependence of solar neutrino 
oscillations is unequivocally demonstrated by future experiments, or if any of the probabilities used 
here change dramatically (or keep the same mean value with smaller error bars), the solutions 
we have found will no longer be valid, and another neutrino generation and/or matter effects on 
oscillations may be implicated. At the present, however, our solution sets fit all of the existing data 
and provide a demonstration of the use of the boxes. 
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CHAPTER VII 



OUTLOOK 



Neutrino oscillations address many current physics questions. Observation of them would present 
the first contradiction of the Standard Model and provide a clue as to the direction that extensions 
of the Model should take. Understanding what role, if any, oscillations play in the solar neutrino 
deficit will help refine theories of stellar evolution. And indications of the magnitude of neutrino 
masses from oscillation experiments would indicate how much neutrinos contribute to dark matter. 
With approximately 115 low-energy "relic" neutrinos of each flavor occupying every cubic centimeter 
of space [^9) , a single neutrino mass on the order of 7 eV (or all neutrino masses summing to 7 eV) 
provides enough mass to close the universe |l2| . 

Once neutrino oscillations are detected and the boxes and neutrino masses determined, neutrinos 
may be used to probe the universe. Neutrino oscillations could help determine the density of different 
regions of the earth |50) and bring information about cosmologically distant objects [[ItJ . Neutrinos 
can act as a long-wavelength telescope, peering further out in space (and therefore further back 
in time) than conventional telescopes. The distance for which neutrino telescopes will be sensitive 
depends on the values of mass-squared differences. This dependence is illustrated in Figure [uj a plot 
popularized by J. G. Learned. The neutrino telescope AMANDA is already taking data, hoping to 
use gamma-ray bursters and supernovae to measure neutrino masses |6l| . Other neutrino telescopes 
are under construction and will soon join the hunt. 

Neutrino physics has entered a golden age of research. New experiments all over the globe promise 
an unequaled amount of data, from the sun, the atmosphere, accelerators, supernovae, and other 
cosmic sources. SuperK in Japan amassed as much data in its first 102 days of running as all of the 
previous solar neutrino experiments did over their combined lifetimes [ |62| . Analyzing such refined 
data requires a consistent, model-independent approach. In this dissertation, we have developed 
such an approach and then illustrated its use by comparison with the low-statistics data currently 
available. 
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logio^). 




lo §lo(ifiv) 

Figure 11: Lengths and energies probed by neutrino oscillations for different mass-squared differ- 
ences, obtained by setting = ^ so that the Earth sits on a maximum of an oscillation length. 
Lengths of interest are labeled on the vertical axis, with representing the radius of Earth and 
Ho representing the present value of the Hubble constant. Shaded regions represent energies and 
distances typical for the indicated cosmological neutrino source. Included are solar neutrinos, both 
up- and down-going atmospheric neutrinos, gamma-ray bursters both within the galactic halo and 
outside our Galaxy, and active galactic nuclei. Double-bang AGNs refer to those AGNs producing 
tau neutrinos energetic enough to produce a "double-bang" signature in detectors. 
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APPENDIX A 



NOTATION KEY 



Greek Lorentz indices /i, v, er, et cetera may take the values to 3. Latin Lorentz indices a, b, et 
cetera refer to the spatial components 1 to 3. 

Greek particle indices a, 0, 7, et cetera represent states of definite flavor. Latin particle indices 
k, et cetera represent states of definite mass. 

We will explicitly write out all flavor sums, since repeated flavor indices are not necessarily 
summed over. 

Repeated Lorentz indices fj,, v, and a are always summed over. 

represents a single sum over the index i, omitting the term for which i = j. Sums over 
both indices will be represented by J^j X^j- 

In this paper, we use the "west-coast metric" in which 500 = L 9aa = — 1, and the other elements 
are zero. 



e = V Ana in this paper, where a here is the fine-structure constant. In standard particle physics 
units (c = h = 1), our e is the electric charge. In other units it is proportional to the electric charge. 



7,7,7,7, and 75 are Dirac gamma matrices, discussed in detail in Appendix BT 



dn = jr~ is the Lorentz-covariant partial derivative. 

" ox ^ 

9w is the Weinberg angle, defined by equation (g). 
Pl = ~o 5 is t ne left-handed projection operator, and 

Pr — 1 ~ l 2 75 is the right-handed projection operator. They are discussed in detail in Appen- 



dix B.3.1 



p = p is the 3- vector portion of the four-vector p. 
per = p • er is the dot product of the vectors p and er. 

ip is used to denote a generic spinor, or the spinor containing all neutrino states, described below. 

t/j* is a complex-conjugated spinor. 

ipl = ip* T is a Hermitian-conjugated spinor. 

The symbols * and t have the same effects on other objects as they do on spinors. 
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if) = is a Dirac-conjugated spinor. 

— T 

%)= Cip is a charge-conjugated spinor. 



C = i7 2 7° is the charge conjugation operator, discussed in detail in Appendix B.4. 
v a L is the left-handed active neutrino of flavor a — e, fj,, r. 
P aR is its CP-conjugate antineutrino of flavor a; it is right-handed. 
N aR is the right-handed sterile partner to v aL . It has not yet been observed. 
N aL is the left-handed sterile partner to P aR - It is also the CP-conjugate of N aR and has not 
yet been observed. 

oll is the left-handed charged lepton of flavor a. 
ur is its right-handed partner. 
(x R is the CP-conjugate of a^- 
(x L is the CP-conjugate of a R . 
| v aL \ 

ijJaL = \ is left-handed lepton doublet of flavor a. 

V f " / 

vl is an nx-dimensional vector containing all ul active left-handed neutrinos. V R similarly 
contains all til of the active right-handed antineutrinos. 

N R is an n/{-dimensional vector containing all ur sterile right-handed neutrinos. N L similarly 
contains all n R sterile left-handed antineutrinos. 

ip, in addition to being a generic spinor, represents the n-dimensional vector containing vl, N r , 
and any additional sterile fields not related to active fields. 

Il,r is the 3-dimensional vector containing all three flavors of charged lepton spinors a^ iR . 

4>d is the Higgs doublet which couples to down- type quarks and charged leptons. 

(j>E> is the Higgs doublet which couples to up-type quarks and neutrinos. 
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APPENDIX B 



THE DIRAC EQUATION, WAVE FUNCTIONS, AND TRANSFORMATIONS 

In this Appendix we review many aspects of field theory as they relate to neutrinos. As in the 
Introduction, the presentation here is synthesized from many sources. For a more in-depth treatment 
of field theory and weak-interaction physics, the reader is encouraged to consult any or all of the 
references E|, & @, El, @, JH . 



B.l Gamma Matrices 

The Dirac gamma matrices 7^ are defined by their appearance in the Dirac equation: 



[ry^ - to] i>{x) = 0. 



(386) 



The spinor ip(x) is a four-component spinor, and the gamma matrices are 4x4. For Dirac's equation 
to give the Klein-Gordon equation when its operation is performed twice, the Dirac matrices must 
obey the anticommutation relations pi[ 



(387) 



where is the metric, defined in Appendix |a| Dirac noted that the matrices 

1 0) l 

and 7 a = 

t / 



o a 
\ -<r a 



(388) 



satisfy the anticommutation requirement. Each element of the matrices in equation (388) represents 



a 2 x 2 submatrix; <r a are the Pauli matrices, 





(0 






f -i\ 




















V i J 



, and (T3 



1 

-1 



(389) 
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and is the 2x2 unit matrix. The product of the four gamma matrices, 75 = i7°7 1 7 2 7 3 , is useful 
also. In the above standard representation it has the form 



75 = 




(390) 



The representation of the gamma matrices chosen by Dirac is not unique. If a set of 7^ satisfy 
equation (|387|) , then so will any transformation of that set by an invertible matrix Aj^A^ 1 . Two 
other popular representations include the Majorana representation |l2j, 







in which Majorana particles have pure real spinors, and the chiral representation 



(391) 



7° = 







% -a a \ 






K o- a J 



, and 75 = 




(392) 



B.2 Neutrino Wave Functions 



Neutrino solutions to the Dirac equation (386) are described by field operators of the form JlS 



ip{x) 



1 4^ £ [e ipx v s (j>)bt(p)+e-^u s (p) as (p)]. (393) 



u s (p) and v s (p) are single-particle spinors, andpx is shorthand for the Lorentz invariant p^x* 1 . Using 
the standard representation for the gamma matrices, the spinors have the forms 



u s {p) 



E-\-m 
2E 



t 



\ E+m ^ s 



and v s (p) 



E+m 
2E 



/ rr • p 

E+m X s 
\ Xs 



(394) 
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with 




(395) 



b£(p) and a s (p) are the appropriate creation operator for antiparticle (negative momentum) fields 
and the destruction operator of particle(positive momentum) fields, cr are the Pauli matrices, given 



in equation (38E). 

For simplicity, we will take the momentum along the +z direction,^] and the mass to be much 

CT ' D 

smaller than the energy of the neutrino. The quantity E+7 ^ appearing in the spinor definitions then 
becomes 

cr -p \p\ 



E + m E + m 



cr 3 w (1 - e)<7 3 , 



(396) 



where e = Using the standard representation for the gamma matrices, the neutrino wave function 
now has the form 



/ 



^(i,x) 



d 3 p 
(2^1 



E + m 
2E 



(1 - e^b^p) + e- lpx a. + (p) 
-(1 - e)e' px bl(p) + e- v P x a-(p) 
e^ x b^(p) + (1 - e)e- vpx a + (p) 
e l P x bl(p) - (1 - e)e- ipx a-(p) 



\ 



(397) 



where the subscripts + and — on the creation and annihilation operators represent the spin states 
+ | and — |, respectively. 

If we reflect the coordinate system, the momentum is in the — z direction and the quantity in 



1 This choice reflects the convention to quantize spin along the z-axis. When discussing kinematics, we follow the 
convention to choose p along the x-axis, which may at first seem to contradict our first assumption. To be consistent 
throughout, we could either choose p always along the z-axis so p • x = pz in our kinematic discussion, or we could 
choose to quantize spin along the x-axis and redefine the Pauli matrices so cr\ is diagonal with ±1 eigenvalues. Either 
consistent treatment would yield the same results as ours, so we will stick with our convention-motivated choices, 
inconsistent though they seem. 
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equation (B9q) changes sign. The wave function for the reflected system is 



-x) = 



(27T) 



cPp I E + m 



2E 



1 -(1 - e)e ipx b\.(-p) +e- l P x a + {-p) 
(1 - e)e lpa: b t _ (-p) + e- ipx a- (-p) 
e^b^(-p) - (1 - e)e-* :c a + (-p) 
i e ^b]l(-p) + (l-e)e- i P a: a_(-p) / 



(398) 



The exponentials do not change under the transformation because x'p' =tE— (— x)(— p) = xp. The 
differential d 3 p changes sign, but the limits of the integration change sign too, and those two effects 
cancel each other. 



B.3 Helicity, Chirality, Parity, and Parity Violation 



B.3.1 A Review of Handedness 



The wave function in equation (397) contains all spin states. We may, however separate the 
states with right-handed, or positive, spin from those with left-handed, or negative, spin through 
the helicity projectors P' R and P' L . The handedness of states is determined by considering the 
component of spin along the momentum axis. When one puts the right thumb along the direction 
of momentum, if the fingers of the right hand curl in the direction of the spin, the state has positive 
helicity. If the spin is instead in the direction of the fingers of the left hand, the state has negative 
helicity. 

The helicity projection operators are given by ||i"o| 



with 

/ 

S = 



(T 
CT 



(400) 
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in the standard representation, and they are functions of the momentum and energy |1C|]: 



,.; II= ^ ±75 ^il5£]. 



where (3 is the particles velocity in terms of c. These operators are not manifestly Lorentz-invariant, 
but in the relativistic limit they become 

^(±75). (402) 

The relativistic limits of the helicity projection operators are Lorentz-invariant and are called the 
chirality projection operators, Pr l : 



^. = ^(+75), and P L = ^(-75). (403) 



The chirality operators have the properties 

PrPr = Pr, P L Pl=Pl, and P L P R = P R P L = 0, (404) 

which are consistent with their roles of projecting out particular states. The helicity is only the 
same as the chirality if the particle mass is zero, which means the e in the wave function equations 
becomes zero. We will restrict ourselves to that case here for simplicity. Due to the smallness of 
the neutrino masses, we also use the words helicity and chirality interchangeably in the text. In the 
standard basis, the chirality projection operators have the forms 




(405) 
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Applying these projections to the wave function of equation (p97j) and taking e to zero, we find 



d 3 p IE 



(2tt)S V 2E 



d 6 p E + m 



(2tt)5 V 2E 



' e ipx bl(p) + e-'P x a. + (p) 


e ip^b^(p) + e~ ipx a + (p) 




-e'^b^pj+e-^a^p) 


I e'Fb[(p)- e -v a -(p) / 



and 



(406) 



(407) 



The only states surviving the right-handed projection are those with spin pointing in the same 
direction as the momentum, and those which survive the left-handed projection have spin pointing 
in the opposite direction. If we were to apply the projection operators to the reflected wave function 
of equation ( |398| ), the right-handed projection would again choose the spin states aligned with the 
momentum, but it those would have negative spin since the momentum is negative: 



P fl V(*, -x) = 



d 3 p E 



(27T)i 



2E 





e ip*bL(-p) + e- l P x a-(-p) 


i e ipx b t_(_ p ) + e -^ a _(-p) J 



(408) 



Operating on the reflected wave function with Pl pulls out the positive spin states in a similar 



manner. One must be especially careful to change the argument of the original wave function (397) 



and form the new wave function (398) before using the projection operators. The expression in 



equation ( |408| ) will never fall out of any operation on the wave function in equation (406), since the 



projection operation has already thrown out the negative spin states in equation (406) 
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B.3.2 Parity Transformations 
The operation of space reflection is called a parity transformation. Two parity operations must 
return a system to its original state, so the reflected annihilation and creation operators a^(— p) and 
bgT(— p) may differ from the original fields only by a phase: 

«4(-p)=a.(-p)=f; B a B (p), and b^(-p) = bj(-p) = ffcbj(p). (409) 



The parity-transformed wave function in the new coordinates x' — (t, — x) and p' = (E, — p) is then 



f{x') = 



d 3 p jE 



(2tt)* 



2E 



d 3 p [E 



2E 



' -(1 - e)e ipx b'l(-p) + e- ipx a' + (-p) ^ 
(1 - eyP x b'l(-p) + e~ ipx a'_(-p) 
e ipx b2(-p) - (1 - e)e-' l P x a' + (-p) 
y e ipx b'_!(-p) + (l-e)e- i J ,x a / _(-p) ) 

-(1 - e)e ipx r ]b b\(p) + e-^a+fc) 
(1 - e)e^ % bL(p) + e- l f^ a a_(p) 
e**T/ 6 b^(p) - (1 - e)e- i ^7 ?a a + (p) 
y e^? 7fc b t _(p) + (l-e)e-^7 ?a a_(p) y 



(410) 



As mentioned above, e % 



The transformed wave function should be related to the original wave function of equation (397) 



by simple matrix operations, so each element of the new wave function should be proportional to the 



corresponding element of the old wave function. Comparison of equations (410) and (397) reveals 
that j] a must equal —rib- Choosing rj a — 1, we have 



^ (l-e)e l P x bl(p) + e- l P x a + (p) ^ 



f(x') = 



(27T)! 



d 3 p IE + ni 



2E 



_(1 _ e )e ipx bl(p) + e- l P x &-(p) 
-e ipx b\{p) - (1 - e)e~ ipx a + {p) 
-e ijKC bL(p) + (1 - e)e- l P x &-{p) , 



(411) 
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So a parity transformation of a wave function multiplies the wave function by 7 , and changes the 
sign of the three-vector position and momentum which also changes the sign on the antiparticle 
creation operators. 

Operating with parity on a state of definite chirality is equivalent to flipping the chirality of the 
state and using the reflected coordinates. For example, the parity conjugate ^j l (x') of a left-handed 
state iph{x) = Pl^{x) is 



b L (x') = j Mx)=J Pl*P{x) 

( 



d 3 p IE 



(2tt)S V 2E 



-e* I bL(p) + e"* ,M a_(p) 


-e l P x bl(p) + e-^a-ip) , 



(412) 



where we have used the result of equation (407). This parity conjugate wave function is identical to 



the one found in equation (408) 



d 3 p IE 



(2tt)S V 2E 



d 3 p IE 



(2tt)5 V 2E 



e ipx hl(-p) + e-* x a_(-p) 


e^b^C-p) + e-^a_(-p) , 


e ipx {~bl(p)) + e-' l P x a-(p) 


y e ipx {-bl(p)) + e-''P x a-{p) J 



(413) 



Thus we may represent the parity conjugate field fy L (x') as tpfi(x'). 

This chirality flip due to parity may also be seen by considering the properties of the gamma 
matrices. 7 commutes with! but anticommutes with 75, so multiplying a projected state by 7 
changes its chirality: 

7 °F L ^ = P fl7 V, (414) 
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which is a right-handed state. Because of the presence of the matrix 7 in the Dirac conjugation 
operation, taking the Dirac conjugate of a particle swaps its chirality: 

^ = ^pLh1° = V> fl ^V = ^7 0l 4p = 4Pr,l. (415) 

Ergo a left-handed outgoing particle is equivalent to a right-handed projection on an incoming 
particle. 



B.3.3 Parity Violation in Weak Interactions 
The weak charged current involves only left-handed particles, so it is explicitly parity-violating. 
Consider a charged-current two-vertex lepton interaction represented by an effective Lagrangian, 
Jcc^CCai as discussed in Section [L~2] . As discussed in Section 3.3, a typical charged-current inter- 
action can be Fierzed to produce an interaction of the form 

£?/C/ = -^pf&c&c* = [^7"(1 " 75>e] [e 7 ,(l - 75)e] . © 

This coupling explicitly violates parity through the presence of the left-handed projection operator 
Pl. If we were to not a priori assume the Standard Model form of the Lagrangian, we would need 



a more general interaction 1 1 



C°f f =-^= [l^O* (d - C; 75 ) up] [/?0,a] , (416) 



where Oi can be any of the complete set of couplings: l(Scalar), 7 M (Vector), a^ v (Tensor), 7^75(Axial 
vector), or 75 (Pseudoscalar) . The expression — ^'75 appears only in the first factor in equation 
( |416| ). Such an expression could be placed in the second factor as well, but it would not introduce any 
new freedom in the equation. The coefficients of the S, T, and P couplings must be zero at tree-level 
to preserve gauge-invariance. Deviations in the A and V couplings from common strength would 
be evidence of some right-handed current. Experiments measure observables, such as interaction 
rates and energy distributions, not effective Lagrangians. Michel and Bouchiat K7[, |ll| showed 
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that the observable results of an interaction of the form (416) depend only on four so-called Michel 



parameters: p, 77, £, and 5. These parameters arc functions of the coefficients Ci and C[, and they 
are predicted to be equal to |, 0, 1, and §, respectively, for the Standard Model's completely parity- 
violating weak current. Experimental measurements yield values consistent with the Standard Model 
predictions Jl0[: 



p = 0.7517 ±0.0026, 

T) = -0.12 ±0.21, (417) 
£ = 0.972 ±0.013, and 
S = 0.7551 ±0.0085. 

Such a marked agreement makes the introduction of right-handed currents in extensions of the 
Standard Model quite difficult. 



B.4 Charge Conjugation 

In this paper, we will distinguish between a charge-conjugate particle and an antiparticle. 
Dirac originally postulated the existence of antiparticles as the negative-energy solutions to the 
Dirac equation having the same mass and opposite charge as their associated particles. The original 
Dirac equation was charge-conjugation-invariant, so the term antiparticle became associated with 
the charge-conjugate statej^J 

1>l,r = C^E T , (418) 

where C is the charge-conjugation operator, discussed below. With the discovery of charged current 
interactions in the late 1950s, however, came the discovery that charge-conjugation and parity are 
not good symmetries independently. The combination of charge-conjugation with parity, CP, seemed 
to be a good symmetry, so the new negative-energy solution to the Dirac equation was thought to 

2 Our notation here differs from that of reference [ [l4[ since they apparently use ip L to represent the charge-conjugate 
field of the right-handed field, while we use it to represent the conjugate of the left-handed field. 
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have the form 

(V>L,fl(x')7 = CfVl^f = C^ RtL (x'f = ty R , L (x') (419) 

until the early 1960s. With the study of kaon systems came the discovery that even CP is not con- 
served exactly. We now believe that CPT is a good symmetry, where the T stands for time reversal. 
The physical antiparticles therefore correspond to the CPT-conjugate states ^(—x). Because the 
violation of CP is so small, and could possibly be zero in the lepton sector, we will loosely use the 
term antineutrino to refer to the CP-conjugate 

V R = «f = CW T . (420) 



Having clarified the terminology, let us return to the definition of the charge- conjugate field in 



equation (418). C has the following properties BA 



C 1 = C* T = C" 



-C, CC t =€, and O^C" 1 



-7/.- 



(421) 



Using Dirac's standard representation of the gamma matrices, C has the form 



C = i 7 2 7°. (422) 

This identification of the matrix C is, however, representation-dependent Jl2| , so care should be 
used when applying it. 

The charge-conjugate 1p R L of a field of definite chirality ipR,L will be composed of the opposite- 
chirality field but transform like the original field |65) . For example, 

^ = Cjf (423) 



is formed from ipn but transforms as a left-handed field since C commutes with the projection 
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operators. In the standard representation, 



p n 1 ± 75 . 2 o -2 1= F 75 • 2 o 1 ± 75 nv 
Pr,lC = - i-i 7' = 17 - 7 U = 17 7 1 — — = CP RtL , 



(424) 



so the charge-conjugate field fy L has the chirality of the Dirac-conjugate field ipR. The latter field 



is shown to be left-handed in equation (415), so the field ip L transforms as a left-handed field, as 
claimed above. 



Using the properties of C and the definition (|419|) , we can find an expression for ^ L R 



= W-yV = -^. (425) 



Applied to the neutrino, we find 



-VL 



(426) 
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APPENDIX C 



A MORE REALISTIC TREATMENT OF NEUTRINO OSCILLATIONS 



The common treatment of neutrino oscillation makes many sometimes contradictory assump- 
tions, yet manages to arrive at an equation for the transition probability, (|93|), that is consistent 
with more realistic treatments to leading order in neutrino masses. These assumptions include (a) 
the treatment of neutrinos as plane waves rather than wave packets, (b) all neutrino 
mass states have a common momentum but different energy (energy is conserved for the 
interaction, but momentum conservation is ignored), (c) neutrinos are relativistic, so t fa |x|, and 

2 

(d) again the neutrino is relativistic, so E{ — p ~ In addition, we define the x-axis to be in 
the direction of the neutrino three- momentum p^, so p^ • x, — PiXi. This final relationship is merely 
a definition of the coordinate system, so we do not lose any information making it. Under these 



assumptions, the phase difference of Chapter III 



becomes 

^historic = _ [{E . _ p)x _ {E . _ p)x]= >_ L X . (428) 

Many authors have addressed the nature of the assumptions (a)-(d). Most notably, the plane- wave 
approximation (a) is addressed in references fl3l| , |6^| , [|69| and |Q. Reference j7^] argues that a 
replacement of assumption (b) with (b') all neutrino mass states have a common energy 
but different momentum is far more logical for most experiments which measure the distance 
traveled by a neutrino, not the time taken to travel. When either p or E is common to all mass 



states, the wave- packet treatment yields the traditional 0^ ls *°" c f equation ( 428 ), as will be shown 



in Section C.l. 



Neither of the assumptions (b) or (b') are correct or even preferable, according to Goldman, the 
author of reference Goldman's work provides a relativistic treatment of oscillations, conserving 
both energy and momentum, and points out that (c) and (d) are inconsistent assumptions, since the 
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difference between t and x is, to first order, the same as the difference between p and E. Goldman 
allows different mass states to have both unique energies and unique momenta, but arrives at the 
common expression (|93|) to leading order in neutrino masses. We examine his work and extend it to 



an arbitrary number of neutrino flavors in Section C.2 



C.l Neutrinos as Wave-Packets 

Nussinov J68| was perhaps the first to examine neutrino oscillations using wave packets rather 
than plane waves to describe the neutrino states. Kayser extended Nussinov's work in reference 
Kayser's work laid the groundwork for many more detailed treatments, such as j7(| and [jFl]. Refer- 



ence 1 69 provides a useful summary of the wave-packet approach. The probability for the transition 



v a —* v@ for wave packets is given by reference j7(| to be 

P (x) =AJ2J2 m Vp 3 cP^" tB C^ D - F . (429) 



i 3 



A is a normalization constant, 

-l 

i k 



IK,- 2 



A =\Y.-i^r) ■ <«»> 



k 



where vj~ represents the velocity of the fcth mass state. 



Vj + Vj _ (ft) ~ (Pj) 



B = - w w (431) 



represents the deviation of the phase from the traditional phase. C and F both result from the 
time-averaging performed in reference |7Q] . 



arises because the probability to find a mass state at a detector is inversely proportional to the speed 
of the mass state. 

((EL) - (EA) 2 , . 

4a2 (vf + vj) 
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where a p is the width of the Gaussian wave packet in momentum space, insures that energy is 
conserved within the uncertainty of the wave packet |7(J . Finally, 

x 2 (Vi — Vj) 2 

D=-^— 2 434 

Aa 2 vf+v] 

where a x is the width of the Gaussian wave packet in coordinate space, is a damping factor which 
accounts for the separation of the mass states as they propagate. As they separate, the overlap 
decreases and the oscillatory behavior disappears. This behavior gives rise to the definition of a 
coherence length, 

/ v 2 + v 2 

L^ h ^4a x J-^ (435) 

V K - vj) 

As long as the coherence length is much greater than the size of the wave packet, the different mass 
states will interfere, and oscillations will occur j7Q] . 

In the traditional approach, AC = B w 1, and D + F w 0. 



C.2 Neutrinos in a Relativistically Correct Light 



Everything in Section 3.1 is relativistically correct, so we may start with equations ( |85| ) and (|8 

v L (t = 0, x = 0) = v aL = ^2 V al v iL , 



and 



—i<f>i 



i=l 



The probability for transition P (x) , as shown in Section 3.2, depends on the phase difference 

u a -> up 



®ij = 2 _ E 3 1 j) - - Pj x i)} ■ ( 436 ) 

To find this phase in terms of the neutrino masses we must consider relativistic kinematics. 

The pi and Ei are of course related by E 2 — pf + mf , and they depend on the other energies 
and momentum of the process which created the neutrino. For simplicity, we will work in the rest 
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frame of the parent particle. The invariant mass squared of the parent particle is M 2 , and the decay 
products other than the neutrino will have total mass squared M 2 and momenta summing to — pj. 
For oscillations to occur, we must not be able to measure the momenta or energies of the other 
particles closely enough to distinguish between the different possibilities of pi . We also must not be 
able to resolve Atjj = U — tj between the arrival times of the different mass states. 
Energy conservation for the production of the neutrino requires that 

M = (E„ + E other ) = sjpj+rf + yJtf + Mj, (437) 



so 



(m - sjfi + m^j = M 2 - 2My/pf + m* + p\ + m? = pj + M 2 , 



(438) 



and 



M 2 -M 2 +m 2 

E = + = 2M~^' 



(439) 



We may extract pi from this equation: 



[M 2 - M 2 ) 2 + 2m? (M 2 - M 2 ) + m\ 4M 2 m 2 



AM 2 



AM 2 ' 



so 



(m 2 - M 2 ) - 2m 2 M 2 - 2m 2 M 2 + mf 



4M 2 



(M 2 - Mf) 2m? [M 2 + M 2 



2M 



^ (M 2 - Mf) (m 2 - Mf) ' 



(440) 



Using the Taylor expansion y/l + x = 1 + § — ^- + • • • and keeping terms to fourth order in the 
neutrino mass m i7 we find 



Pi = 



[M 2 - M) 
2M 



' M 2 - Mj 
2M 



1 2mj (m 2 + Mf) +mf m 4 4 (m 2 + M 



+ 0(m?) 



(m 2 -M 2 ) 8 (m 2 -M 2 

m\ ( M 2 + M 2 ) a ( M 2 — M 2 ) 2 — ( M 2 + M 2 Y 
[M 2 - M 2 ) Z (M 2 - M 2 ) 
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M 2 



Mj 



2M 



1 - mi 



M 2 



M 2 



2M 2 M] 



2 ~~ m i 

(m 2 - M 2 ) (M* Mf) 



■ OK 6 ) 



(441) 



We will later want the differences between two energies and between two momenta: 



Ei — Ej 



Pr ~ Pj 



M 2 



M 2 



i 2 M 2 - M 2 



2 9 

mi — mn 

l j 



M 2 



2M 
-Mj 



2M 



2M 
mi (M 2 



2M 



and 



(A/ 2 + Mj) 



4M 2 M 2 



M 2 - M 2 



(m 2 



M 2 



m 2 [M 2 + M 2 



AM 2 M 2 



(m 2 -m 2 ) 2 ( 



M 2 - Mj 



M 



Mf 



2M 



(m 2 - mi) 



M 2 + M'j 
(m 2 - Mj 



(442) 



(443) 



(m 2 — m 2 ) (m 2 + m 2 ) 



2M 2 Mj 
(M*-Mj) 



<f>ij does not, however, depend only on E and p] it also depends on x and t, and deciding what 
values to use for these variables is rather tricky. As pointed out in reference |n|, experiments 
measure the distance a; between the neutrino source and the detector, so all mass states must have 
traveled the same distance One could then find the individual times from the relativistic 

relation ti = ^-x and arrive at the phase 



J 2 



Ei p, 

Pi 



Pi 



mi 
Pi 



2p 3 



(444) 



To leading order in mass-squared, this phase is twice the phase obtained by the traditional approach: 



m, — m. 

"Pave 



(445) 



Goldman reproduces the historic phase to leading order in reference |32| by using the average 
time of travel, 



iave 



pj 

,2 







mf / 













2p 2 2p 2 
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4M 2 



4M 2 



(m 2 -M 2 Y 4 (hP-Nlfj 



Q(m\ 



(446) 



for all mass states rather than the individual travel times This choice may be justified on the 
grounds of uncertainty since we cannot resolve between the different travel times if oscillations are 
to occur. A combined treatment of wave packets and correct relativistic kinematics would of course 
illustrate the correct phase, but that is a topic for a different paper. In this paper, Gentle Reader, 
we will stick with Goldman's t ave . 



Now we may calculate $jj using equations (442), (443), and (446): 



{(Ei - E 3 )t ave - fa ~Pj): 



(E t - 



2 2 

mf — rrij 



4M 



(mf + mf) M 2 
I (AP-Af 2 ) 2 

^ M 2 + Mj (mf + rrij) M 2 2M 2 M] (mf + mf 
(m 2 -A/ 2 ) 2 {m 2 -M 2 



M 2 - M 2 



2 2 

mf — m^ 



4M 



2M 2 



M 2 - M 2 f 



2 2 

m~ — m- 

1 j 



2M 



M 2 - M 2 



(Pi-Pj) 



mf + m 2 
* j_ 

(m 2 - M 2 



[M 2 (M 2 - Mf) + 2M 2 M]] 



(mf + mf) M (M 2 + M] 



M 2 - M 2 



(447) 



The first term in the brackets, M 2I f M 'i is just the reciprocal of the leading term in all pi, so we 



2 2 

mf — m~ ... 

1 J_ ^historic 

X — 



4p 



(448) 



for small rrii, as promised above. 
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L2[, and p3| for a more detailed study of neutrinoless double-beta decay. 
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APPENDIX E 



A GRAPHICAL REPRESENTATION OF BOXES 



Many of the relationships between boxes developed in Chapter [V were originally derived using 
a graphical representation. In this method, boxes in the numerator of a product are represented by 
two vertical lines; boxes in the denominator are represented by two horizontal lines. Lines exit the 
locations of the matrix elements which are not complex conjugated and enter the locations of the 
complex-conjugated matrix elements. For example, the box n D22 = Vu ^2^22^2* is represented 
by a vertical line pointing from V\\ to V21 and a vertical line pointing from V22 to V12, as shown in 
Figure [l4|a. The inverse box — = ( 11 CI22 ) 1 is represented by a horizontal line pointing from 



V\\ to V12 and one pointing from V22 to V21, as shown in Figure 14 3. The complex-conjugated box 
11 ^22 is equal to 12 ^21 , so one just reverses the arrows to complex conjugate a box, as shown in 



Figure 14 



II I J 

• • • *-* • • • • • 

• •• ••• ••• 

(a) (b) (c) 

Figure 14: The graphical representation for (a) 11 CI22 , (b) ( 11 ^22 )~ , and (c) lx n| 2 ■ 

To multiply boxes together graphically, one merely draws the lines corresponding to each factor 
on the same grid. Horizontal arrows entering or leaving a point cancel out vertical arrows enter- 
ing or leaving, respectively, that point. Uncanceled arrows entering a point signify the survival 
of the complex-conjugated matrix element associated with that point. Those leaving a point sig- 
nify the survival of the ordinary matrix element of that point. Figure |lqa represents the product 



n n 22 12 D 33 ( 23 D 32 ) . The horizontal arrows do not cancel vertical arrows at any point, so no 
simplification may occur. The upward arrow at V22 represents that element in the numerator. The 
horizontal arrow represents the element V 2 \ in the denominator, which does not cancel. Counting 
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off the arrows at each vertex, we find the expression 



ll n 12 D / 23 n \-! = Vl 1^12^2^1*3^21^22^*2^33 
22 33 1 32 ' ^1/23^32^33 



(450) 



which agrees with the definitions of boxes. 



• 



(a) 



(b) 



Figure 15: The graphical representation for the products a) ^22 Q 

b) 11 ^ 12 n 33 ( 22 n 3 3) _1 . 



12 n ( 23 n 
J 33 { u 32 



and 



Figure |T^b represents n D22 12 ^33 ( 22 ^33 ) > a product in which some canceling does occur. 
Picking out the uncanceled arrows at each vertex, we are left with 



"□22 12 n 



33 



( 22 n 3 3) = 



VnVizV&V&VZ 



(451) 



23 



The graphical method is a powerful tool for finding relationships between boxes. Consider the 
degenerate boxes m O ai = |V Q i| 4 . Trying to obtain the equations (143) and (144) presented in 



Chapter IV without graphs required quite a bit of running down dead ends. Using the graphical 
method, we need only find a series of arrows which cancel for every point except (a, i) and leave 
two incoming and two outgoing vertical arrows at that point. Consider IV21I as an example. We 
choose to draw all of the arrows involving V21 pointing downward, as shown in Figure |l6|a. These 
arrows must be part of boxes, so in Figure |l6| b we add the arrows to finish those boxes. Next we 
draw two horizontal boxes in Figure |l6| c to cancel the extra arrows in the first column of the matrix. 
This still leaves V22 and V23 with two sets of uncanceled arrows apiece. In Figure [l^ d we draw two 
more horizontal boxes to compensate. This adds arrows to our previously clean V12, V13, V32, and 
V33; drawing the final vertical box in Figure |l6|e cancels those. Recapping what we have done, we 



see that step (b) completes ^22, n 



23 



21r 



□32 , and O33 in the numerator. Step (c) divides 
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(d) 



(e) 



Figure 16: The steps to obtaining | V21 1 4 as a function of ordered, non-degenerate boxes. The 
additions in each step are designated by the thicker arrows. 

by 11 D 3 2 and 11 n 33 , and step (d) divides by 12 D 2 3 and 22 33 . Step (e) multiplies by 12 D 33 , 
leaving only the point V21 with uncanceled arrows. It has two vertical arrows coming in and two 
leaving, so our graph represents the equation 



V21I = 



2i n 11 n 21 n Hn 12 n 

u 32 u 22 u 33 u 23 u 33 

n n 3 2 n n 33 22 n 33 12 n 23 



(14£) 



of Chapter IV. Other examples of this representation are included in that Chapter. 
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